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ABSTRACT 


f 

Experiments  were  carried  out,  over  a  period  of  two  years, 
in  an  Adaptable-Wall  wind  tunnel  configured  for  testing  of 
airplane  models  at  very  large  lift.  The  program  was  intended 
especially  to  demonstrate  Proof  of  Concept  for  this  type  of 
wind  tunnel,  in  which  the  simulated  stream  vector  is  inclined 
appreciably  to  the  tunnel  axis.  The  measured  inner  flow  is 
matched  to  the  computed,  updated  outer  flow  by  an  iterative 
process. 

Wal 1 -adaptati on  controls  in  this  tunnel  are  vaned  panels 
in  the  floor  and  ceiling  of  the  working  section  and  a 
variable-angle  inlet  nozzle.  Velocity  components  are  measured 
by  a  Laser-Doppl er  system  using  a  fixed  laser  and  movable 
optical  components. 

The  test  model  used  in  these  experiments  was  a  high-wing 
V/STOL  configuration  having  full-span  wing  flaps  with  lower- 
surface  blowing  of  their  inboard  portions.  In  all  of  the 
experiments  reported  here,  the  combination  of  angle  of  attack, 
flap  setting,  and  flap  blowing  was  such  as  to  produce  large 
flow  deflection  and  severe  wall  interference  in  a  conventional 
tunnel . 

The  model  conf i gurati on  was  always  laterally  symmetrical, 
and  most  runs  were  made  under  the  assumption  of  symmetrical 
flow. 


The  iterative  procedure,  based  on  measured  control  mat¬ 
rices,  typically  led  to  minimum  matchi ng-di screpanci es  (root- 
mean-square  values)  of  about  three  percent  of  stream  speed 
after  about  six  iterations.  It  is  estimated  that  this  reflects 
residual  errors  at  the  model  of  about  one  percent  of  stream 
speed . 

It  is  concluded  that  these  results  constitute  successful 
Proof  of  Concept.  Suggestions  are  made  regarding  the  direc¬ 
tions  of  further  development  of  this  type  of  wind  tunnel. 


BACKGROUND 


The  Adaptable-Wall  Wind-Tunnel  idea  dates  -from  the 
early  1970’s  when  it  was  independently  invented  by  Ferri 
(Re-ference  1)  and  Sears  (Reference  2).  It  was  an  invention 
intended  to  solve  a  persistent  problem  of  wind-tunnel  tes¬ 
ting,  namely  boundary  interference.  The  effects  of  bounda¬ 
ries  can  be  corrected  for,  in  many  wind-tunnel  tests  at  low 
or  moderate  speeds,  but  this  cannot  be  done  in  nonlinear 
regimes,  such  as  the  transonic,  especially  in  ventilated 
wind  tunnels. 

In  briefest  outline,  the  Adaptable-Wall  scheme  consists 
of  modifying  the  tunnel’s  geometry  so  as  to  produce  the  cor¬ 
rect,  unconfined  flow  within  the  working  section,  in  the 
presence  of  the  model,  by  adjustment  of  "wall  controls". 

This  is  accomplished  by  measuring  two  independent  flow- 
variable  distributions  at  an  interface  that  encloses  the 
model,  calculating  the  outer  flow  field  exterior  to  this 
interface  (based  on  the  measured  values),  and  progressively 
reducing  the  matching-discrepancy  at  the  interface  by  actua¬ 
tion  of  the  wall  controls. 

It  is  clear  that  when  the  experimental  (inner)  and 
calculated  (outer)  fields  are  matched  at  the  interface  the 
whole  combined  field  is  correct,  for  the  outer  flow  is 
always  made  to  satisfy  the  far-field  boundary  conditions; 


own  boundary  conditions  in  the  tunnel. 

Thus,  an  Adaptable-Wall  Tunnel  needs  three  special 
features:  (1)  wall  controls,  (2)  instrumentation  to  measure 
two  distributions  at  the  interface,  and  <3)  computer  capabi¬ 
lity  to  calculate  the  outer  flow  for  given  boundary  values 
at  the  interface.  Several  such  tunnels  have  been  construc¬ 
ted  in  various  countries. 

In  the  process  of  pursuing  this  development  in  the 
1970s,  the  senior  author  of  this  report  became  aware  that 
the  Adaptable-Wall  scheme  also  replaces  the  customary  "cali¬ 
bration"  of  a  wind  tunnel;  the  determination  of  the  stream 
vector  (i.e.,  the  negative  of  the  flight-velocity  vector)  is 
made  by  choosing  this  vector  in  the  outer-flow  calculation 
and  carrying  out  the  elimination  of  the  matching-discre¬ 
pancy.  This  may  very  well  turn  out  to  be  the  most  important 
feature  of  the  Adaptive-Wall  Tunnel,  since  the  traditional 
"calibration",  which  is  carried  out  with  an  empty  tunnel  and 
then  used  for  tests  of  models,  is  notoriously  unreliable  in 
both  principle  and  practice  —  a  very  dubious  procedure. 

Besides  the  transonic  flight  regime,  another  nonlinear 
regime  in  which  difficulties  in  wind-tunnel  testing  are 
encountered  is  that  of  very  large  lift  coefficients,  as 
typified  by  the  low-speed  V/STOL  regime  and  powered  lift. 
Here  it  is  principally  the  f 1 oor  and  ceiling  of  the  working 
section  that  interfere  with  the  highly  deflected  airstream 
and,  especially,  the  powered,  vortical  wake.  At  the  Univer- 


sity  of  Arizona,  a  program  of  research  was  undertaken  to  see 
whether  this  persistent  problem  of  the  low-speed,  high-lift 
regime  could  be  solved  by  exploitation  of  the  Adaptable-Wall 
strategy . 


The  first  three  years  of  this  program  were  devoted  to 
numerical  simulations  of  wind-tunnel  tests  of  high-lift 
configurations.  The  first  breakthrough  came  when  it  was 
shown  that  the  interface  could  be  a  five-sided  box  with  open 
downstream  end  and  with  the  vortical  wake  extending  down¬ 
stream  through  the  open  end.  The  top,  bottom,  and  sides  of 
the  interface  were  modelled  as  di str i buted-vortex  panels, 
which  have  trailing  vortices  running  back  from  their  lateral 
edges,  so  that  the  interface  is  essentially  semi -i nf i ni te  in 
length.  This  made  it  unnecessary  to  try  to  locate  the  wake 
and  to  model  it  in  the  outer-flow  calculation  —  a  major 
simplification. 


Next,  it  was  ascertained  that  the  free-stream  direction 
could  be  chosen  to  make  a  large  angle  to  the  nominal  "tunnel 
axis"  —  such  as  40  degrees  —  and  the  scheme  of  iteration 
successfully  reduced  the  matching-discrepancy  at  the  inter¬ 
face  to  negligible  values.  The  results  of  these  studies 
were  published  in  References  3  and  4,  where  a  radical  new 
type  of  wind  tunnel  was  proposed,  which  seemed  to  offer  a 
solution  to  the  high-lift  testing  problem:  It  would  be  an 
Adaptable-Wall  Tunnel  in  which  the  free-stream  vector  would 


would  be  controlled  by  the  Adaptable-Wall  principle;  the 
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test  model  would  be  mounted  nose-down  to  give  the  desired 
angle  of  attack;  and,  consequently,  the  powered  vortical 
wake  would  trail  harmlessly  downstream,  essentially  horizon¬ 
tally. 

I 
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The  resulting  wind  tunnel  was  called  the  "Arizona  Wind 
Tunnel"  and  was  proposed  as  a  solution  of  the  l  lgh-lift 


i 

i 

testing  problem  described  above,  characteristic  of  the 
V/STOL  class  of  aircraft.  The  present  program  is  an  effort 
to  demonstrate,  in  a  laboratory  embodiment,  that  the  concept 
is  sound. 


The  principle  of  the  Arizona  Tunnel,  which  i s  to 
divorce  the  flight  direction  from  the  architecture  of  the 
laboratory  and  the  wind  tunnel,  would  also  seem  to  have 
applications  other  than  testing  of  high-lift  aircraft  confi¬ 
gurations.  There  are  other  categories  where  it  may  be  more 
convenient  to  change  angles  of  pitch  and  yaw  by  rotating  the 
free-stream  vector  than  by  rotating  the  model.  For  example, 
there  are  many  tests  where  models  are  connected  to  fuel 
lines,  exhaust  collectors,  etc.,  and  are  therefore  difficult 
to  move.  There  is  also  the  case  of  vehicle  models  being 
tested  in  the  presence  of  a  moving  ground  plane  representing 
a  roadway  or  runway;  it  may  be  important  to  simulate  cross- 
winds  in  such  tests.  Nevertheless,  the  investigations 
reported  here  are  confined  to  the  testing,  at  zero  yaw,  of  a 
typical  high-lift  airplane  model  having  blown  flaps  (de¬ 
scribed  later  in  the  report). 


OBJECTIVES 


Since  the  major  goal  of  this  research  has  been  Proof  of 
Concept,  no  effort  was  made  to  produce  a  facility  suitable 
for  routine  testing.  Neither  the  instrumentation  nor  the 
wall  controls  is  "on-line".  The  working  section  is  20  in. 
by  20  in.  in  cross-section.  Tunnel  speed  can  be  about  70 
knots  (maximum),  but  its  magnitude,  within  reason,  is  not 
considered  to  be  important.  Moreover,  three  important 
simplifications  have  been  introduced,  which  have  greatly 
reduced  both  cost  and  complication  of  this  project: 

1.  No  balance  system  has  been  provided.  This  decision 
was  based  on  the  conviction  that  the  basic  principles 
of  the  Adaptable-Wall  scheme  have  been  verified  in 
other  embodiments  (principally  in  the  transonic 
regime).  Thus,  success  in  iterating  to  unconf ined-f low 
conditions  is  measured  by  the  magnitude  of  the  match¬ 
ing-discrepancy  distribution,  rather  than  by  any 
attempt  to  compare  test  results  with  results  obtained 
with  the  same  model  in  a  very  large  tunnel  or  any  other 
"exact"  or  "target"  results.  Clearly,  it  would  be  a 
major  extension  of  this  work  —  an  extension  that  may 
be  desired  at  some  later  date  —  to  obtain  such  "tar¬ 
get"  results. 

2.  Since  balance  readings  are  not  available,  lift  coef¬ 
ficients  produced  in  our  tests  can  only  be  estimated. 


H 

I 

I 

E 

! 

B 


I: 

I 

I 

£ 


"a 


I 

i. 


ii 


Moreover,  operating  values  of  the  overall  "jet-flap 
coef f i cient"  for  the  model  are  not  precisely  known. 
Instead  of  accurate  values  of  these  parameters,  we 
propose  that  experimentally  observed  jet-wake  geomet¬ 
ries  be  understood  to  characterize  the  high-lift  regime 
of  the  tests.  Geometries  were  measured  and  are  presen¬ 
ted  here.  It  is  our  contention  that  this  characteriza¬ 
tion  is  more  significant  to  the  reader  in  evaluating 
the  success  of  the  tunnel  scheme  than  precise  knowledge 
of  lift  and  jet  coefficients  would  be.  Our  tunnel  is 
intended  to  permit  testing  of  1 ar ge-def 1 ecti on  flows, 
especially  those  with  powered  wakes  —  rather  than 
flows  at  any  particular  lift  coefficient  or  jet-flap 
coefficient . 

3.  The  working  section  was  provided  with  only  a  rela¬ 
tively  small  number  of  wall  controls  —  16  vaned 
panels,  of  which  the  vane  angles  can  be  varied  between 
closed  and  full  open,  a  variable-angle  nozzle  leading 
from  the  settling  chamber  to  the  working  section,  and  a 
throttling  device  downstream  of  the  working  section, 
which  controls  the  pressure  differences  across  the  top 
and  bottom  vaned  panels  —  IS  controls  in  all.  No 
controls  were  provided  on  the  side  walls.  According  to 
numerical  simulations,  this  array  of  controls  should  be 
adequate  to  reduce  the  discrepancy  function  to  satis¬ 
factorily  small  values. 

Besides  (1)  Proof  of  Concept,  major  goals  of  this  re- 
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search  have  been: 


(2)  Demonstration  o-f  Laser — Doppler  Anemometry  as  an  instru¬ 
mentation  system  -for  Adaptive-Wall  Wind  Tunnels  in  the  fully 
three-dimensional  application.  Success  in  this  demonstra¬ 
tion  confirms  the  positive  results  already  achieved  at  Ames 
Research  Center. 

<3)  Preliminary  determination  of  the  performance  of  wall- 
control  panels  consisting  of  adjustable  vanes. 

< 4 )  Preliminary  determination  of  the  attainable  accuracy  of 
simulation  of  unconfined  flow  about  high-lift  models  at 
selected  free-stream  vectors  in  wind  tunnels  of  the  Arizona 
type. 

(5)  Acquisition  of  preliminary  experience  in  practical  oper¬ 
ating  procedures  with  tunnels  of  the  Adaptabl e-Wal 1  type. 
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EQUIPMENT 


Wind  Junnei 

The  tunnel  used  in  these  investigations  is  shown  in 
Figure  1  and  in  photographs  reproduced  at  the  end  of  this 
report.  It  is  an  open-return  tunnel,  powered  by  a  three-phase, 
induction  motor,  rated  25  h.p.  at  1765  r.p.m.,  which  drives  a 
ventilating  blower  through  a  belt  drive.  The  blower  is  rated 
at  18200  cubic  -feet  per  minute  at  1695  r.p.m.  Nominal  tunnel 
airspeed  is  controlled  by  adjustment  o-f  an  inlet  valve  up¬ 
stream  o-f  the  blower.  This  motor /blower  combination  is 
observed  to  run  accurately  at  constant  speed  (1640  blower 
r.p.m.,  plus  or  minus  4)  in  all  the  experiments  reported  here, 
regardless  o-f  inlet-valve  position,  model  con-figuration,  or 
settings  o-f  the  wall  controls. 

The  blower  output  goes  through  a  coarse  screen  (0.5  x 
0.5-inch  mesh)  directly  into  a  wide-angle  diffuser  with  split¬ 
ter  plates,  then  into  the  settling  chamber.  There  are  four 
screens  in  the  settling  chamber,  as  shown  in  Figure  Is  these 
are  plastic  window  screening,  14-mesh,  seamless. 

As  seen  in  Figure  1,  the  nozzle  leading  air  from  the 
settling  chamber  into  the  working  section  is  adjustable  in 
angle:  this  is  one  of  the  " wal 1 -control "  organs.  Through  an 
extension  of  the  settling  chamber,  air  is  also  admitted  to  the 


working  section  through  control  panels  in  its  -floor.  There  are 


eight  panels  in  the  floor  and  eight  in  the  ceiling;  these  are 
rectangular  panels  made  up  of  vanes  of  0.5-inch  chord.  The 
vane  angle  of  each  panel  is  variable  in  angle  from  fully 
closed  to  about  75  degrees  open.  At  the  downstream  end  of  the 
working  section,  as  shown  in  Figure  1,  there  is  an  eighteenth 
control  organ,  namely  a  valve  made  up  of  flat  aluminum-alloy 
plates  of  3-inch  chord,  which  are  turned  in  alternating  direc¬ 
tions  so  as  to  throttle  the  airflow. 


Each  of  the  vaned  panels,  as  well  as  the  downstream 
valve,  is  controlled  by  a  vernier  dial/rotator  that  drives  the 
vane  angle  through  a  worm  drive.  These  are  accurately  control¬ 
lable  to  approximately  1/10  degree. 


m 


In  the  first  runs  made,  total-head  surveys  were  made  in 
the  empty  working  section  with  controls  set  to  produce  30-  to 
40-degree  flow  inclination.  Not  surpr i si ngl y ,  appreciable 
differences  were  found  between  the  total  pressures  of  air 
coming  through  the  nozzle  and  through  the  floor  panels.  Trials 
were  then  made  with  various  combinations  of  screening  over  the 
nozzle  entrance  and  over  the  entrance  to  the  settling-chamber 
extension  that  feeds  the  floor  panels.  A  combination  was  found 
that  minimized  the  total -pressure  difference  and  also  made  the 


flow  in  the  working  section  more  steady.  These  screens  were 
permanently  installed. 


A  total -head  survey  was  carried  out  over  a  transverse 
plane  just  upstream  of  the  test  model,  with  wall -control 


settings  determined  for  a  high-lift  model  conf  igurati  on  (viz.. 
Experiment  91).  Total  pressure  was  found  to  be  uniform  to 
within  one  percent. 

Instrumentation 

Instrumentation  consists  primarily  of  a  Laser-Doppl er- 
Anemometer  system  (L.D.A.)  based  upon  an  argon-ion  laser, 
nominally  2-watt,  manufactured  by  Lexell,  TSI  optical  compo¬ 
nents  and  counter,  and  a  traversing  frame  that  moves  the 
focussing  lens  and  appropriate  mirrors  in  three  mutually  per¬ 
pendicular  directions. 

This  system  is  shown  in  Figure  2. 

Being  a  key  element  of  this  research  project,  the  L.D.A. 
system  naturally  occupied  much  of  our  time  and  attention  as 
the  project  was  planned  and  put  into  operation.  It  could 
occupy  much  of  this  report,  but  we  will  confine  it  to  a  rather 
brief  account  and  some  conclusions. 

It  is  a  single-component  back-scatter  system.  The  laser 
is  fixed  in  position,  together  with  the  collimator,  polariza¬ 
tion  rotator,  beam-steering  module,  beam  splitter,  signal 
receptor,  and  photo-multiplier.  As  shown  in  the  diagram,  the 
rest  of  the  optical  system,  namely  three  mirrors  and  the 
focussing  lens  are  mounted  on  the  traversing  frame.  Thus  the 
measuring  volume  is  translated  in  three  mutually  perpendicular 
directions,  say  x,  y,  and  z. 

As  has  already  been  mentioned,  the  adaptab 1 e-wal 1  scheme 


requires  measurement  of  two  velocity  components.  In  our  exper¬ 
iments  this  is  accomplished  by  manually  rotating  the  beam¬ 


splitter  assembly.  In  these  experiments,  involving  as  they  do 
large  inclinations  of  the  undisturbed-flow  vector,  it  has  been 
found  that  the  recognized  ambiguity  of  L.D.A.s  near  zero 
velocity  does  not  pose  a  problem.  (The  ambiguity  arises 
because  the  system  times  the  passage  of  seed  particles  through 
the  measure  volume  and  does  not  indicate  their  direction, 
unless  special  equipment  is  provided.)  When  the  direction  of  a 
component  is  ambiguous,  as  in  a  few  runs  made  near  zero  flow 
inclination,  components  at  plus  and  minus  45  degrees  were 
measured,  from  which  horizontal  and  vertical  components  were 
computed;  in  all  other  cases  these  were  measured  directly. 

The  process  of  aligning  the  optical  components  of  the 
L.D.A.,  so  that  a  good  signal  was  obtained  at  all  the  points 
of  the  working  section,  was  tedious.  It  was  facilitated  by 
removing  the  back  (West)  wall  of  the  working  section  and 
placing  a  large  vertical  board  there,  upon  which  the  two  laser 
beams  impinged  and  could  be  monitored  —  translated  in  the  x, 
v,  and  z  directions  and  rotated  about  the  optical  axis.  The 
assistance  and  advice  of  Dr.  Ari  Glezer  was  invaluable. 

After  some  experimentation,  it  was  found  that  a  liquid 
seed  consisting  of  80-percent  water  and  20-percent  glycerine, 
introduced  into  the  flow  at  the  blower  intake,  produced  a 
satisfactory  back-scatter  signal  with  indicated  laser  output 
of  1  watt.  It  will  be  appreciated  that  ours  is  an  open-return 
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tunnel,  so  that  there  were  severe  limitations  on  the  seed  that 
could  be  used.  Pure  water  was  tried,  but  was  -found  to  evapor¬ 
ate  before  reaching  the  working  section.  The  water-glycerine 
mixture  used  makes  no  noticeable  deposit  in  the  laboratory, 
even  when  the  tunnel  is  operated  all  day.  Exposed  optical 
components,  namely  mirrors,  need  to  be  cleaned  after  about  two 
days  of  operation.  The  plate-glass  East  wall  of  the  working 
section  is  the  exception:  it  is  found  advisable  to  clean  its 
inner  side  after  every  two  hours  of  tunnel  running. 

From  the  photo-multiplier,  our  back-scatter ed  signal  goes 
into  a  counter  manufactured  by  TSI  Inc.;  thence,  via  an  inter¬ 
face  designed  and  constructed  for  this  project  by  H.E.  (Dutch) 
Hal  deman  of  our  department,  into  our  Osborne  1  micro-compu¬ 
ter  . 

Other  instrumentation  includes  a  counter  that  monitors 
blower  r.p.m.  and  a  pressure  transducer  (Setra) .  with  asso¬ 
ciated  digital  voltmeter,  that  monitors  the  static-pressure 
difference  between  settling  chamber  and  laboratory. 

Jest  Model 

Except  for  some  empty-tunnel  runs,  described  later,  all 
runs  reported  here  were  made  with  our  "generic"  STOL  airplane 
model  installed.  This  model,  sketched  in  Figure  3,  represents 
a  high-wing  transport  airplane  with  lower-surf ace-bl own  wing 
flaps.  Although  its  flaps  are  full-span,  they  are  blown  by 
flat  nozzles  mounted  on  the  wing  lower  surface,  which  directly 


affect  only  the  inboard  part  of  the  wing.  Flap  blowing  air  is 


piped  to  the  model  -from  the  laboratory's  compressed-air  system 
("shop  air")  through  the  model -support  strut.  Pressure  gauges 
are  installed  in  this  air  line  to  monitor  both  the  supply 
pressure  and  the  pressure  drop  across  an  ori-fice  in  the  line. 

Tuft  surveys  showed  clearly  that,  at  high  angles  of 
attack,  with  flaps  deflected  to  60  degrees,  the  wing  is  com¬ 
pletely  stalled,  without  flap  air.  When  flap  air  is  provided, 
the  wing  and  flaps  are  unstalled  inboard  —  at  least  60  per¬ 
cent  of  the  half-span  —  but  still  violently  stalled  outboard. 


PROCEDURES 


Measurement  of  Vel.gci.ty  Corngonen ts 

The  -flow  in  the  working  section  of  our  demonstration 
tunnel  is  turbulent  and,  with  the  model  at  the  high-lift 
conditions  described  above,  grossly  deflected  (as  will  be 
shown)  and  unsteady.  Consequently,  all  readings  of  velocity 
components  must  be  mean  values;  we  presume  that  this  is  true 
of  all  wind-tunnel  tests  in  the  high-lift,  partially  stalled 
regime.  Moreover,  we  also  think  that  the  flow  conditions 
observed  in  our  tunnel  are  probably  character! sti c  of  flight 
in  this  regime  and  that,  therefore,  the  measurement  of  time- 
averaged  quantities  is  necessary  to  obtain  meaningful  data. 


Our  test  procedures  have  therefore  been  developed  to 
measure  mean  values,  averaged  over  a  long  enough  time-span  to 
represent  statistically  steady  values.  At  the  beginning  of  the 
tests,  after  some  trials,  we  decided  that  an  average  over  100 
"samples"  was  adequate  to  meet  this  requirement.  The  TSI 
counter  is  designed  to  accept  only  those  signals  (a  "signal" 
being  the  trace  of  a  particle  passing  through  the  interference 
bands  at  the  field  point)  that  meet  minimum  specifications 
regarding  number  of  full  cycles  and  uniformity  of  cycle  fre¬ 
quency.  In  our  tests  we  set  these  requirements  at  eight  cycles 
and  not  more  than  seven  percent  difference  between  successive 
cycles  within  the  signal.  The  word  "sample"  is  used  here  to 


mean  a  signal  accepted  by  the  counter. 


It  required,  typically,  about  45  seconds  to  obtain  such 
an  average  over  100  samples. 

As  the  research  progressed,  we  became  aware  that  the  -flow 
in  the  working  section  was  much  steadier  at  some  -field  points 
than  at  others.  In  particular,  the  flow  at  points  closer  to 
the  wing  wake  was  more  unsteady.  Accordingly,  we  devised  a 
procedure  of  extended  averaging  that  would  take  more  samples 
at  points  where  the  flow  was  more  unsteady. 

The  following  procedure  was  arrived  at: 

(1)  A  "reading"  consists  of  the  arithmetic  mean  value  of 
a  number,  NSPL,  of  samples  accepted  by  the  TSI  counter. 
Typically,  NSPL  was  put  equal  to  50. 

(2)  The  computer  was  programmed  to  repeat  the  process  of 
taking  readings  and  to  calculate  the  running  average  of 
the  readings  taken.  This  process  was  repeated  until  four 
successive  averages  fell  within  a  specified  tolerance, 
TOL.  The  last  average  value  was  then  accepted  and  recor¬ 
ded  as  the  value  of  the  velocity  component.  Typically, 

TOL  was  set  at  less  than  one  percent  of  the  nominal 
resultant  stream  speed. 

It  will  be  seen  that  the  result  of  this  procedure  is  that 
each  recorded  value  is  the  mean  of  at  least  200  samples,  and 
that  the  number  of  samples  and  the  time  averaged  over  is 
larger  at  points  where  the  flow  is  unsteady  —  i.e.,  where  an 
increased  number  of  readings  is  required  to  bring  four  succes- 
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sive  average  values  within  the  specified  limit  TOL.  In  our 


experience  to  date,  the  greatest  number  of  readings  required 
by  this  procedure  was  ten;  i.e.,  in  this  case  the  recorded 
experimental  value  was  the  average  of  500  samples. 

Typically,  15  to  20  seconds  are  required  to  obtain  a 
reading  with  NSPL  equal  to  50;  hence  this  procedure  of  exten¬ 
ded  averaging  consumes  60  to  80  seconds  at  field  points  where 
the  flow  is  relatively  smooth  and  twice  as  much  time  at  our 
"rough  points"  —  at  most,  in  a  few  cases,  about  three 
minutes . 


We  are  sure  that  the  accuracy  of  our  test  data  was  im¬ 
proved  when  we  introduced  the  extended  averaging  procedure, 
but  that  the  data  obtained  previously,  using  simply  an  average 
over  100  samples,  are  also  acceptable. 


The  Adaptabie^Wal 1  Algorithm 

The  logic  of  the  adaptable-wall  wind-tunnel  scheme  has 
been  set  forth  under  BACKGROUND,  above.  It  has  been  converted 
into  a  detailed  procedure  and  presented  in  several  papers 
since  it  was  first  proposed  in  References  1  and  2.  (See,  for 
example.  References  3  and  4).  The  version  used  here  follows 
Reference  4: 

Let  f  and  g  be  two  independent  velocity-component 
distributions  at  the  inner/outer-f low  interface,  and  let 
the  subscript  m  denote  measured  values.  Now  let  the 
notation  fCg]  denote  that  a  component  distribution  f,  at 
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the  inter-face,  is  calculated  for  the  outer  field  by  using 
the  component  g,  at  the  interface,  as  inner  boundary  data 
and  satisfying  the  far-field  boundary  condition  (uniform 
flow  at  a  prescribed  angle,  in  our  case).  For  given  f  and 
g,  the  "mismatch"  distribution  at  the  interface  is 

f  Cgm3  -  fm  =  Df ,  sav.  (1) 

The  adaptable-wall  strategy  is  to  adjust  the  tunnel 
walls  and  tunnel  speed,  so  as  to  add  a  fraction  k*Df  to 
the  measured  values  fm  and  to  repeat  the  process  ltera- 
tivelv  until  the  mismatch  is  reduced  to  acceptable  magni¬ 
tude. 


The  factor  k  is  called  a  "relaxation  factor",  and 
experience  shows  that  it  must  be  less  than  one  to  avoid 
overshoot . 

Introducing  a  superscript  notation,  namely  letting 
superscript  (p)  denote  values  measured  in  the  pth  itera¬ 
tion  (p  =  1,2,...),  we  have 

d(p)  f  =  fCgaCp)]  -  fm(p)  (2) 

and  f  m  (p+1  )  is  adjusted,  as  closely  as  possible,  to 

fm(p+1)  =  fm(p)  +  k*D(p)f  <3> 

The  procedure  used  to  set  the  wall  controls  as 
required  by  Equation  (3)  will  be  explained  later. 

It)1?  HEigyE®  of  Merit"^  RMS 

In  actual  experiments  it  will  always  be  impossible  to 
drive  the  matching  d i scr epanc l es ,  or  mismatch  values,  at  the 
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interface  to  zero,  because  the  number  o-f  control  organs  is 
finite,  the  modelling  of  the  outer  flow  is  approximate,  and 
there  are  other  sources  of  experimental  error.  The  operator  of 
an  Adaptable  Wind  Tunnel  must  therefore  decide  upon  a  “Figure 
of  Merit"  that  is  a  measure  of  how  good  a  matching  has  been 
achieved,  overall,  at  the  interface,  and  must  undertake  to 
reduce  the  Figure  of  Merit  to  a  minimum. 


Several  suggestions  have  been  made  of  possible  choices 
for  this  Figure  of  Merit,  such  as  the  mean  matching  error,  the 
mean  absolute  matching  error,  various  weighted  averages  of  the 
matching  error,  etc.  In  this  investigation  we  have  chosen  the 
r oot-mean-sguare  matching  error;  viz.,  in  the  notation  used 
above , 

RMS  =  root-mean-sguare  value  of  Df 

N 

=  SOR<  (1/N)  *  'ST  (Dfi>2)  (4) 

where  N  is  the  number  of  field  points. 

?!  Calculation  of  Outer  Flow 

Our  interface  is  similar  to  the  one  defined  in  Reference 
4,  namelv  a  rectangular  box,  semi -i nf i ni te  in  length.  (See 
Figure  4  .  The  coordinate  system,  x,y,z,  is  also  shown.) 

To  calculate  the  outer  flow  field,  singularity-panels  are 
arraved  on  the  interface.  These  are  di str l buted-vortex  panels 
on  the  top,  bottom,  and  sides  of  the  interface  and  source 
osnels  on  the  front.  The  use  of  singularity-panels  requires 
that  the  velocity  components  f  and  g  be  redefined  as 


L/V  vs  v  v  -r  -  .  -  .  .  •  .  -r. 


19 


perturbat i gns  measured  -from  the  far— field  values  of  the 
appropriate  velocity  components.  This  does  not  involve  any 
changes  in  Equations  (1)  to  (4);  it  means  that  the  far-field 
boundary  values  of  f  and  g  are  zeros. 

These  perturbation  velocity  components,  for  the  experi¬ 
ments  reported  here,  have  been  defined  as  follows: 

q  is,  in  all  cases,  the  tangential  perturbation  at 
the  center  of  the  respective  panel;  thus,  at  the  four 
front  panels  g  is  the  vertical  (z>  component,  and  at  all 
other  panels  it  is  the  x  component.  This  component  is 
usually  referred  to  as  "Vfc"  or,  in  some  computer  pro¬ 
grams  ,  " Vj "  . 

f  is,  in  all  cases,  a  combination  of  normal  and 
tangential  perturbations.  At  the  four  front  panels  it  is 
the  normal  (x)  component  at  panel  center;  at  the  16  top 
and  bottom  panels  it  is  the  normal  (z)  component  at  panel 
center;  at  the  12  side  panels  it  is  the  x  component  at  a 
point  one  inch  outboard  of  the  panel  center.  This 
component  is  referred  to  as  "Ve  "  or,  in  some  compu¬ 
ter  programs,  "V2" ■ 

In  Appendices  A  and  B  of  this  report  are  reproduced  our 
BASIC  proqrams  to  calculate  the  32-bv-32  matrices  that  relate 
the  distributions  Ve  and  to  the  32  panel  strengths;  these 
matrices  are  called  "IntTMatE"  and  " IntTMatB" ,  respecti vel v. 
The  programs  reproduced  contain,  of  course,  panel  coordinates 


and  dimensions,  as  well  as  subroutines  that  calculate  the 
velocity  -fields  of  source  panels  and  distributed-vortex 
panel s. 

For  the  outer  flow  field  defined  by  these  singularity- 
panels.  the  calculation  designated  by  fCg],  above,  consists 
merely  of  a  matrix  multiplication;  namely, 

Ve  CVt3  =  E  *  (B^y  )  *  Vt  (5) 

where  E  and  B  are  the  square  matrices  described  above,  rela¬ 
ting  V  e  and  to  the  panel  strengths,  and  B^y  is  the  inverse 
of  B.  Our  procedure,  of  course,  has  been  to  construct  the 
square  matrix  £  *  <B  j^y)  ,  as  in  Appendix.  C,  and  provide  it  on 
our  working  disks;  it  is  a  function  of  interface  geometry 
only.  The  result  is  that  the  outer-flow  calculation  becomes  a 
simple  and  rapid  step  in  any  of  our  runs. 

?9Qtrgl_  Matrix  &  Qal_cuiat ign  of  Controi  Settings 

A  crucial  step  in  the  adaptable-wall  algorithm  is  the  one 
represented  by  Equation  (3),  namely  the  process  of  altering 
the  measured  velocity  components,  by  adjustment  of  the  wall 
controls,  so  as  to  reduce  the  matching  discrepancy  at  the 
interface.  Our  procedure  for  this  step  is  to  measure  the 
matrix  of  control  effects  and  undertake  to  calculate  the  array 
of  control  increments  required  to  satisfy  Equation  (3).  Speci¬ 
fically,  since  D  (p)f  in  Equation  (2)  represents  D^p)ve,  we 
measure  the  matrix  that  relates  increments  of  Ve  to  increments 
of  control  settings. 


In  general  ,  one  must  expect  that  control  effects  are 


dependent  upon  the  existing  -flow  in  the  working  section  and 
therefore  upon  both  the  model  conf i gur ati on  and  the  existing 
control  settings.  To  employ  a  control -effect  matrix  (or,  pre¬ 
sumably,  any  other  computati onal  procedure)  at  this  point  must 
imply  that  the  increments  involved  will  be  small  enough  to 
permit  a  local  1 inearization  of  the  process.  In  this  report  we 
hope  to  present  data  that  will  cast  light  on  this  subject. 


Let  us  proceed  under  the  assumption  that  this  local 
linearization  is  possible.  An  experiment  is  then  carried  out 
in  wh.ch  each  control  organ,  in  turn,  is  given  a  small  incre¬ 
ment  and  the  resulting  increments  of  Ve  at  all  the  field 
points  of  the  interface  are  measured.  The  result  is  an  N-by-M 
matrix,  where  N  is  the  number  of  field  points  and  M  the  number 
of  controls.  Let  this  matrix  be  called  C;  its  members  are  the 
increments  of  resulting  from  unit,  positive  increments  of 
control  settings. 


For  the  vaned-panel  controls,  a  unit  deflection  is  about 
one-tenth  the  available  vane  angle,  or  about  7.5  degrees. 


Our  procedure  for  the  measurement  of  the  control  matrix 
involved  the  same  technique  of  measuring  time-averaged 
velocity  components  as  has  been  described  above.  The  matrix 
members  were  measured  by  rows;  i.e.,  with  the  L.D.A.  set  to 
measure  the  appropriate  component,  usually  Ve,  at  a  given 
field  point  of  the  interface,  the  wall  controls,  each  in  turn, 
were  first  given  a  negative  increment  (typically  -1.0)  and 
then  an  equal  positive  increment,  and  the  component  was  mea- 
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sured.  The  difference,  from  which  the  matrix  member  is  com¬ 


puted,  is  therefore  centered  on  the  nominal  control  setting. 
Strictly,  each  control  matrix  measured  in  these  experiments 
pertains  to  a  given  model  configuration,  including  angle,  and 
a  given  array  of  control  settings,  but  it  is  expected  that 
such  a  matrix  has  validity  over  a  reasonable  range  of  control 
settings,  as  mentioned  above. 


In  Appendix  D  we  reproduce  one  of  our  BASIC  operating 
programs  for  the  measurement  of  a  control  matrix,  C.  For  the 
full  tunnel,  N  is  32  and  M  is  18,  as  has  already  been  pointed 
out;  however,  most  of  our  experiments  were  made  with  laterally 
symmetrical  model  configurations  and  under  the  assumption  of 
laterally  symmetrical  flow,  as  will  be  explained  below,  so 
that  N  and  M  were  16  and  10,  respectively.  In  either  case  N  is 
greater  than  M;  hence  M  control  settings  are  over-determined 
by  N  values  of  k*DV@.  The  control  settings  are  then  calculated 
to  give  the  best  fit  —  i.e.  the  least  mean-square  error  —  to 
the  desired  values.  This  "inversion"  of  matrix  C,  which  we 
refer  to  as  the  "best- mean-square  inversion"  of  C,  is  program¬ 
med  as  a  subroutine  in  the  operating  program  in  Appendix  D. 


If  the  result  of  this  "inversion",  an  M-by-N  matrix,  is 
called  Y,  the  formula  for  the  array  of  control  increments 


required  by  Equation  <3>  becomes 


C.  S. 


Y  *  k  *  DVe 


(6) 
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Bowel  Us  ?<  "Bt3s?-Gr  adi.ent "  Methods! 

There  are  at  least  two  different  methods  of  carrying  out 
the  iterative  procedures  of  the  adaptive-wall  scheme,  besides 
the  one  presented  above,  which  is  based  on  Equation  (3).  They 
have  both  been  tried,  at  least  minimally,  in  our  experiments 
and  therefore  will  be  sketched  here. 

A.  Dowell’s  Method  (See  also  Appendix  G.) 

The  method  described  above  will  be  seen  to  be  based  on 
the  assumption  that,  if  component  f  can  be  adjusted  as  sugges¬ 
ted  by  the  matching  discrepancy  (Equation  (3)),  the  discre¬ 
pancy  will  be  reduced,  even  though  component  q  will  also  be 
changed.  Dowell,  in  Reference  5,  has  pointed  out  that  the 
change  of  g  can  be  estimated  in  the  same  linear  approximation 
if  its  control  matrix  is  measured  as  well  as  the  control 
matrix  of  component  f.  This  algorithm  is  derived  in  Appendix  G. 

In  one  of  our  experiments  we  measured  both  control  mat¬ 
rices  and  tried  Dowell’s  method.  The  results  were  not  particu¬ 
larly  impressive  and,  in  view  of  the  tediousness  of  measuring 
two  control  matrices,  rather  than  one,  we  returned  to  the 
conventional  procedure  of  iterating  on  a  single  component, 
namely  Ve . 

B.  "RMS-Gr adi ent "  (See  also  Appendix  F. ) 

This  is  the  name  we  have  given  to  a  procedure  suggested 
by  Dr.  J.  C.  Erickson  of  Calspan  Corp.  at  A.E.D.C.  It  is  a 
very  direct  attack  upon  the  convergence  problem:  For  any  given 
model  and  control  configuration,  a  series  of  runs  is  made. 
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similar  to  those  made  to  measure  a  control  matrix;  i.e.,  in 
each  run  one  control  organ  is  given  a  unit  increment;  both 
velocity  components  are  measured  and  the  calculation  of  DV^  is 
carried  out.  If,  again,  M  denotes  the  number  of  controls,  this 
procedure  obviously  leads  to  a  set  of  M  values  of  RMS,  the 
root-mean-square  matching  error  on  the  interface  (Equation 
(4)),  or  to  a  set  of  M  values  of  D(RMS),  say,  the  RMS  incre¬ 
ments  resulting  from  individual  control  increments.  From 
these,  assuming  local  linearisation,  the  best  array  of  control 
increments  to  minimise  RMS  can  be  deduced,  as  in  Appendix  F. 

It  will  be  seen  that  the  process  of  making  M  runs,  each 
with  a  single  control  increment,  is  the  same  as  the  process  of 
measuring  both  control  matrices  in  Dowell's  method  —  the 
labor  is  identical.  Again  the  results  did  not  appear  to 
justify  the  additional  expenditure  of  time  and  work,  and  this 
procedure  was  not  adapted. 

(The  method  of  "RMS-grad 1 ent "  also  requires  recalculation 
of  the  outer  flow  for  each  run  (each  control  increment).  This 
is  not  a  serious  chore  for  low-speed  testing,  as  in  the  pre¬ 
sent  experiments,  but  can  be  time-consuming  where  more  diffi¬ 
cult  computations  are  involved,  such  as  in  transonic  flow,) 

Summary  of  "RUN"  ETocedyres 

For  practical  reasons,  all  of  the  computations  and  mani¬ 
pulations  described  in  this  report  are  carried  out  in  terms  of 
V  i  and  V£  (viz.  and  Ve )  .  per  tur  bat  ions  of  the  far-field 
values  of  the  respective  velocity  components,  as  has  already 


been  pointed  out,  whereas  the  components  measured  by  the  LDA 
are  actual  total  velocity  components.  Qur  procedure  has  there¬ 
fore  been  to  enter  the  values  of  the  unperturbed  (stream) 
velocity  components  into  the  computer  at  the  beginning  of  each 
run.  When  both  components  were  recorded  at  the  last  field 
point,  the  program  proceeded  to  calculate  the  outer  perturba¬ 
tion  flow  and  to  evaluate  and  print  out  the  distribution  DV 

e 

and  RMS.  Then,  unless  interrupted,  it  used  the  inverted  con¬ 
trol  matrix  called  C,  above,  and  then  calculated  and  printed 
the  control  increments  and  resulting  control  settings  for  a 
next  iteration. 

A  typical  computer  program  used  in  our  experiments  — 
iterations  and/or  other  runs  —  is  presented  in  Appendix  E.  In 
it  the  reader  will  find  the  "extended  averaging"  procedure, 
calculation  of  outer  flow  and  DVe ,  "Search"  procedures  (See 
below),  determination  of  new  control  settings,  etc.,  which 
have  been  described  above. 

"SEARCH"  Procedures 

Clearly,  any  set  of  measured  data  (velocity  components) 
can  be  interpreted  as  perturbations  from  any  chosen  unper¬ 
turbed  stream  components.  In  the  original  adaptable-wall 
scheme,  the  undisturbed  stream  vector,  say  (U,W),  is  to  be 
chosen,  and  the  iterative  procedure  is  to  lead  to  exact  simu¬ 
lation,  in  the  wind  tunnel,  of  flow  at  this  stream  vector.  (It 
seems  clear  that  in  high-speed  testing,  for  example,  one  would 
want  to  iterate  as  close  as  passible  to  flow  at  specified  Mach 


Number  and. angle.)  In  our  experiments,  however,  it  was  ob¬ 
served  that  a  given  set  of  experimental  data  often  matched 
better  —  i.e.,  gave  a  smaller  RMS  —  when  interpreted  as 
perturbations  from  a  stream  vector  other  than  the  one  origi¬ 
nally  specified.  In  other  words,  a  search  could  be  carried 
out,  varying  the  two  parameters  U  and  W,  and  a  smallest  RMS 
found,  for  any  given  set  of  data  (<=iny  given  run).  Inspection 
of  the  operating  program  in  Appendix  E  will  disclose  the 
details  of  this  search,  which  consisted  simply  of  adding 
constant  increments  to  V^s  and  V0s  and  recalculating  DV0  and 
RMS. 

For  most  runs,  such  searches  were  carried  out,  usually 
the  procedure  called  "SEARCHW"  (See  Appendix  E)  ,  in  which  the 
§09l®  of  the  stream  vector,  viz.  arctan(W/U),  was  held  con¬ 
stant.  Our  reasoning  was  that  in  low-speed  tests  the  speed  is 
of  less  interest  than  the  angle  of  the  stream,  and  that  a 
search  for  the  best-fit  speed,  at  constant  angle,  might  there¬ 
fore  be  an  acceptable  procedure.  In  some  cases  a  full  two- 
parameter  search,  called  SEARCHUW  (See  Appendix  E)  was  also 
carried  out. 

The  ^ssumpt  i_on  9i  Lateral  Symmetry 

If  the  wind  tunnel  and  the  model  configuration  are  both 
laterally  symmetrical,  it  might  be  expected  that  the  resulting 
flow  would  exhibit  the  same  symmetry,  in  which  case  measure¬ 
ments  would  only  have  to  be  made  in  one  half  of  the  tunnel  and 
the  wall  controls  would  be  operated  in  pairs  —  excepting,  of 
■ourse.  the  y ar 1 ab 1 e-anq 1 e  nozzle  and  the  downstream  valve. 


which  are  themselves  laterally  symmetrical.  This  situation 
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would  afford  considerable  simplification  o-f  the  experimental 
procedures  and  the  processing  of  data.  For  example,  control- 
effect  matrices  would  become  16-by-10  instead  of  32-by-18, 
with  resulting  great  reduction  of  time  required  to  measure 
them. 


This  symmetry  was  assumed  to  exist,  and  the  procedural 
modifications  mentioned  here  were  made,  in  the  early  experi¬ 
ments  of  this  program.  When  the  assumption  was  tested  by 
carrying  out  full  32-by-18  measurements  and  calculations,  the 
results  were  somewhat  surprising  and  disappointing:  It  ap¬ 
peared  that  the  basic  tunnel  flow  was  not  accurately  symmetri¬ 
cal  ,  and  that  this  asymmetry  was  amplified  by  the  high-lift 
model,  even  in  its  (nominally)  laterally  symmetrical  configu¬ 
ration.  Enough  iterations  were  carried  out  in  the  32-by-18 
mode  to  show  that  the  iteration  procedure  worked  —  i.e.,  the 
overal 1  RMS  was  substantially  reduced  by  iteration. 
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At  that  point  in  the  program  it  was  decided  that  major 
emphasis  should  be  put  on  studying  our  experimental  technique 
and,  if  possible,  learning  how  to  improve  the  accuracy  of  the 
matching;  i.e..  to  reduce  minimum  F.'MSs.  It  seemed  unlikely 
that  this  effort  required  staying  in  the  32-by-18  mode.  Our 
reasoning  was  that  the  16-by-10  half-tunnel,  running  with  a 
symmetrical  model  conf 1 gurati on  and  with  controls  operated 
symmetrical ly,  constituted  an  acceptable  demonstration  vehicle 
for  the  high-lift  adaptable-wall  scheme.  Moreover,  its  economy 
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of  time  and  labor  would  allow  us  to  carry  out  substantially 
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more  experiments. 

We  therefore  returned  to  that  mode.  Our  test  volume  in 
that  mode  is  the  East  half  of  the  tunnel;  our  measurements  and 
calculations  are  directed  toward  control  of  the  flow  in  that 
test  volume;  the  West  half  of  the  tunnel  possesses  the  same 
control  and  model  configurations,  both  of  which  affect  the 
East  half,  but  we  make  no  claims  for  the  accuracy  of  the 
simulation  in  the  West  half.  It  is  our  conviction  that  the 
same  techniques  that  succeed  in  our  half-tunnel  would  succeed 
likewise  in  the  full  32-by-18  embodiment. 

W=>ke  Surveys 

Our  procedure  for  "wake  surveys"  consisted  simply  of 
using  a  square-ended  total-head  tube  to  locate,  approximately, 
the  locus  of  maximum  total  head  in  a  vertical  plane  6.25 
inches  from  the  wall  of  the  working  section.  It  was  easy  to 
identify  the  maximum  at  locations  back  to  about  11  inches 
behind  the  wing  trailing  edge;  the  jet-wake  was  more  diffuse 
farther  aft.  The  L.D.A.  and  its  traverse  system  were  used  to 
locate  the  nose  of  the  total-head  tube. 
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PRESENTATION  OF  RESULTS 


Essentially  -four  different  kinds  o-f  experiments  were 
carried  out  in  the  course  o-f  these  investigations: 
Measurement  o-f  control  matrices 
Interations  toward  desired  -flow  conditions 
Runs  at  constant  control  settings,  varying 
angle  o-f  attack 
Wake  surveys 


In  this  section  o-f  the  report,  results  will  be  presented 
in  all  ot  these  categories,  as  well  as  some  miscellaneous  data 
and  comparisons  that  seem  to  be  pertinent  to  the  goals  of  the 
research . 


Some  typical  results  o-f  control  -matr  i  >:  measurements  are 
presented  first.  These  are  graphical  presentations  o-f  several 
repr esentat 1 ve  control  matrices;  i.e.,  bar  graphs  showing  all 
of  their  members. 


Next,  results  of  several  series  of  iterations  are  pre¬ 
sented.  These  series  have  been  chosen  to  show  results  in 
different  ranges  of  the  und i sturbed-str earn  vector  (angle  and 
magnitude)  and  angle  of  attack.  What  are  presented  in  this 
category  are  (a)  graphs  of  control  settings  and  RMS  plotted 
against  1  ter at i on-number ,  (b)  tabulated  values  of  these  data, 

(c)  tabulated  details  of  matching  error  at  the  interface,  and 

(d)  sketches  of  measured  wake  loci.  For  the  runs  made  at 
varying  angle  of  attack  and  constant  control  settings,  the 
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same  kind  of  data  is  presented  (but  of  course  there  was  no 
i  terati on )  . 

The  wake  surveys  were  made  to  determine,  at  least  rough¬ 
ly,  the  qeometry  o-f  the  energized  jet— wake  of  the  wing  with 
its  blown  flaps.  According  to  the  argument  presented  under 
OBJECTIVES,  above  (See  "S",  pages  6-7),  the  wake  qeometrv  best 
characteri zes  the  test  reqime  of  our  tunnel.  Wake-survey  data 
are  therefore  presented  for  the  same  cases  as  results  of 
iterations.  They  are  plots  of  wake  position  in  a  side-view 
sketch  of  tunnel  and  model . 

Matrix.  Measurements 

The  collected  results  of  all  the  matrix  measurements 
carried  out  in  these  experiments  constitute  a  large  volume  of 
data.  Here  we  attempt  to  present  enough  of  these  data  to  give 
some  idea  of  the  effectivenesses  of  the  various  control  ele¬ 
ments  and  how  they  are  altered  by  changes  of  parameters  such 
as  tunnel  speed,  inclination  of  the  simulated  stream,  and 
model  conf lquration.  To  this  end,  firstly  a  series  of  bar 
charts  (Figures  5.1  to  5.10)  is  presented,  showing  all  160 
matrix  members  measured  in  four  different,  typical  experiments 
under  the  assumotion  of  lateral  svmmetrv  —  640  data.  These 
data  are  collected  in  such  a  way  as  to  facilitate  comparison 
between  cases.  Secondly,  similar  bar  charts  (Figures  5.11  to 
5. SB)  are  presented  for  a  single  " asymmet i cal "  (32-bv-18)  case 

--  5^6  data.  Here  the  data  have  been  collected  in  such  a  wav 
to  show  the  svmmetrv  or  asymmetry  of  the  flow;  viz.,  matrix 


members  -for  each  field  point  and  control  organ  are  shown 


together  with  those  pertaining  to  the  laterally  opposite  field 
point  and  control  organ. 

Iter at i ons 

The  results  of  five  experiments  are  presented.  Four  of 
these  were  carried  out  under  the  assumption  of  laterally 
symmetric  flow,  as  explained  above  (See  pages  27  -  29).  Plots 
of  control  settings  and  RMS  are  presented,  showing  how  these 
quantities  varied  as  the  iteration  was  carried  out  —  from 
left  to  right  in  each  plot.  (Some  further  explanation  of  these 
graphs  is  given  below.) 

Details  of  the  resulting  match i ng-error  distribution. 

PV  .  are  presented  in  the  accomoanvina  tables. 

e 


The  fallowing  table  is  an  index  to  the  Fiqures  and  Tables 
related  to  this  series  of  experiments. 
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INDEX  TO  GRAPHS  &  TABLES  OF  ITERATION  DATA 


Exp.  Run 
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in  the  accompanying  figure  are  connected  bv  a  line. 
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In  a  number  of  cases,  on  the  contrary,  the  measured  data 
were  "searched"  by  changing  the  stream-yector  components  U  and 
W,  as  described  under  PROCEDURES,  and  the  next  iteration 
carried  out  by  calculation  from  the  resulting  velocity  compo¬ 
nents.  In  these  cases  another  paint  is  plotted  in  the  graph  of 
RMS  vs.  iteration  number,  but  it  is  not  connected  to  the 
previous  point.  If  the  iteration  was  continued  by  calculation 
of  new  control  settings  from  the  "searched"  values,  this  is 
indicated  by  the  appropriate  connecting  line  to  the  next  data 
point.  It  will  be  seen  that  occasionally  more  than  one  calcu¬ 
lation  was  made  from  a  given  data-set,  using  different  stream 
vectors  and/or  (in  at  least  one  case)  a  different  value  of  k . 

There  is  also  one  instance  (in  Figure  6.4)  where  a  run 
was  repeated  without  change  of  settings:  therefore  two  values 
of  RMS  are  shown,  with  the  same  symbol,  at  one  iteration 
number.  Data  are  tabulated  for  both  runs  (Runs  6  and  8  of 
Experiment  T70)  . 

ability  of  Data 

In  pages  1 62  to  164  the  raw  data  (L.D.A.  readings)  of 
these  two  runs  'Experiment  79,  Runs  6  and  8)  are  presented, 
together  with  data  from  two  other  pairs  of  "identical"  runs  as 
evidence  of  the  repeatability  of  our  experiments. 


The  following  is  an  index  of  these  data: 


TABLE  II 

INDEX  TO  TABLES  RELATING  TO  REPEATABILITY  OF  DATA 

Page  Exp.  Run  Model  Date 

Number  Angle 

(deg . ) 
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Runs  it  CQQitiQt  QQQtESl  Settings^.  Varying  Angie  of  Attack 
The  purpose  of  experiments  in  this  category  was  to  ex¬ 
plore  the  possibility  that,  having  arrived  at  a  "best"  control 
con-figuration  by  iteration  with  a  certain  model  con-figuration, 
one  could  vary  the  angle  of  attack  without  carrying  out  new 
iterations,  i.e.,  with  unchanged  control  settings,  and  still 
obtain  satisfactory  simulations  as  measured  by  RMS.  This  might 
indicate  that  the  control  array  was  not  overly  dependent  upon 
angle  of  attack,  for  given  stream  angle,  flap  configuration, 
etc.  Such  a  conclusion  would,  of  course,  afford  important  time 
savings  in  practical  wind-tunnel  testing. 


Two  series  of  runs  were  made  in  this  category. 
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TABLE  III 

INDEX  TO  TABLES  OF  DATA  :  RUNS  AT  CONSTANT  CONTROL  SETTINGS 

S<  VARYING  ANGLE  OF  ATTACK 


Exp.  Run  Page 
Number 


Model 
Ang  1  e 
(deg .  ) 


Free-stream  Vector 
Speed  Inclination 

(m/s)  (deg. ) 


Angle  of 
Attack 
(deg . ) 
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DISCUSSION  AND  CONCLUSIONS 


Let  us  begin  this  section  of  the  report  by  commenting  on 
the  data  presented,  under  the  same  subheadings  as  in  the 
preceding  section. 


Matrix  Measurements 

The  bar  graphs  o-f  normalized  control -matri  x  members  show 
rather  clearly,  by  their  general  similarity  of  shape,  that 
control  effects  are  not  greatly  influenced  by  changes  of  model 
and  test-section  conf i gurat i ons.  This  conclusion  must  be  tem¬ 
pered,  however,  because  it  does  appear  that  the  data  of 
"MtxMdllO"  differ  from  the  other  matrices  appreciably  at  some 
points  of  the  interface.  The  configuration  for  this  case  is 
extreme:  both  the  undi sturbed-stream  angle  and  the  jet-flap 
parameter  were  at  maximum  values.  To  maximize  the  latter,  of 
course,  the  tunnel  speed  was  at  a  minimum. 


The  graphs  also  show  that  our  control  device  Number  1, 
the  var i ab 1 e-ang 1 e  nozzle,  is  not  very  effective  at  any  field 
points.  An  ineffective  control,  of  course,  is  a  fault  in  an 
Adaptable  Tunnel.  Control  Number  3,  the  farthest  upstream 
ceiling  panel,  is  also  of  limited  effectiveness. 


i  ons 

The  general  conclusions  to  be  drawn  from  the  records  of 
iterations  are  that  the  process,  based  on  a  measured  matrix  of 

control  effects  on  V  ,  succeeds  in  reducing  the  RMS  matching- 

0 

error  and  that  the  minimum  attainable  value  of  that  Figure  of 


Merit  is  about  three  percent  of  the  simulated  stream  speed, 
usually  reached  in  about  six  to  eight  iterations  using  a 
relaxation  constant  of  0.15.  The  process  is  accelerated  by 
using  the  "Search"  technique. 

Much  of  our  time  was  devoted  to  studying  the  iteration 
process  and  the  sources  of  the  apparent  limitation  of  matching 
accuracy.  Such  sources  would  presumably  include  (1)  inaccuracy 
of  measurement  of  velocity  components,  (2)  limited  number  of 
controls,  (3)  inaccuracies  in  the  determination  of  control 
settings  to  achieve  desired  Ves,  introduced  either  by  the 
approx i mat i on  of  local  linearity  or  by  the  inaccuracy  of 
"best-mean-square  inversion"  of  the  control  matrix. 

Some  information  about  our  accuracy  of  measurement  was 
obtained  by  carrying  out  repeated  runs,  the  details  of  which 
have  been  presented  above.  These  data  show  that  our  accuracy 
is  about  one  percent  of  stream  speed.  In  somewhat  greater 
detai  1  : 

(a)  The  average  absolute  differences  in  measured  values 
in  Runs  5  and  8  of  Experiment  68  (page  l62>  amount  to  0.79 
percent  and  0.84  percent  of  nominal  stream  speed,  for 
and  ,  respectively.  These  runs  were  made  before  the 
introduction  of  the  procedure  called  "extended  avera¬ 
ging",  as  defined  in  pages  16  and  17. 

(b)  These  values  are  0.78  percent  and  0.73  percent, 
respectively,  for  Runs  6  and  3  of  Experiment  79.  These 
runs  were  made  using  "extended  averaging". 


m 


<c)  Our  third  pair  of  repeated  runs  consists  of  Run  7  of 


Experiment  71  and  Run  1  of  Experiment  79,  which  were 


carried  out  three  months  apart.  The  differences  in  these 


data  amount  to  1.78  percent  and  1.56  percent.  We  do  not 


know  the  reason  that  the  differences  are  larger  in  this 


case,  but  suspect  that  the  control  settings  may  not  have 


been  identical.  Our  records  show  that  the  control -sett i ng 


zeros  were  reset  during  the  intervening  period. 


It  appears  that  the  probable  error  in  our  measurement  of 


velocity  components,  say  one  percent,  is  large  enough  to 


account  for  a  substantial  part  of  our  residual  RMS  (about 


three  percent). 


It  is  certainly  true  that  the  availability  of  only  18 


controls  in  an  asymmetric  case,  or  10  in  a  symmetric  case. 


limits  the  minimum  RMS.  The  decision  to  use  such  a  modest 


number  of  controls  was  based  upon  the  results  of  the  computer 


simulations  that  preceded  this  experimental  investigation 


(References  3  and  4),  which  indicated  that  this  limitation 


would  not  be  too  serious. 


Evidence  of  linearity  or  nonlinearity  of  control  effects 


was  available  during  each  matrix  measurement:  viz.,  it  could 


be  observed  whether  the  increments  of  velocity  due  to  positive 


and  negative  control  increments  (See  pages  22  -  23)  were  of 


nearly  the  same  magnitude.  This  criterion  was  well  satisfied. 


lithin  the  limits  of  accuracy  discussed  above,  for  most  field 


points,  when  increments  of  plus  and  minus  one  unit  were  used 


(or,  -for  the  single  case  of  "Matrix  10",  where  the  tunnel 
speed  was  unusually  low,  plus  and  minus  two  units).  Attention 
was  given  to  the  magnitudes  o-f  control  increments  made  during 
iterations,  and  an  e-f-fort  made  to  be  sure  that  this  one-unit 
limit  was  not  exceeded.  The  reader  will  find,  by  reference  to 
the  data  pages,  that  it  seldom  was.  We  do  not  think  that 
departure  from  local  linear  behavior  was  a  serious  problem  in 
this  research. 


The  errors  due  to  "best-mean-square"  fitting  in  calcula¬ 
ting  control  deflections  were  evaluated  in  several  typical 
cases:  the  measured  control  matrix  was  used  to  calculate  the 
expected  Ve  increments  that  should  result  from  our  calculated 
( approx i mate)  control  settings,  and  these  were  compared  with 
the  values  desired.  In  every  case  the  differences  were  very 
small  —  well  within  our  experimental  error.  In  other  words,  a 
best-mean-square  fit  to  32  points  by  18  controls  (or  to  16 
points  by  10  controls)  is  a  good  approximation.  This  detail  of 
our  procedure  is  not  one  of  the  serious  limitations. 

Runs  at  Constant  Control  Settings*.  Varying  BdqI®  of  Attack 

These  data  (pages  165-172  >  probably  need  no  discussion. 

It  appears  that  the  matching  criterion  RMS  is  relatively 
insensitive  to  angle  of  attack,  for  a  given  model  configura¬ 
tion. 


Effects  of  Matching  Discrepancy 

Up  to  this  point,  the  success  of  the  Arizona  tunnel  has 


been  measured  in  terms  of  a  rather  arbitrary  Figure  o-f  Merit, 
the  RMS  matching  discrepancy  at  the  interface.  It  is  interes¬ 
ting  to  explore  the  relationship  between  this  Figure  o-f  Merit 
and  the  accuracy  o-f  simulation  at  the  model. 

Some  data  on  this  relationship  can  be  obtained  by  a 
simple  numerical  calculation,  at  least  for  representati ve 
cases. 

When  a  distribution  of  matching-error  is  present  at  the 
interface,  the  correct  boundary  conditions  are  satisfied,  both 
in  the  far  field  and  at  the  model,  but  there  is  a  discontinu¬ 
ity  at  the  interface.  If  an  equal -but-opposi te  discontinuity 
is  introduced  at  the  interface,  the  flow  field  will  be  conti¬ 
nuous  and  will  satisfy  the  far-field  conditions,  but  condi¬ 
tions  at  the  model  will  be  disturbed  by  a  field  of  extraneous 
velocities  equivalent  to  changes  of  the  models  geometry  and/or 
angle  of  attack. 

Thus  an  evaluation  of  the  tunnel’s  accuracy  can  be  ob¬ 
tained  by  introducing  these  discontinuities  at  the  interface 
and  calculating  their  field  of  velocities  at  the  model  posi¬ 
tion. 

In  the  present  case,  when  DV0  is  known  at  the  field 
points,  the  compensating  discontinuity  distribution  consists 

of  <1>  source-sink  panels  of  strength  -DV  at  the  top,  bottom, 

6 

and  front  field  points  and  (2)  distributed-vortex  panels  of 
strength  DV  at  points  one  inch  outboard  of  the  side  panels. 

9 

Computer  oroarams  have  been  constructed  to  calculate  the  velo- 
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city  -field  of  this  array  of  singularity-panels,  and  these  have 
been  used  to  calculate  the  vertical  and  horizontal  extraneous 
velocity  components  at  the  following  locations: 


mmetr 1 c 

Cases: 

Point  No 

.  > c 

y 

z  Comment 

1 

5 

0 

0 

Wing 

9 

it 

-2.5 

0 

Locati on 

ii 

-5 

0 

4 

i 

0 

0  .  .  Nose  of  Body 

5 

10 

0 

0  .  .  Tail  of  Body 

6 

•1 

0 

4  .  .  Empennage 

7 

II 

0 

-4  .  .  Wing  Wake 

iymmetr  i  c 

Cases: 

1 

5 

5 

0 

<-} 

ii 

2.5 

0  Wing 

3 

ii 

0 

0 

4 

ii 

-2  •  5 

0  Location 

5 

ii 

-5 

0 

6 

i 

0 

0  .  .  Nose  of  Body 

7 

10 

0 

0  .  .  Tail  of  Body 

8 

11 

0 

4  .  .  Empennage 

9 

II 

0 

-4  .  .  Wing  Wake 

The  results  of  these  calculations  are  presented  in  pages 


173to  177-  The  extraneous  velocity  components  at  the  model 
appear  to  be,  in  most  cases,  less  than  one  percent  o-f  stream 
speed;  there  are,  however,  a  -few  larger  values.  The  values  at 
"Point  4"  (or  "Point  6  in  the  asymmetric  case),  at  the  nose  o-f 
the  model,  are  consistently  larger;  two  values  are  greater 
than  three  percent.  This  observation  suggests  that,  in  this- 
tunnel,  it  might  be  advisable  to  locate  the  test  model  -farther 
downstream. 

Concluding  Remarks 

The  authors  o-f  this  report  believe  that  their  experiments 
have  succeeded  in  proving  that  the  concept  o-f  the  Arizona 
High-Li-ft  Wind  Tunnel  is  sound.  The  demonstration  tunnel  used 
was  designed  economically;  the  wall  controls  employed,  as  well 
as  the  configuration  of  the  blower,  diffuser,  and  settling 
chamber  reflect  this  economy.  Nevertheless,  the  performance  of 
this  tunnel,  measured  by  the  matching  accuracies  obtained,  is 
probably  adequate  for  V/STOL  testing. 

The  use  of  vaned  panels  as  wall -control  devices  seems  to 
be  quite  successful.  In  further  development  of  this  type  of 
wind  tunnel,  attention  should  be  given  to  (1)  redesign  or 
replacement  of  the  variable-angle  nozzle,  (2)  reconsi derati on 
of  the  geometry  and  arrangement  of  the  working  section  and  its 
controls,  (3)  provision  of  a  steadier,  more  uniform  basic 
tunnel  airflow,  which  may  require  redesign  of  diffuser  and 
settling  chamber  control  devices. 
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Our  success  in  using  the  L.D.A.  to  acquire  the  interface 
data  required  -for  this  type  o-f  wind  tunnel  seems  to  con-firm 
the  experience  o-f  NASA  personnel  at  Ames  Research  Center.  The 
L.D.A.,  however,  is  a  relatively  slow  instrument  —  at  least 
as  it  has  been  used  in  these  experiments.  It  appears,  to  the 
present  authors,  that  production  wind  tunnels  o-f  the  Adapt¬ 
able-Wall  category  must  be  provided  with  on-line  instrumenta¬ 
tion;  we  presume  that  L.D.A.s  can  be  adapted  to  meet  this 
requirement.  Such  tunnels  will  also  require  on-line  computa¬ 
tions  and  on-line  setting  o-f  wall  controls;  these  requirements 
do  not  seem  difficult  to  achieve  —  in  fact,  these  features 
already  exist  in  various  forms. 


1  .  Ferr  i ,  A. 
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APPENDIX  A 


PROGRAM  TO  CALCULATE  OUTSIDE-FLOW  MATRIX 
OF  COMPONENT  Ve 


1  REM  *************  IntTMatE. BAS  ******************* 

5  REM  This  pgm  uses  new  numbering,  vis.  1-4,  5—20,  21-32.  3/12/83 

10  REM  This  is  calc,  of  matrix  members;  vis.  A(32,32)  (normal  /H-outbd )  3  fid.  p 

s.  of  Interface  "T"  due  to  4  source  panels  and  28  vo 

rtex  panels  on  T.  Constructed  40985  from  IntRMatE.BAS 


20  DIM  X<3,32>,  R(3,32),  D(2),  E(2,32) ,  Y(3>,  W(3),  L(3,3,32),  Z (32) 

30  DIM  A  (32, 32)  '  norm,  /outbd .  comp. 

40  DIM  F (4) ,  G(4),  H(4)  'functions  used  in  panel  pgms.  Symbol  F()  is  used  in  bo 
h  panel  pgms. 

50  INPUT  "Input  N  =  no.  strips  for  calc.  W(3)";  N 
60  INPUT  "Input  NN  =  no.  strips  for  vortex  panels” ;NN 
70  CC  =  .0795775 
75  GOSUB  3000 


80 

FOR  Q  =  5  TO  32 

90 

L(l, 1,0)  =  1  ;  L (2, 2, 0)  =  1  : 

3,0) 

=0  :  L( 1,2,0)  =0  ;  L(2,1,Q>  =  0 

1 00 

NEXT  Q 

1 10 

FOR  Q  =  21  TO  32 

130 

L (2, 2, Q)  =  0  ;  L(3,3,Q)  =  0  : 

140 

NEXT  Q 

150 

FOR  Q  =  1  TO  4 

160 

L ( 1 , 1,0)  =0  :  L (2, 2, Q)  =  1  ; 

U 

O 

II 

o 

L ( 1 , 2, Q)  =0  ;  L(2, 1,0)  =0 

170 

NEXT  Q 

180 

FOR  Q  =  1  TO  32 

190 

L (3, 1,0)  =  — L (1,3,0)  ;  L (3, 2, 

200 

NEXT  0  :  GOTO  600  ’as  p.o.  of 

210 

FOR 

J  =  1  TO  32 

220 

LPRINT  J ; ”  " ; L ( 1 , 1 , J ) ; "  " ; L (2, 2, J ) ; " 

) ;  " 

"  ;  L (3,  1, 

J);"  " ; L  ( 1 , 3, J ) ; "  " ; L (2, 1 , J )  ;  "  "; 

) ;  " 

" ;  J 

230 

NEXT  J 

240 

STOP 

600 

FOR 

Q  =  1  TO  4  ’Front  panels 

620 

D  (  1  ) 

=  E  ( 1 ,  Q >  :  D  ( 2 )  =  E  (2,  Q) 

630 

FOR  P  =  1  TO  4  ’Front  points 

650 

GOSUB  6  0  O  0 

660 

GOSUB  2000 

670 

NEXT  P 

680 

FOR  P  =  5  TO  20  ’Ft s  on  top  it 

700 

GOSUB  6000 

710 

GOSUB  2300 

720 

NEXT  P  47 

wow 


7  3u 

FOR  P  =  21  TO 

32  ’  F'ts  on 

750 

GOSUB  6000 

760 

GOSUB  2600 

770 

NEXT  P.Q 

800 

FOR  0 

=  5  TO  32  ’All 

vortex  panel 

820 

D  ( 1  >  = 

E  ( 1 ,  Q)  :  D  (2) 

=  EC2.Q) 

830 

FOR  P  =  1  TO 

4 

840 

GOSUB  4000 

850 

GOSUB  2000 

860 

NEXT  P 

870 

FOR  P  =  5  TO 

20 

880 

GOSUB  4000 

390 

GOSUB  2300 

900 

NEXT  P 

910 

FOR  P  =  21  TO 

-rrr-l 

920 

GOSUB  4000 

930 

GOSUB  2600 

940 

NEXT  P, Q 

950  GQSUB  1097 

’As  complete  print-out  of 

1 000 

PRINT  "CHANGE  DISK  IN  B" 

:  STOP 

1010 

RESET 

1020 

OPEN  "0" , 

#1,  "B: IntTMatE 

.DAT" 

1030 

FOR  0=1 

TO  32  :  FOR  P 

=  1  TO  32 

1 040 

PRINT  #1, 

A  (P,  Q)  , 

1050 

NEXT  P.Q 

1060 

CLOSE  #1 

1070 

STOP 

sides 


desired. 


1097 

1098 

1099 

1 1 00 
1 105 
1110 
1 160 


LF'RINT  "Complete  P.Q.  af  A(P,Q)  =  IntTHatE.DAT" 

LF'RINT  "  Read  down  columns  For  Q" 

LPRINT  :  LPRINT 

F0F  P  =  1  TO  32  :  LPRINT  :  LPRINT  "P  =";P 

FCR  Q  =  1  TO  8  :  LPRINT  A(P,Q),  A(P,Q+S),  A(P,Q+16),  A(P,Q+24) 

NEXT  Q.P 

RETURN 


2000 

A (P.Q)  =  0 

2010 

FOR  L  =  1  TO  3 

2020 

A  ( P ,  Q )  =  A  ( P ,  Q )  + 

LQ.L.Q)  *W(L) 

2030 

NEXT  L 

2040 

RETURN 

2300 

A  (P,  Q)  =  0 

2310 

FOR  L  =  1  TO  3 

2320 

A  (P,  Q)  =  A  (P ,  Q ) 

+  L (3, L, Q> *W(L> 

2330 

NEXT  L 

2340 

RETURN 

2600 

A  (P,  Q) 

=  0 

2610 

FOR  L  = 

;  1  TO  3 

2620 

A  (P,  Q) 

=  A (P, Q)  +  L(l.L.Q) 

2630 

NEXT  L 

2640 

RETURN 

ww 


WW 


3000  REM  Pgm  to  input  coords,  of  panels  on  interface  T. 

3010  FOR  0  =  1  TO  4 

3020  R  d ,  0)  =0  :  E(1,Q>  =7.5  :  E(2,Q)  =7.5  :  Z(Q)  =  1 
3030  R (2, Q)  =  —7. 5*  INT ( (Q— 1 ) /2> 

3040  R (3, Q)  =  —7. 5* ( I NT (Q/2)  -  I NT ( (Q-l ) /2) ) 

3050  NEXT  Q 

3055  FOR  Q  =  5  TO  20 

3060  Z (Q)  =  -1  +  2# ( I NT ( (Q— 5) /4)  -2* INT ( (Q-5> /S) > 

(  (Q-5) /8)  :  R(3.Q>  =  -6.5*Z(Q> 

3065  E  <2.  0)  =  7.5  :  NEXT  0 


R  ( 2 . 0 )  =  -7. 5*  IN 


3072 

FOR 

0  =  5 

TO  17  STEP 

4  : 

R(1 

,0)  = 

0  : 

E  ( 1 , 0)  =  5  : 

NEXT  0 

C*U7 

FOR 

0  =  6 

TO  IS  STEP 

4  : 

Rd 

,0)  = 

5  : 

E(1,Q)  =  5  : 

NEXT  0 

3074 

FOR 

0  =  7 

TO  19  STEP 

4  : 

R  ( 1 

,0)  = 

10 

:  E  <1 ,  Q)  =  5  : 

NEXT  0 

3075 

FOR 

0  =  0 

TO  20  STEP 

4  : 

Rd 

,0)  = 

15 

:  E(1,Q>  =  10 

:  NEXT  0 

3080 

FOR  Q  = 

21  TO  32 

3090 

Z  (Q) 

=  1  - 

2*INT( (0-21) /6 ) 

:  R  (2 

0) 

=  6.5*Z(0)  : 

R  ( 3 , 0 )  = 

)  /3) 

-  2* INT ( (0 

!-21)/6))  : 

E(1 

,0) 

=  5  + 

5* ( INT ( (0+1 ) / 

3) -INT (0/3))  :  EC2.Q)  =7.5 
3100  NEXT  D 

3110  FOR  Q  =  21  TO  30  STEP  3  :  R(1,Q)  =  2.5  :  NEXT  Q 

3120  FOR  D  =  22  TO  31  STEP  3  :  R(1,0)  =  7.5  :  NEXT  Q 

3130  FOR  Q  =  23  TO  32  STEP  3  :  R(1,Q)  =  12.5  :  NEXT  Q 

3140  ’GOTO  3330  ’  as  p.o.  of  panel  coords  not  needed. 

3200  FOR  J  =  1  TO  32 

3210  LPRINT  J ; "  " ; R ( 1 , J ) , R (2, J ) , R (3, J ) , Z ( J ) , E ( 1 , J ) ; "  ";E(2.J> 

3220  NEXT  J 

3230  STOP 


3330  REM  This  is  pgm  to  input  coords  of  fid  pts  on  "T". 

3340  FOR  P  =  1  TO  4 

3350  X(1,P)  =  0  :  X(2.P)  =  R(2,P)  +  3.75  :  XC3,P)  =  R(3,P>  +  3.75 
3360  NEXT  P 

3500  FOR  P  =  5  TO  20 

3510  X(1,P)  =  R ( 1 , P )  +  .5*E(1,P)  :  X(2,P)  =  R(2,P>  +  3.75  :  X(3.P) 

3520  NEXT  P 

3530  FOR  p  =  21  TO  32 

3540  Xd,P)  =  R ( 1 , P )  +  . 5*E  d , P)  :  X(3.P>  =  R(3,P)  +  3.75 

3543  FOR  P  =  21  TO  26  :  X(2,P)  =  R(2,P)  +  1  :  NEXT  P 

3545  FOR  P  =  27  TO  32  :  X ( 2 , P )  =  R(2,P)  -  1  :  NEXT  P 

3551  ’  RETURN  ’  as  p.o.  of  pt.  coords  not  needed. 

3552  FOR  I  =  1  TO  32 

3553  LPRINT  I,  X(1,I>,  X(2,I>,  X (3, I),  I 

3554  NEXT  I 
3560  RETURN 


=  R  (3,  P) 


NEXT  P 


4000  ’  Di  str  i  buted-vorte::  panel  program. 

4010  C  =  10000 
4020  FOR  L  =  1  TO  3 
4030  V (L)  =0 
4040  FOR  K  =  1  TO  3 

4050  Y  (L)  =  Y  (L)  +  L  <K,  L,  Q)  *  (  X  <K,  P)  -R  (K,  D)  ) 

4060  NEXT  K,L 

4070  W(l)  =0  :  W (2)  =0  :  W(3)  =  0 

4000  IF  Y (3)  =  0  AND  Y(l)/D(l)  <  1  AND  Y ( 1 ) /D ( 1 )  >0  AND  Y(2)/D(2>  <  1  AND  YC2) 

D (2)  >  O  THEN  GOTO  4290  __ 

4090  H ( 1 )  =  Y(3)#Y(3)+Y(2)*Y(2>  :  H<2>  =  Y  (3) *Y (3) +C*C  :  H(3>  =  Y  (  3 ) * Y  ( 3 )  +  ( Y ( 2 ) 

D  (2)  )  *  <Y (2) -D  (2)  ) 

49 
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4100 
H  (2) 
4110 
4120 
4 1 30 
4140 
4150 
4160 
)  ) 

4 1 70 
4180 
4 1 90 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 


6000 
6010 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6090 
6 1 00 


V 

ii 


s 

I 


F  ( 1 )  = 
IF  H  ( 1 ) 
IF  H (2) 
IF  H (3) 
IF  H ( 4 ) 
IF  Y (3) 
GOTO  42 
W  ( 1 )  = 

W  (2)  = 

W  (3)  = 
NEXT  L 
RETURN 
W  ( 1 )  = 
IF  Y (2) 
W  (3)  = 
RETURN 
PRINT  " 
STOP 


0 


0 


F (2)  =  0  :  F (3)  =  0  :  F<4>  = 

<>0  THEN  LET  F ( 1 ) =Y (2) /H ( 1 ) 

<>0  THEN  LET  F(2)=C/H(2) 

<>0  THEN  LET  F  (3)  =  (D  (2) -Y  (2)  > /H  (3) 

<>0  THEN  LET  F (4) =Y ( 1 ) /H (4) 

<>0  THEN  GOTO  4240 
60 

W(l)+Y(3) *CC/NN* (G(2)/H(2>  — G  (4)/H(4)> *SGN(D(2) ) 

W(2)+Y(3) *CC/NN*(G(3) /H (3) -G ( 1 ) /H ( 1) ) *SGN(D(2> ) 

W (3) +CC/NN* (F(l) *G(1)+F(2> *G(2> +F (3) *G(3>+F(4) *G(4> ) *SGN(D(2) ) 


.5*2 (0) /D(l)  :  W (2)  =0 
<>  . 5*D (2)  OR  X (1,P)-R(1,Q> 
4*CC/ABS (D (2) ) 


. 5*D  ( 1 )  THEN  GOTO  4330 


W  ( 1 )  ="  ;  W  ( 1 )  ,  ,,W(2)="  ;W<2)  ,  "W (3) =?" 


S3  6110 
6120 
gpj  6130 
fX.  6140 
0 

6150 
6160 
6170 
) 

6130 
6190 
6200 
6210 
6220 
6230 


REM  Source-panel  program 
FOR  K  =  1  TO  3 

Y  (K)  =  0 

FOR  L  =  1  TO  3 

Y  <K)  =  Y(K) +L(L,K,Q> * (X (L,P)-R(L,Q) ) 

NEXT  L,K 

W(l)  =  0  :  W (2)  =  0  :  W(3)  =  0 
F  ( 1 )  =  SQR  ( Y ( 1 ) *Y ( 1 )  +Y  (2)  *Y  (2)  +Y  (3)  *Y  (3>  ) 

F<2)  =  SQR< (Y(l)-D(l) ) * (Y(l)-D(l) )+Y(2> *Y(2) +Y(3) *Y(3) ) 

F (3)  =  SQR ( Y ( 1 ) * Y ( 1 ) + ( Y (2) — D (2) )* (Y(2)-D(2) ) +Y (3) *Y (3) ) 

F (4)  =  SQR ( (Y(l)-D(l) ) * (Y(l)-D(l) )  +  (Y(2)-D(2) ) *  <Y(2) -D(2) ) +Y (3) *Y(3> ) 

IF  ABS  ( Y  ( 1 ) )  +ABS ( Y (3) ) <  >0  OR  Y(2>  >0  OR  Y (2) -D (2) >0  THEN  GOTO  6140 
W(l)  =  — CC*LOG  ( (Y(2)-D(2)+F(4) ) *ABS  (Y (2) -D (2) ) / ( Y (2) +F (2) ) /ABS  <Y(2)) 
GOTO  6180 

IF  ABS  (Y(l) — D ( 1 ) ) +  ABS  (Y(3))<>0  OR  Y(2)>0  OR  Y(2)-D(2)>0  THEN  GOTO  61 

W(l)  =  -CC*LOG  (ABS  (Y(2))*(Y(2) +F ( 1 ) ) /ABS (Y (2) -D (2) ) / (Y (2) -D (2) +F (3) ) ) 
GOTO  6180 

W(l)  =  — CC*LOG  (  (Y  (2)  — D  (2)  +F  (4)  )  *  (Y  (2)  +F  ( 1 )  )  /  (Y  (2)  -D  (2)  +F  (3)  )  /  (Y  (2)  +F  (2) 

IF  ABS  (Y  (2)  -D  (2)  )  +ABS  (Y  (3)  )  <  >0  OR  Y  ( 1 ) -D  ( 1 )  >0  OR  Y  ( 1 )  >0  THEN  GOTO  6210 
W (2)  =  -CC*LOG (ABS (Y(l ) ) * ( Y ( 1 ) +F ( 1 ) ) /ABS (Y( 1) -D ( 1 ) ) / (Y ( 1 ) -D ( 1 ) +F (2) ) ) 
GOTO  6350 

IF  ABS (Y (2)  ) +ABS ( Y (3) ) <  >0  OR  Y ( 1 )  >  0  OR  Y  ( 1 ) -D ( 1 )  >0  THEN  GOTO  6240 
W (2)  =  -CC*LOG ( ( Y ( 1 ) — D ( 1 ) +F (4 ) ) *ABS (Y(l)-D(l) ) / ( Y ( 1 ) +F ( 3 ) ) /ABS (Y( 1 ) ) ) 
GOTO  6350 


50 


G (2)  =  Y(2)/SQK(Y(2)*Y(2)+H(2) )  -  (Y (2) -D (2) ) /SQR ( ( Y (2) -D (2) ) * ( Y (2) -D (2) ) 

) 

FOR  L  =  1  TO  NN 

Y ( 1 )  =  X (1,P)-R(1,Q)  -  (L-.5) *D(1)/NN 
H  ( 4 )  =  Y  (3)  *Y  (3)  +  Y  ( 1 )  * Y  ( 1 ) 

G ( 1 )  =  Y ( 1 ) /SQR ( Y ( 1 ) *Y ( 1 ) +H ( 1 ) >  +  1 
G (3)  =  Y(1)/SQR(Y(1)*Y(1)+H(3))+1 

G (4)  =  Y (2) /SQR (Y (2) *Y (2) +H (4) ) - (Y (2) -D (2) )/SQR( ( Y (2) -D (2) ) * ( Y (2) -D (2) ) +H ( 


U 

0 

v" 

V 


f 


tlnalia  'a',  i.  - 


W  (2) 


6240 
) 

$  6250  IF  Y (3)  =  0  THEN  GOTO  6270 

6260  GOTO  6300 

*  6270  IF  Y(2)/D(2)  >  1  OR  Y(2)/D(2)  <  O  OF:  Y  ( 1 ) /D  ( 1 )  >  1  OR  Y  ( 1  > /D  <  1  >  <  0  THE 

GOTO  6350 

6290  W (3)  =  .5*2(0) 

6290  GOTO  6350 

(V*  6300  FOR  L  =  1  TO  N 

6310  E  =  Y<1  5)  *D<1> /N 

6320  W ( 3 )  =  W(3)-Y(3) *D<1) *CC/N/ <E*E+Y(3> *Y<3) ) * ( (Y(2)-D(2) ) /SQR <E*E+ ( Y (2) -E 

?2) ) * (Y (2>-D<2) )+Y(3) *Y<3> >-Y<2> /SQR <E*E+Y (2> *Y (2>+Y ( 

3)  * Y (3)  )  ) 

6330  NEXT  L 

6340  W (3)  =  W<3) *SGN(D<1) >*SGN(D<2) ) 

M  6350  W (2)  =  W(2)*SGN(D<1) )*SGN(D<2) > 

$  6360  W<1)  =  W(1 ) *SGN(D< 1) ) *SGN<D (2) ) 

6370  RETURN 


=  -CC*LOG< (Y(1)-D(1)+F(4))*(Y!1!+F(1)!/(Y(1)-D(1)+F(2))/(Y(1)+F(3 


m  6290 
6290 
6300 
6310 
6320 

?2) )*< 
3)*Y( 

ATTA 


WWW 


g 


APPENDIX  B 


600 

FOR  Q  =  1  TO  4  ’Fran 

620 

D  ( 1 )  =  E(1,Q)  :  D (2) 

630 

FOR  P  =  1  TO 

650 

GOSUB  6000 

660 

GOSUB  2000 

670 

NEXT  P 

680 

FOR  P  =  5  TO 

700 

GOSUB  6000 

710 

GOSUB  2300 

720 

NEXT  P 

bottom 


52 


PROGRAM  TO  CALCULATE  OUTSIDE-FLOW  MATRIX 
OF  COMPONENT  V  t 


1  REM  *********#:***  IntTMatB.BAS  ******************* 

5  REM  This  pgm  uses  new  numbering,  viz.  1-4,  5-20,  21-32.  3/12/83 
10  REM  This  is  calc,  of  matrix  members;  viz.  A  (32, 32)  (tangential)  at  fid.  pts. 
of  Interface  "T"  due  to  4  source  panels  and  28  vorte 
>:  panels  on  T.  Constructed  40985  from  IntRMatB.BAS  . 

20  DIM  X (3, 32) ,  R (3, 32) ,  D(2),  E (2,32) ,  Y (3) ,  W (3) ,  L (3, 3, 32) ,  Z  (32) 

30  DIM  A (32, 32)  ’tangential  comp. 

40  DIM  F  ( 4 )  ,  G  (4)  ,  H(4)  ’functions  used  in  panel  pgms.  Symbol  FO  is  used  in  bo 
h  panel  pgms. 

50  INPUT  "Input  N  =  no.  strips  for  calc.  W(3)“;  N 
60  INPUT  "Input  NN  =  no.  strips  for  vortex  panel s";NN 
70  CC  =  .0795775 
75  GOSUB  3000 

80  FOR  Q  =  5  TO  32 

90  L ( 1 , 1 , Q)  *  1  :  L (2, 2, Q)  =  1  :  L(3,3,Q)  =  1  ;  L(1,3,Q)  a  0  :  L (2 

3, Q)  =  0  :  L ( 1 , 2, Q>  =  0  :  L(2,1,Q)  =  0 

100  NEXT  Q 

110  FOR  Q  =  21  TO  32 

130  L (2, 2, Q)  =  0  :  L(3,3,Q)  =  0  :  L(2,3,Q)  =  -1 

140  NEXT  Q 

150  FOR  D  =  1  TO  4 

160  L ( 1 , 1,0)  =  0  :  L (2, 2,0)  =  1  :  L(3,3,Q)  =  0  :  L(1,3,Q)  =  -1  :  L( 

,3,0)  =  0  :  L ( 1 , 2 , Q >  =  0  :  L(2,1,Q)  =  0 

170  NEXT  0 

180  FOR  Q  =  1  TO  32 

190  L (3, 1 , Q)  a  -L (1,3,0)  :  L(3,2,Q)  =  -L(2,3,0) 

200  NEXT  0  :  GOTO  600  'as  p.o.  of  L’s  not  needed. 

210  FOR  J  =  1  TO  32 

220  LPRINT  J ;  "  " ;  L  ( 1 ,  1 ,  J)  :  "  '* ;  L  (2,  2,  J)  ;  “  "  ;  L  (3,  3,  J )  ;  "  "  ;  L  ( 1 , 2,  J )  ;  "  ";L(2,3, 

);"  " ; L (3, 1 , J ) ; "  " ; L ( 1 , 3, J) ; "  " ; L ( 2, 1 , J) ; "  ";L(3,2,J 
)  ;  "  "  ;  J 

230  NEXT  J 

240  STOP 


730 

FOR  P  =  21 

750 

GOSUB  6000 

760 

GOSUB  2600 

770 

NEXT  P.Q 

800 

FOR  0=5  TO  32  ' A1 

820 

D  ( 1 )  =  E(1,Q)  :  D  (2 

830 

FOR  P  =  1  T 

840 

GOSUB  4000 

850 

GOSUB  2000 

860 

NEXT  P 

870 

FOR  P  =  5  T 

880 

GOSUB  4000 

890 

GOSUB  2300 

900 

NEXT  P 

910 

FOR  P  =  21 

920 

GOSUB  4000 

930 

GOSUB  2600 

940 

NEXT  P.Q 

1000  PRINT  "CHANGE  DISK  IN  B"  :  STOP 
1010  RESET 

1020  OPEN  "0",#1,  "B: IntTMatB.DAT" 

1030  FOR  Q  =  1  TO  32  :  FOR  P  =  1  TO  32 
1040  PRINT  #1,  A (P, Q) , 

1050  NEXT  P,Q 
1060  CLOSE  #1 
1070  STOP 


1097  LPRINT  "Complete  P.O.  of  A (P.Q)  =  IntRMatB.DAT" 

1093  LPRINT  "  Read  down  columns  for  Q" 

1099  LPRINT  :  LPRINT 

1100  FOR  P  =  1  TO  32  :  LPRINT  ;  LPRINT  "P  =  ";P 

1105  FOR  Q  =  1  TO  8  :  LPRINT  A(P.Q),  A(P,Q+S),  A(P,Q+16>,  A(P.Q+24) 

1110  NEXT  Q,  P 

1160  RETURN 


2000  A  (P,  Q)  =  0 
2010  FOR  L  =  1  TO  3 

2020  A (P, Q)  =  A  (P,  Q)  +  L  (3,  L,  Q)  *W  <L) 

2030  NEXT  L 
2040  RETURN 
2300  A  (P,  Q)  =  0 

2310  FOR  L  =  1  TO  3 

2320  A (P, Q)  =  A(P.Q)  +  L ( 1 , L, 0) *W (L) 

2330  NEXT  L 

2340  RETURN 

2600  A (P.O)  =  0 

2610  FOR  L  =  1  TO  3 

2610  A (P, Q)  =  A (P, Q)  +  L(1,L,Q) *W(L) 

2630  NEXT  L 

2640  RETURN 


mmm  s 


} 


3000  REM  Pgm 
3010  FOR  0  = 
3020  R ( 1 , 0)  = 
3030  R (2, Q)  = 
3040  R (3, Q)  = 
3050  NEXT  D 
3055 
3060 

<<Q— 5>/8)  :  R 
3065  £ (2, Q)  = 

3072  FOR  Q  = 

3073  FOR  O  = 

3074  FOR  Q  = 

3075  FOR  Q  = 

3080  FOR  Q 

3090  2(0)  =  1 
>/3>  -  2* INT ( 
3) — INT (Q/3) ) 
3100  NEXT 

3110  FOR  0 

3120  FOR  Q 

3 1 30  FOR  0 

3140  ’GOTO  33 
3200 

3210 

3220 

3230  STOP 


to  input  coords.  of  panels  on  interface  RS. 

1  TO  4 

0  :  E ( 1 , Q)  =  7.5  :  E(2,Q>  =  7.5  :  Z(Q>  =  1 
—7. 5* INT ( (Q— 1 ) /2) 

—7. 5* ( INT (0/2)  -  INT( (0-l)/2) ) 

FOR  Q  »  5  TO  20 

Z(Q)  =  -1  +  2* ( INT ( (0-5) /4)  -2* INT ( (0-5 ) /S) )  :  R(2,Q)  =  -7. 5* IN 

(3,0)  =  -6. 5*Z (0) 

7.5  :  NEXT  0 

5  TO  17  STEP  4  :  R(1,Q)  =  0  :  E(1,Q)  =  5  :  NEXT  0 

6  TO  IS  STEP  4  :  R(l.Q)  =5  :  E(1,0)  =5  :  NEXT  0 

7  TO  19  STEP  4  :  R(1,Q)  =  10  :  E(1,Q)  =5  :  NEXT  O 

8  TO  20  STEP  4  :  R(1,Q)  =  15  :  E(1,Q)  =  10  :  NEXT  0 

=  21  TO  32 

-  2*INT ( (0-21 ) /6)  :  R(2.Q)  =  6.5*Z(Q)  :  R(3,Q)  =  -7. 5* ( INT ( (0-2 


E  ( 1 , 0)  =  5  : 
E(l.Q)  =  5  : 
E(l.O)  =  5 


NEXT  0 
NEXT  0 
NEXT  O 


E ( 1 , 0)  =  10  :  NEXT  0 


R ( 2 . 0 )  =  6.5*Z(G)  :  R(3,Q)  =  -7. 5* ( INT ( (0-2 


(0—21 ) /6>  >  :  E ( 1 , Q>  =  5  +  5* ( INT ( (0+1 ) / 
:  E (2, 0)  =  7.5 


=  21  TO  30  STEP  3  :  R(1,Q)  =  2.5  :  NEXT  0 

=  22  TO  31  STEP  3  :  R ( 1 , 0)  =7.5  :  NEXT  0 

=  23  TO  32  STEP  3  :  R(1,Q)  =  12.5  :  NEXT  0 

30  ’  as  p.o.  of  panel  coords  not  needed. 

FOR  J  =  1  TO  32 

LPRINT  J ;  "  ";R(1,J)  ,R(2,J)  ,  R  ( 3 ,  J )  ,Z(J)  ,  E  ( 1 ,  J )  ;  "  ,,;E(2,J) 
NEXT  J 


3330  REM  This  is  pgm  to  input  coords  of  fid  pts  on  "RS". 

3340  FOR  P  =  1  TO  4 

3350  X (l.P)  =  0  :  X  < 2, P)  =  R (2, P)  +  3.75  :  X  (3, P)  =  R(3,P>  +  3.75 
3360  NEXT  P 

3500  FOR  P  =  5  TO  20 

3510  X(1,P)  =  R ( 1 , P )  +  .5*E(1,P)  :  X(2,P>  =  R(2,P)  +  3.75  :  X(3,P)  =  R(3,P) 

3520  NEXT  P 

3530  FOR  P  =  21  TO  32 

3540  X(1,P)  =  R(1,P)  +  .5*E(1,P)  :  X(2,P)  =  R(2,P)  :  X(3,P)  =  R(3,P)  +  3.75  :  N 
XT  P 

3551  ’  RETURN  ’  as  p.o.  of  pt.  coords  not  needed. 

3552  FOR  I  =  1  TO  32 

3553  LPRINT  I,  X(1,I),  X(2,I>,  X(3,I>,  I 

3554  NEXT  I 
3560  RETURN 

4000-4340  Distributed-vortex  panel  program  :  See  Appendix  A 
6000-6370  Source-panel  program:  See  Appendix  A 


APPENDIX  C 


PROGRAM  TO  CALCULATE  OUTER-FLOW  MATRIX 


Combined  Matrix 


ExBinvT$.BAS 


10  REM  Calc  of  E;:BinvT$ 

20  INPUT  "Input  N  " ,  N 
30  DIM  A(N.N),  B (N, N) ,  C(N.N) 

50  OPEN  "I",  #1;  "B: IntT*MtE. DAT 
60  FOR  J  =  1  TO  N  :  FOR  I  =  1  TO  N  ’  <<<<<< 
70  INPUT  #1,  A ( I , J )  :  NEXT  I,J  :  CLOSE  #1 
75  STOP 


<<<  BY  COLUMNS 


<  BY  ROWS  !  !  ! 


100  OPEN  "I",  #1,  "B: BinvrsT*. DAT 

110  FOR  I  =  1  TO  N  :  FOR  J  =  1  TO  N  '  <<<<<<<<<• 

120  INPUT  #1,  B  ( I , J )  :  NEXT  J, I  :  CLOSE  #1 

125  STOP 


150  FOR  I  =  1  TO  N  s  FOR  J  =  1  TO  N  :  C<I,J)  =0 
160  FOR  K  =  1  TO  N  s  C(I,J)  =  C(I,J)  +  A ( I , K) *B (K, J) 
170  NEXT  K, J, I 


200  OPEN  "0",  #1,  "B:  E>:Bi  nvTt.  DAT 

210  FOR  I  =  1  TO  N  s  FOR  J  *  1  TO  N  ’  <<<<<<  BY  ROWS 

220  PRINT  #1,  C<I,J)  :  NEXT  J, I  :  CLOSE  #1 

225  STOP 


250  LPRINT  :  LPRINT  "Complete  P.O.  of  C(I,J)  =  ExBinvT*  by  rows" 
i.  260  LPRINT  :  LPRINT  "Read  columns  down  for  increasing  J. " 

If;  270  FOR  I  =  1  TO  N  :  LFRINT  :  LPRINT  "I  =";I 

“  2S0  FOR  J  =  1  TO  N/4  :  LPRINT  C(I,J),  C(I,J+4>,  C(I,J+8),  C(X,J+12) 
290  NEXT  J, I 
W  300  PRINT  "END"  :  END 


A??2iDIX  D 


.  PROGRAM  FOR  MEASUREMENT  AND  "INVERSION"  OF  CONTROL  MATRIX  (S2-D3TRICA L  CASE) 

I 

S^PRINT  “>>>>>>>>>>>>>>>>  MeasMtxS.  BAS  <<<<<<<<<<<<<<<<<<<<<<<" 

b  PRINT  :  PRINT  -  >>>>>>>>>>>>>  DOUBLE  INCREMENTS,  CENTERED  ««««<" 
>  10  PRINT  :  PRINT  "  PROGRAM  FOR  MAKING  MATRIX  RUNS" 

15  PRINT  :  PRINT  "  FULLY  SYMMETRICAL  CASE:  16  POINTS.  10  CONTROLS” 

M  16  PRINT  :  PRINT  “  Extended  averaging  240186“ 


17  PRINT  :  INPUT  “  Enter  NSPL  “ , NSPL 
13  GOSUB  12000  :  TRANSFER=S<HCOOO  :  NETS  =  1 

19  PRINT  :  INPUT  “  Enter  TOL  " ,  TOL 

20  BC=TF:ANSFER+3 : POKE  BC.O 
30  CALL  TRANSFER 

40  KKKK= 156115! 

70  FILES  "B:*.*“  :  PRINT 
30  INPUT“Enter  experiment  No.  “ ,  EX 
90  INPUT” Enter  date  (DDMMYY)  ".DATE 
100  DIM  V (3, 17) ,  1(17) 


110  V ( 3 ,  11)  =  EX  :  V ( 3 , 13)  =  DATE 

120  PR I NT” ENTER  INITIAL  CONTROL  SETTINGS” : PRINT 

130  FOR  K= 1  TO  10 

140  PRINT "Control  setting  #"K; : INPUT  V(3,K) 

150  NEXT  K: PRINT: PRINT 

160  LFRINT  "Initial  Control  Settings"  :  LPRINT 


AD  =  49403 1 


■V  170 

FOR  K 

■K  130 

LPRINT 

'  190 

PRINT 

m  200 

PRINT 

1  205 

PRINT 

"  210 

PR  I  NT 

300  PRINT  :  PRINT  :  INPUT  “What  kind  of  run  do  you  want  to  make? 

Enter  1  for  CHECK,  enter 

3  for  MATRIX.  ",  SEL 

310  PRINT  :  PRINT 

320  IF  SEL  =  1  THEN  GOTO  500 

340  IF  SEL  ■  3  THEN  GOTO  3900  ELSE  GOTO  300 

500  REM  *»*t*  SUBROUTINE  FOR  "CHECK"  ****:M****************** 

510  INPUT  "Enter  Run  No.  (between  .001  .999)  ",  F:X 

520  INPUT  "What  is  Number  of  Field  F'oint" ;  I 
530  PRINT 

540  INPUT  "Which  component  are  you  measuring? 

If  I  <  11,  enter  1  for  horiz.,  2  for  vert 
.  If  I  >  10,  enter  2  for  outboard,  1 

nboard.  " , CO 


for  i 


550  K  =  CO 

555  INPUT  "Enter  Chamber  pressure  SCP1. 
560  IF  I  <  3  THEN  K  =  K  MOD  2+1 
565  GOSUB  5000 

570  GOSUB  7000  _/• 

530  V  ( K ,  I )  =  CF*  AVER 


",  SCP1 


590  LPRINT  :  LPRINT  "Exp .  No.  ** ;  EX;  "  ,  Run  No.";RX;",  Comp .  "  ;  K;  "  ,  Fld.Pt.";I 
595  LPRINT  "SCP  =";SCP,  "SCP1  =  " :SCP1,  "CF  =  "  ;CF 

600  LPRINT  "  V  (  "  ;  K ;  "  ,  "  ;  I ;  " )  =" ; V (K, I ) , "  ","Standard  Dev.=";SD  :  LPRINT 
620  INPUT  "Enter  1  if  you  want  to  take  another  reading. 

Enter  0  if  you’re  through.  ",  NTHR 
630  IF  NTHR  =  1  THEN  GOTO  520 

635  IF  NTHR  =  0  THEN  GOTO  640  ELSE  GOTO  620 

640  LPRINT  :  LPRINT  :  GOTO  300 

3800  REM  **********  PROCEDURE  TO  MEASURE  MATRIX  OF  Ve  <*ymm.  case)  ********** 
3900  FOR  Z  =  1  TO  17  :  1 (Z)  *  Z  :  NEXT  Z 

3910  FOR  Z  -  7  TO  10  :  I ( Z  >  =  17  -  Z  :  NEXT  Z 

3920  FOR  Z  =  14  TO  16  :  I (Z)  =  30  -  Z  :  NEXT  Z 


3940  PRINT 
3950  PRINT 
3960  PRINT 
3970  PRINT 
3980  PRINT 
3990  PRINT 
3995  PRINT 
4000  PRINT 


:  PRINT  "If  run  must  be  i nterrupted , " 

"enter  GOTO  13000." 

"REMEMBER  THIS!!" 

:  PRINT  " _  : 

:  PRINT  "If  this  is  continuation  of  an  interrupted  run," 
"enter  GOTO  14000.  Otherwise  enter  "'1."  :  STOP 
:  PRINT 


PRINT 


STOF 


4010  INPUT  "Enter  Run  No.  (above  1000)  ", 

4015  INPUT  "Enter  SCP1  ",  SCP1 
4020  DIM  C ( 17, 13) 

4030  FOR  J  =  1  TO  13  s  C(17,J)  *  V(3,J)  : 
4040  Z  =  1  :  I  =  I <  Z )  :  J  =  0 
4050  OPEN  "O",  #1,  "B: TempMtx.DAT" 

4055  ON  ERROR  GOTO  13000 


NEXT  J 


4060  PRINT  :  PRINT  "  Go  to  Fid.  Pt.  #";I;"  &  set  optics" 

4070  PRINT  "  to  measure  HOF:IZ  component."  :  PF:INT 

4071  J=J+1  :  PRINT  :PRINT  "DEcrease  setting  of  Control  #":J  "by  an  amount  ARB." 

4072  PRINT 

4073  INPUT  "  What  is  this  amount,  ARB?  ",  ARB  :  IF  AF:B=0  THEN  GOTO  4073 


4074 

4075 
4080 
4090 
4110 
4115 
4120 

4125 
4  1 30 
4 1 40 
4145 
4147 
4 1 50 
4 1 60 


PRINT  :  PRINT  "Setting  of  Control  #";J"  should  be"; V (3, J) -ARB  :  PRINT 
GOSUB  5000 
GOSUB  7500 
VV  =  AVER 

PRINT  :  PRINT  " INcrease  setting  of  Control  #";J;"by  2  times  ARB," 
PRINT  "namely  to" ; V (3, J > +ARB  :  PRINT 


INPUT  "  Proceed,  by  entering  123,  to 

re-measure  HQRIZ  component. 

IF  ZW  =  123  THEN  GOTO  4130  ELSE  GOTO  4120 


IF  ZW 
GGSUB 
F  R I  NT 
INPUT 


"Return  Control  #";J;"to  original  setting, 
"then  enter  345  to  continue.  ",  ZW 


vie. " ; V ( 3 , J ) ", " 


IF  ZW  =  345  THEN  GOTO  4150  ELSE  GOTO  4145 
IF  J  <  10  THEN  GOTO  4071 
GOSUB  4570 


0 


4170  2  =  Z  +  1  :  I  =  I (Z)  :  J  =  0 

4180  IF  I  <  3  THEN  GOTO  4060  . 

4190  PRINT  :  PRINT  "  Go  to  Fid.  Pt.  #";I;"  ?-  set  optics" 

4200  PRINT  “  to  measure  VERT  component."  :  PRINT 

4201  J=J+1  :  PRINT  :  PRINT  "DEcrease  setting  of  Control  #";J  "by  an  amount 

4202  PRINT 

4203  INPUT  "  What  is  this  amount,  ARB?  ",  ARB  :  IF  ARB=0  THEN  GOTO  42C 

4204  PRINT  :  PRINT  "Setting  of  Control  4";J"  should  be" ; V (3, J) -ARB  :  PRINT 

4205  GOSUB  5000 
4210  GOSUB  7500 
4220  VV  =  AVER 

4240  PF:INT  :  PRINT  "INcrease  setting  of  Control  #";J;"by  2  times  ARB," 

4245  PRINT  “namely  to" ; V (3, J > +ARB  :  PRINT 
4250  INPUT  "  Proceed,  by  entering  123,  to 

re-measure  VERT  component.  ",  ZW 
4255  IF  ZW  «  123  THEN  GOTO  4260  ELSE  GOTO  4250 
4260  GOSUB  9500 

4270  PRINT  "Return  Control  #";J;"to  original  setting,  vi z . “ ; V (3, J ) " , " " 

4275  INPUT  "then  enter  345  to  continue.  ",  ZW 

4277  IF  ZW  =  345  THEN  GOTO  4280  ELSE  GOTO  4275 
4280  IF  J  <  10  THEN  GOTO  4201 
4290  GOSUB  4570 

4300  Z  =  Z  +  1  :  I  =  I (Z)  :  J  =  0 

4310  IF  I  11  THEN  GOTO  4190 

4320  PRINT  :  PRINT  "  Go  to  Fid.  Pt.  #  "  ;  I ;  "  &  set  optics" 

4330  PRINT  "  to  measure  OUTboard  horiz.  component." 

4331  J=JM  :  PRINT  :  PRINT  "DEcrease  setting  of  Control  #“:J"  by  an  amount 

4332  PRINT 

4333  INPUT  "  What  is  this  amount,  ARB?  ",  ARB  :  IF  ARB=0  THEN  GOTO  42 

4334  PRINT  :  PRINT  "Setting  of  Control  #";J"  should  be" ; V (3 , J ) -ARB  :  PRINT 

4335  GOSUB  5000 
4340  GOSUB  7500 
4330  VV  =  AVER 

4370  PRINT  :  PRINT  "INcrease  setting  of  Control  #“;J;"by  2  times  ARB," 

4375  PRINT  "namely  to " : V ( 3 , J ) +ARB  :  PRINT 
4380  INPUT  "  Proceed,  by  entering  123, 

to  re-measure  OUTbd.  horiz.  comp. 

",  ZW 

4385  IF  ZW  =  123  THEN  GOTO  4390  ELSE  GOTO  4380 
4390  GOSUB  9500 

4400  PRINT  "Return  Control  #";J;"to  original  setting,  vi z . ” ; V (3, J ) ” , " 

4405  INFUT  "then  enter  345  to  continue.  ",  ZW 
4407  IF  ZW  =  345  THEN  GOTO  4410  ELSE  4405 
4410  IF  J  <  10  THEN  GOTO  4331 
4420  GOSUB  4570 

4430  Z  =  Z  +  1  :  I  =  I  (  Z )  s  J  =  0 
4440  IF  I  <  17  THEN  GOTO  4320 
4450  GOSUB  4570 

4460  CLOSE  #1  :  ON  ERROR  GOTO  0 


ARB .  " 


ARB .  " 


S446 

ttt 


5  ’  *********1!**  COMPLETE  PRINT-OUT  OF  MATRIX  OF  Ve  ***  **  ****  *  *  *  *  *  ******** * 


4468  DIM  A ( 17, 13) 

£  4470  OPEN  "I",  #1,  "B: TempMt:: .  DAT" 


4480 

4490 


FOR  Z  =  1  TO  17  :  FOR  J  =  1  TO  13 

INPUT  #1,  A ( Z , J )  :  NEXT  J,Z  :  CLOSE  41  :  STOP 


B 


B 


assts 


4491 

4492 
4494 

4496 

4497 

4498 
4500 
4510 
4520 
4530 
4540 
4550 
4555 
4557 


FOR  I  =  1  TO  17  :  FOR  J  =  1  TO  13  :  C(I,J>  =  A(I,J)  :  NEXT  J.I 

FOR  I  =  7  TO  10  :  FOR  J  =  1  TO  13  :  C(I,J>  =  A(17-I,J>  :  NEXT  J.I 

FOR  I  =  14  TO  16  :  FOR  J  =  1  TO  13  :  C(I.J)  =  A (30-1. J)  :  NEXT 
OPEN  “0",  #1,  "B: LtStNtx.DAT 

FOR  I  =  1  TO  17  :  FOR  J  =  1  TO  13  ’  <<<<<<<<<  0(17,13)  BY  ROWS 

PRINT  #1,  C(I,J)  :  NEXT  J.I  :  CLOSE  #1 

LPRINT  :  LPRINT  “  P.O.  of  C(I.J)  =  LtStNtx .  DAT*' 

LF'RINT  "  vis.  latest  control  matrix"  :  LPRINT 

FOR  I  =  l  TO  17  :  LPRINT  :  LPRINT  "  I  =  ";  I 

FOR  J  =  1  TO  4  :  LPRINT  ,  C ( I . J ) , C  ( I , J+4) , C ( I , J+8)  :  NEXT  J 

LPRINT  ,  C ( I . 13)  :  LPRINT  :  NEXT  I 

PRINT  "End  o-f  complete  print-out  o-f  matrix  LtStNtx . DAT" 

INPUT  "Enter  123  to  proceed  with  'inversion".  ",  ZW 
IF  ZW  =  123  THEN  GOTO  4560  ELSE  GOTO  4555 


C(I,J)  =  A (30-1. J) 


0(17,13)  BY  ROWS ! ! ! 


NEXT  J, I 


4560  GOTO  6000 


4570  0(1,11)  =  EX  :  C(I,12>  =  RM  :  0(1,13)  =  DATE 

4580  FOR  J  =  1  TO  13  :  PRINT  #1,  C(I,J)  :  NEXT  J 

4585  LCTR  =  Z 

4590  LPRINT  j  LPRINT 

4600  LPRINT  :  LPRINT  "I  =“;I  :  LPRINT 

4610  FOR  J  =  1  TO  13 

4620  LPRINT  C(I,J),  :  NEXT  J 

4630  LPRINT  :  LPRINT 

4640  RETURN 


5000 
50 1 0 
5020 


INPUT  "When  LDV  is  ready,  enter  123.  ",  ZW 

IF  ZW  =  123  THEN  GOTO  5020  ELSE  GOTO  5000 
RETURN 


6000  REN  *********  BEST-NSQ  "INVERSION"  OF  0(16,10)  ******************** 
6010  DIN  B(10,10),  U(10,10),  Y(10,10) 

6020  FOR  Q  =  1  TO  10  :  FOR  P  =  1  TO  10  ;  B(F',Q>  =  0 

6030  FOR  K  =  1  TO  16  :  B(P,Q)  =  B(F',Q)  +  C  <K,  P)  *C  (K,  Q)  :  NEXT  K,P,Q 

6080  FOR  K  =  1  TO  10  :  FOR  J  =  K  TO  10  :  U(K,J)  =  B(K,J> 

6090  IF  K  =  1  THEN  GOTO  6110 

6100  FOR  P  =  1  TO  K— 1  :  U(K,J)  =  U(K,J)  -  U (K, P) #U (P, J )  :  NEXT  P 
6110  NEXT  J 

6120  FOR  I  =  K+l  TO  10  :  U(I,K>  =  B ( I , K) /U (K, K) 

6130  IF  K  =  1  THEN  GOTO  6160 

6140  FOR  P  =  1  TO  K-l 

6150  U  ( I ,  K)  =  U  ( I ,  K>  -  U  ( I ,  P >  *U  ( P ,  K )  /U  ( K. ,  K )  :  NEXT  F' 

6160  NEXT  I , K 

6170  FOR  J  =  1  TO  10  :  FOR  I  =  J  TO  10 

6180  IF  I  =  J  THEN  D  =  1  ELSE  D  =  O 

6190  Y ( I , J )  =  D 

6200  IF  J  >  1-1  THEN  GOTO  6220 

6210  FOR  P  =  J  TO  1-1  :  Y(I,J>  =  Y(I,J>  -  U  ( I ,  P)  *Y  (P,  J )  :  NEXT  F' 

6220  NEXT  I,J 

6230  FOR  J  =  1  TO  10  :  FOR  H  =  J  TO  1  STEP  -1 
6240  IF  J  *  H  THEN  D  =  1  ELSE  D  *  0 

6250  Y ( H , J )  =  D/U (H, H)  :  IF  J  <  H+l  THEN  GOTO  6280 
6260  FOR  P  =  H+l  TO  J 

6270  Y  ( H ,  J  )  =  Y  ( H ,  J  )  -  U  (H,  P)  *Y  (P,  J  ) /U  (H,  H)  :  NEXT  P 
6280  NEXT  H.J  55 


6290  DET  =  U<1,1>  :  FOR  I  =  2  TO  10 

6300  DET  =  DET *U (1,1)  :  NEXT  I 

6310  PRINT  "DET  =";DET 

6320  FOR  J  »  1  TO  10  i  FDR  I  =  1  TO  10  :  U(I,J>  =  0 

6330  FOR  K  =  1  TO  10 

6340  IF  K  =  J  THEN  YY  =  1  ELSE  YY  =  Y(K,J> 

6350  IF  K  <  I  OR  K  <  J  THEN  GOTO  6370 

63o0  U  <  I ,  J  >  =  Li  ( I ,  J )  +  Y  ( I ,  K)  *YY 

6370  NEXT  K,I,J 

6380  ’  GOSUB  8000  it  p.o.  o -f  B*Binv  is  needed 
6390  ERASE  B,  Y  :  DIM  Y(13,16) 

6400  FOR  I  =  1  TO  16  :  Y  ( 1 1 , I )  =  EX  :  Y<12,I>  =  RM  :  Y(13,I>  =  DATE 

6410  FOR  J  =  1  TO  10  :  Y(J, I)  =  0 

6420  FOR  K  =  1  TO  10  :  Y<J,I>  =  Y<J,I)  +  U  <  J  ,  K>  *C  <  I ,  K) 

6430  NEXT  K.J,I 


6440  OPEN  “0",  #1,  "B: VEintoC5.DAT"  ’  <<<<<<<<<  Y(13, 16)  BY  COLUMNS!!! 

6450  FOR  1=1  TO  16  :  FOR  J  =  1  TO  13 

6460  PRINT  #1,  Y  ( J ,  I )  :  NEXT  J,I  :  CLOSE  #1 

6470  LPRINT :  LPRINT  :  LPRINT  "  Complete  p.o.  o-f  Y  =  VEintoCS.DAT 

6480  LPRINT  ”  £*p.  EX;  "  Run  #" ;  RM{  "  Date";  DATE 

6490  FOR  J  =  1  TO  10  :  LPRINT  :  LPRINT  "  Control  #,  J 

6500  FOR  1=1  TO  4  :  LPRINT  ,  Y  ( J ,  I ) , Y  (J , 1+4) , Y < J , I +8) , Y < J , 1  +  1 

6510  NEXT  I,J 

6540  PRINT  "End  o-f  complete  p.o.  o-f  VEi ntoC* .  DAT" 

6550  PRINT  :  FRINT  " VEi ntoC* . DAT  is  saved  on  B. " 

6560  PRINT  ;  INPUT  "Enter  123  to  return  to  Line  300, 

viz.,  What  kind  o-f  run  ..?  ", 

ZX 

6570  IF  ZX  =  123  THEN  GOTO  300 
6580  FRINT  "END  OF  RUN"  :  END 


7000  REM  *************  Subroutine  -for  CHECK  ***************** 
7005  GOSUB  10000 

7010  INPUT  "Enter  SCP  ",  SCP 

7020  CF  =  SQR ( SCP 1 / SCP ) 

7025  AVER  =  W(CR> 

7030  PRINT  :  PRINT  s  PRINT  "Average  VEL  =";AVER,  "Std.  Dev.  =";SD 


7040  PRINT  ! 
7050  PRINT  ' 
7060  FRINT  : 

",  F'Q 

7070  IF  PQ  = 
7075  IF  PQ  = 
7080  PRINT  : 
7090  RETURN 


:  PRINT  "SCP  =" 
"Corrected  Vel . 
:  INPUT  " 


SCP,  "SCP1  =" ;SCP1,  "CF  =" ; CF 
=";  CF  *  AVER 

To  proceed,  enter  9;  to  repeat  measurement 
o-f  vel.  comp.,  enter  8 


9  THEN  GOTO  7080 

8  THEN  GOTO  7005  ELSE  GOTO  7060 
PRINT  " - 


PR  I  NT 


/.v.v.v 


/.V ■/ 


/ 


.v  .-  .-.-.-.v;.- 


r  j  <-> 


7500  ’  Subr.  for  MATRIX  a************* ********** 

7505  GOSUR  10000 
7507  AVER  =  W(CR) 

7510  PRINT  "Average  VEL  =";AVER,  "Std.  Dev.  =";SD 
7515  PRINT  "  _ " 

7520  PRINT  :  INPUT  "  To  repeat  measurement  of  vel.comp., 

enter  8;  to  proceed,  ent 

er  9 .  ",  PQ 

7530  IF  PQ  =  9  THEN  GOTO  7540 

7535  IF  PQ  =  8  THEN  GOTO  7500  ELSE  GOTO  7520 

7540  PRINT  :  PRINT  “ - 


•"  :  PRINT 


7550  RETURN 


8000  ’  *************  Check:  of  inversion  by  multiplication  *************** 

8010  DIM  Z(10,10>  :  FOR  J  =  1  TO  10  :  PRINT  "J  =";J 

8020  FOR  I  =  1  TO  10  :  Z(I,J)  =0  :  FOR  K  =  1  TO  10 

8030  Z(I,J)  =  Z(I,J)  +  B ( I , K) *U (K, J)  :  NEXT  K 

8040  PRINT  Z(I,J),  :  NEXT  I  :  PRINT  :  NEXT  J 

8050  RETURN 


.‘•I  9500  GOSUB  7500 

'■  9510  C(I,J)  =  (AVER  -  VV) /ARB/2 

9520  RETURN 

|  FOR  SUBROUTINE  TO  MEASURE  VELOCITY  &  CALCULATE  AVER  &  STD  DEV,  SEE  APP2JDH  E 
(  Lines  9990  -  12640) 


13000  ’  *************  Error  Handling  ************** 

13010  CLOSE  #1 

13020  PRINT  :  PRINT  "INTERRUPTED  MATRIX  RUN" 

13030  PRINT  "Currant  I  and  J  are"  I  "  and  “  J 

13040  PRINT  :  PRINT  "The  last  counter  Z  for  which  row" 

13050  PRINT  "C(Z,J>  was  stored  on  disk  is";LCTR 
13060  END 

14000  REM  *******  Procedure  for  restarting  interrupted  MATRIX  run  *********** 
14010  INPUT  "Enter  value  LCTR  of  aborted  run.  " , Z 
14020  I  =  I ( Z ) 

14030  DIM  C  < 1 7 , 13) 

14040  OPEN  "I",  #1,  "B: TempMtx . DAT 

14050  FOR  ZZ  =  1  TO  Z  :  FOR  J  =  1  TO  13 

14060  INPUT  #1,  C(ZZ,J)  :  NEXT  J,ZZ  :  CLOSE  #1 

14070  PRINT  "Check  some  values  of  C(ZZ,J),  then  use  "'1  to  continue."  :  STOP 

14080  OPEN  "0",  #1,  "B:  TempMt:: .  DAT 

14090  ON  ERROR  GOTO  13000 

14100  FOR  ZZ  =  1  TO  Z  :  FOR  J  =  1  TO  13 

14110  PRINT  #1,  C(ZZ,J)  :  NEXT  J,ZZ 

14120  Z  =  Z+l  :  I  =  I ( Z )  :  J  =  O 

14130  IF  Z < 3  THEN  GOTO  4180 

14140  IF  Z<11  THEN  GOTO  4310 

14150  IF  Z < 1 7  THEN  GOTO  4440  ^ 


&&&&&& 


mm 


APPENDIX  E 


1 


PROGRAM  FOR  "CHECK"  READINGS  AND  FOR  RUNS 


JV* 

5  REM  > > >  >  >  >  >  >  > >  >  WTXp  $m 1 8 . BAS  <  <  <  <  < < <  <  <  < <  < < <  <  < <  <  <  < < 

■  6  PRINT  "  This  is  WTXp$ml8.BAS  -  -  230186"  :  PRINT 

•7  7  PRINT  "PROGRAM  FOR  WIND-TUNNEL  RUNS:  DATA  ACQUISITION  ETC."  :  PRINT 

°  8  PRINT  "  SYMM.  CASE.  16  FLD  F'TS.  10  CONTROLS"  :  PRINT 

9  PRINT  "  with  extended  averaging  of  VEL  readings."  :  PRINT 


10  PRINT  :  PRINT  "For  re-calc.  w.o.  measurements,  GOTO  70"  :  STOP 
r.  15  GOSUB  12000:  TRANSFER=&HC000:  NBTS=1 :  AD=49408  ! 

17  PRINT  :  INPUT  "  Enter  NSPL  ",  NSPL 

13  PRINT  :  INPUT  "  Enter  TOL  ",  TOL 

v.  20  BC=TRANSFER+3 :  POKE  BC,0 
2  30  CALL  TRANSFER 
“  40  KKKK= 1 56 115! 

70  PRINT  .-FILES  "B:*.*“  :  PRINT 
y  SO  PRINT  :  INPUT "Enter  experiment  No.  ",  EX 
V  90  INPUT "Enter  date  (DDMMYY)  " , DATE 
100  DIM  V(3,16),  1(17),  W (25) 

,  105  PRINT  :  PRINT  "For  re-calc  of  data  GOTO  3710  or  3000"  :  STOP  :  PRINT 

i 

110  V(3,ll)  =  EX  :  0(3,13)  =  DATE 
..  120  PRINT-ENTER  INITIAL  CONTROL  SETTINGS" :  PRINT 

•X  125  PRINT  "(To  skip,  GOTO' 300)"  :  PRINT  :  STOP 
130  FOR  K=1  TO  10 

140  PR I NT "Control  setting  #"K; : INPUT  0(3, K) 

■  150  NEXT  K: PRINT: PRINT 

*  160  LPRINT  "Initial  Control  Settings"  :  LF'RINT 

170  FOR  K  =  1  TO  13  :  LPRINT  K,  0(3, K) ,  K  :  NEXT  K 

v.  130  LPRINT  :  LPRINT 

190  PRINT  "Please  check  initial  control  settings." 

200  PRINT  "To  change  Kth  C.S.,  write  ’  0  (3,  K)  =.  .  .  '  " 

210  PRINT  "To  proceed,  use  1 .  "  :  STOP 

300  INPUT  "What  kind  of  run  do  you  want  to  make? 

Enter  1  for  CHECK,  enter  2  for  RUN  " , SEL 
310  PRINT  :  PRINT 
“  320  IF  SEL  =  1  THEN  GOTO  500 
.  330  IF  SEL  =  2  THEN  GOTO  1000  ELSE  GOTO  300 


500  REM  *****  SUBROUTINE  FOR  "CHECK"  ************************ 
510  INPUT  "Enter  Run  No.  (between  .001  •!<  .999)  ",  RX 

520  INPUT  "What  is  Number  of  Field  Point";  I 
530  PRINT 


540  INPUT  "Which 
If  I 

",C0 


component  are  you  measuring? 

11,  enter  1  for  horiz.,  2  for  vert 

If  I  >  10,  enter 


for 


outboard , 


nboard 


62 


SCP1 


IT. 

i 


V 

.v 


R 


% 

1L* 


8 


.V 


•>: 

v, 


550 

555 

560 

565 

570 

580 

590 

595 

600 

610 

620 


I 


5 


V 


630 

635 

640 


K  =  CQ 

INPUT  "Enter  Chamber  Pressure  SCP1.  ", 

IF  I  <  3  THEN  K  =  K  MOD  2  +  1 
GOSUB  5000 
GOSUB  7000 
V ( K , I )  =  CF*AVER 

LPRINT  :  LPRINT  "Exp . No. " ; EX ; " ,  Run  No.";RX;",  Comp.";K;",  Fld.Pt.":I 
LPRINT  "SCP  =" ;SCP,  "SCF'l  =";SCP1,  "CF  =";CF  ,  "CR  =  ";CR 
LPRINT  " V ( " ; K; " , " ; I ; " )  =  " ; V (K, I ) ,  "  "."Standard  Dev.=";SD  :  LPRINT 
LPRINT 

INPUT  "Enter  1  if  you  want  to  take  another  reading. 

Enter  0  i-f  you're  through.  ",  NTHR 
IF  NTHR  =  1  THEN  GOTO  520 
IF  NTHR  =  0  THEN  GOTO  640  ELSE  GOTO  620 
LPRINT  :  LPRINT  :  GOTO  300 


UU 

WW 


1000  REM  ******  SUBROUTINE  FOR  "RUN"  ************** 

1010  INPUT  "Enter  Run  No.  minus  1  ",  RN 

1020  INPUT  "Enter  UU  =  x  comp,  of  stream  vel.  ", 

1030  INPUT  "Enter  WW  -  2  comp,  of  stream  vel.  ", 

1040  INPUT  "Enter  relax,  factor  kk.  ",  KK 
1045  INPUT  "Enter  SCP1  ",  SCP1 

1050  RN  =  RN  +  1  :  0(3,12)  =  RN  ;  V(3,14>  =  UU  : 

1060  LPRINT  "This  is  Exp.  #";EX;",  Run  #";RN 

LPRINT  "UU  =" ; UU,  "WW  =";WW,  "k  =";KK,  "SCP1  =";SCP1 


0(3, 15)  =  WW 


106J 


LPRINT 


1070 

1071 

1072 

1073 

1074 
1 080 


FOR  CTR 
FOR  CTR 
FOR  CTR 
FOR  CTR 
FOR  CTR 


1  TO  3  : 
4  TO  6  : 
7  TO  8  : 
9  TO  12 
13  TO  16 


I (CTR)  =  CTR  +  10 
I (CTR)  =  20  -  CTR 
I (CTR)  *  9  -  CTR  : 
I (CTR)  =  CTR  -  6 


;  NEXT  CTR 
:  NEXT  CTR 
NEXT  CTR 
:  NEXT  CTR 


I (CTR) 


23  -  CTR 


NEXT  CTR 


1(17)  =  17 


1085 

CTR  = 

1  :  I  =  11  :  K 

1090 

PRINT 

"Go  to  fid.  pt 

1 100 

GOSUB 

5000 

1110 

GOSUB 

9000 

1 120 

IF  K  = 

>  1  THEN  K  =  2 

1 1 30 

PRINT 

"Move  optics  1 

1140 

GOSUB 

5000 

1150 

GOSUB 

9000 

1160 

CTR  = 

CTR  +1:1= 

1170 

PRINT 

"Go  to  fid.  pt 

1  180 

GOSUB 

5000 

1190 

GOSUB 

9000 

1200 

IF  K  = 

'  1  THEN  K  =  2 

1210 

PRINT 

"Move  optics  1 

1220 

GOSUB 

5000 

1230 

GOSUB 

9000 

1240 

CTR  = 

CTR  +1:1= 

1250 

IF  I  < 

>  2  THEN  GOTO 

1 255 

IF  I  = 

=  2  THEN  K  =  2 

1260 

IF  I  = 

:  3  THEN  K  =  1 

1270 

PRINT 

"Go  to  Fid.  Pt 

1 280 

GOSUB 

5000 

1290 

GOSUB 

9000 

1 

1 


to  measure  HORIZ  component" 


63 


HT» 

Bi 1 


Jfc  142C 
S  143C 


IF  K  =  1  THEN  K  =  2  ELSE  K  =  1 

PRINT  "Rotate  optics  to  measure  VERT  component." 

GOSUB  5000 
GOSUB  9200 

CTR  =  CTR  +  1  :  1=1 (CTR) 

PRINT  "Go  to  fid.  pt.  #";I;"to  measure  VERT  component." 
GOSUB  5000 
GOSUB  9200 

IF  K  =  1  THEN  K  =  2  ELSE  K  =  1 

PRINT  "Rotate  optics  to  measure  H0RI2  component." 

GOSUB  5000 
GOSUB  9000 

CTR  =  CTR  +1:1=  I (CTR) 

IF  CTR  <  17  THEN  GOTO  1255 


1435  V (3, 14)  =  UU  :  V ( 3 , 1 5 )  =  WW 

1440  LPRINT  “I",  "Vt  =  V(1,I)",  "Ve  =  V(2,I>",  "CS  =  V(3,I)",  "I"  :  LPRINT 

1450  FOR  I  =  1  TO  16  :  LPRINT  I,  V(1,I>,  0(2,1),  V(3, I),  I  :  NEXT  I 

1453  BETA  =  ATN ( WW/UU) *57. 3 

1460  LPRINT  :  LPRINT  :  LPRINT  "UU  =";UU,"WW  =";WW,"BETA  =" ; BETA,  "k=  " :  KK 
1465  LPRINT  :  LPRINT  :  PRINT  "To  skip  saving  V ( I , j ) s, GOTO  2000"  :  STOP 

1470  OPEN  "0",  #1,  "B: TempVel *. DAT" 

1400  FOR  K  =  1  TO  3  :  FOR  I  =  1  TO  16  :  PRINT  #1,  V(K,I) 

1490  NEXT  I , K  :  CLOSE  #1 


m 

1500 

PRINT 

"  V  ( I  „  J  ) 

1505 

PRINT 

:  PRINT 

jg 

1510 

PRINT 

:  PRINT 

2000 

REM  *******  1 

P 

2010 

DIM  L ( 16, 16) , 

7- 

2020 

OPEN  " 

I",  #1, 

2030 

FOR  I 

=  1  TO 

« • 
> 

2040 

INPUT 

#1,  L ( I 

$:• 

2050 

LPRINT  "I" 

2060 

FOR  I 

=  1  TO 

2000  REM  *******  CALCULATION  OF  NEW  Ve  ************** 

2010  DIM  L(16,16),  D(16),  X (13, 16) ,  F(16),  K(13) 

2020  OPEN  "I",  #1,  "A:  ExBinvT* .  DAT" 

2030  FOR  I  =  1  TO  16  :  FOR  J  =  1  TO  16 
2040  INPUT  #1,  L ( I , J )  :  NEXT  J, I  :  CLOSE  #1 
2050  LPRINT  "I",  "D ( I ) " ,  "Dsrd  Ve  (k  = " ; KK ; " ) "  :  LPRINT 
2060  FOR  I  =  1  TO  16  :  D(I)  =  -V (2, I ) 

2070  FOR  J  =  1  TO  16 

2000  D ( I )  =  D ( I )  +  L  ( I , J ) *V  ( 1 , J )  :  NEXT  J,I 

2090  MSO  =  O  :  FOR  P  =  1  TO  16 

2100  MSQ  =  MSQ  +  D(P)*D(P)/16  :  NEXT  P  :  RMS  =  SQR  (MSQ) 

2110  FOR  I  =  1  TO  16  :  F(I>  =  V(2,I)  +  KK*D(I) 

2115  VV  =  SQR(UU*UU+WW*WW) 

2120  LPRINT  I,  D ( I ) ,  F(I),  I  :  NEXT  I  :  LPRINT 

2130  LPRINT  "RMS  =";  RMS,  "VV  =";VV,  "RMS  percent  =" ; RMS/VV* 100 
2135  LPRINT  :  LPRINT 

2140  PRINT  :  PRINT  "To  continue  iteration,  GOTO  2500." 

2150  PRINT  :  PRINT  "To  exercize  '*EARCHVV;,  GOTO  3800." 

2155  PRINT  "To  exercise  ’*EARCHUW',  GOTD  3810 
2160  STOP 


mmm 


2500  REM  **********  CALCULATION  OF  CONTROL  INCREMENTS  ************* 
2516  OPEN  "I",  #1,  "B:  CM;: Mod  19.  DAT" 

2520  FOR  I  =  1  TO  16  :  FOR  J  =  1  TO  13  :  INPUT  #1,  X(J,I> 

2530  NEXT  J, I  :  CLOSE  #1 

2540  LPRINT  "J",  "Cntrl . Incr . ” ,  "New  C.S.",  "J"  :  LPRINT  :  LPRINT  " 

=" ; KK; " ) "  :  LPRINT 

2550  FOR  J  =  1  TO  10  :  K(J>  =  0  :  FOR  I  =  1  TO  16 
2560  K  ( J )  =  K ( J )  +  X (J, I) *KK*D (I)  :  NEXT  I 
2570  V  ( 3 ,  J )  =  V(3,J)  +  K(J>  :  NEXT  J 

2580  FOR  J  =  1  TO  13  :  LPRINT  J,  K(J),  V(3,J>,  J  :  NEXT  J  :  LPRINT 
2585  GOSUB  8000 

2590  PRINT  "Reset  Controls;  then,  to  run  next  iteration," 

2591  PRINT  "GOTO  2610."  :  PRINT 
2595  STOP 

2610  ERASE  L,D,X,K,F  ’  ««<««  These  are  re-DIMed  in  2010 
2615  GOTO  1050 


3000  ’  **********  To  input  known  V(I,J)s  ******************************* 

3010  FOR  K  *  1  TO  3  ;  FOR  1  =  1  TO  16 

3020  PRINT  " V ( " j K; " , “ ; I j " )  =";:  INPUT  V(K,I> 

3030  NEXT  I , K 
3040  STOP 
3700  DIM  V (3, 17) 

3710  OPEN  "I",  #1,  "B: TemOl 120. DAT 

3720  FOR  K=1  TO  3  :  FOR  1=1  TO  16  :  INPUT  #1,  VCK,I> 

3730  NEXT  I.K  ;  CLOSE  # 1 

3745  PRINT  :  PRINT  "Present  UU.WW,  &  k  are";UU;"  ";WW;"  S<  " ;  KK 
3750  PRINT  j  PRINT  "To  exercise  ’*EARCHW’,  GOTO  3800."  :  PRINT 
3752  PRINT  "To  exercise  ’*EARCHUW’,  GOTO  3810 
3755  PRINT  ;  PRINT  "To  print-out  VCI,J)s,  GOTO  1435." 

3760  STOP 


3800  INPUT  "Input  new  UU  ",  UUU 
3802  WWW  =  UUU*TAN (BETA/57. 3) 

3804  GOSUB  3820 
3806  GOTO  3800 

3810  INPUT  "Input  new  UU  ",  UUU 

3812  INPUT  "Input  new  WW  ",  WWW 

3813  BETA  =  57. 3*ATN ( WWW/UUU) 

38 1 4  GOSUB  3820 

3816  GOTO  3810 

3820  AVHQR  *  UUU  -  UU  ;  AVVRT  =  WWW  -  WW 

3830  GOSUB  4400 

3850  UU  =  UUU  :  WW  =  WWW 

3860  LPRINT  :  LPRINT 

3870  LPRINT  "UU  =";UU,  "WW  =";WW,  "BETA  =";  BETA 
3880  W  =  SQR(UU*UU+WW*WW) 

4000  FOR  I  =  1  TO  16  :  D(I)  =  -V(2,I) 

4010  FOR  J  =  1  TO  16 

4020  D ( I )  =  D ( I )  +  L ( I , J ) *V  < 1 , J )  :  NEXT  J, I 
4030  MSQ  =  O  ;  FOR  P  =  1  TO  16 

4040  MSO  =  MSQ  +  D(P)*D(P)/16  :  NEXT  F  :  RMS  =  SQR(MSQ) 

4050  LPRINT  "RMS  =";RMS,  "W  =";VV,  "RMS  percent  =";RMS/W*100 

4060  RETURN 


4400  FOR  I 

4410  FOR  I 

4420  FOR  I 

4430  RETURN 
4600 
4610 


I  TO  2  :  A  =  1  i  B  =  2  :  GOSUB  4600  :  NEXT  I 
3  TO  10  :  A  =  2  :  B  =  1  :  GOSUB  4600  :  NEXT  I 

II  TO  16  :  V(l,  I)=V(1,  D-AVHOR  :  V(2,I>  =  V(2,I)-AVH0R 

V  ( B ,  I )  =  V(B,  D-AVHOR  :  V(A,I>  =  V(A,I)-AVVRT 
RETURN 


NEXT  I 


5000  INPUT  "When  LDV  is  ready,  enter  123.  ",  ZW 

5010  IF  ZW  =  123  THEN  GOTO  5020  ELSE  GOTO  5000 

5020  RETURN 


7000  REM  **************  Subr.  -far  CHECK  8<  RUN  ******************* 

7005  GOSUB  10000 

7010  INPUT  "Enter  SCP  ",  ZYX 

7015  IF  ZYX  *  O  THEN  GOTO  7020 
7017  SCP  =  ZYX 

7020  CF  =  SQR ( SCP 1 / SCP ) 

7025  AVER  =  W(CR) 

7030  PRINT  :  PRINT  :  PRINT  "Average  VEL  =" ; AVER,  "Std.  Dev.  =";SD 

7040  PRINT  :  PRINT  "SCP  =";SCP,  "SCP1  =";3CP1,  "CF  =";CF 

7050  PRINT  :  PRINT  "  >>>>>>  Corrected  Vel .  =";CF*AVER"  <<<<<" 

7055  PRINT  "  - " 

7060  PRINT  :  INPUT  "  To  repeat  measurement  of  vel.  comp.,  enter  S; 

to  proceed,  enter  9.  " 

,  PQ 

7070  IF  PQ  =  9  THEN  GOTO  7080 

7075  IF  FQ  =  8  THEN  GOTO  7000  ELSE  GOTO  7060 

7080  PRINT  :  PRINT  " - "  :  I 

7090  RETURN 


■"  :  PRINT 


9000  GOSUB  7000 
9010  V ( K , I )  =  CF*AVER 
9020  RETURN 

9200  GOSUB  7000 
9210  V (K, I )  =  CF*AVER 
9220  RETURN 


9990  REM  ***********  SUBR  TO  MEASURE  VEL  &  CALC  AVER  ETC.  ************* 

10000  CR  =  0  :  SUM3  =  0 

10010  CR  =  CR  +  1 

10012  SUM  =  0  :  SUM2  =  0 

10015  FOR  LOOP=AD  TO  < AD+  <3*NSPL) -1 )  STEP  3 

10016  HB=INT (LOOP/256) : LB=LOOP- (HB*256 ) 

10017  POKE  49185! , LB: POKE  49186 ! , HB 
10020  POKE  BC, NBTS 

10030  CALL  TRANSFER 
10035  NEXT  LOOP 


10040  FOR  LQOP=AD  TO  ( AD+(3*NSPL) -1 )  STEP  3 
1 0050  MANT=256  * ( PEEK ( LOOP+ 1 )  AND  1 5  >  +PEEK ( LOOP ) 

1 0060  XFON= I  NT  ( PEEK < LOOP+ 1 ) / 1 6 ) 

1 0070  C=PEEK ( L00P+2 ) 

100S0  VEL= (0/ (WANT* (2A (XPON-3) ) *2E-09> ) /KKKK 
10085  ’  PRINT  VEL 
10090  SUM=SUM+VEL 
10100  SUM2=SUM2+(VEL*VEL) 

10110  NEXT  LOOP 
10120  A VER=SUM / NSF'L 

10122  SUM3  =  SUM3  +  AVER 

10123  W(CR)  =  SUM3/CR 

10124  PRINT  CR,  AVER,  W(CR) 

10125  IF  ABS(W(CR)-W(CR-1) XTOL  THEN  GOTO  10126  ELSE  GOTO  10010 

10126  IF  ABS  (W(CR)-W(CR-2)  XTOL  THEN  GOTO  10127  ELSE  GOTO  10010 

10127  IF  ABS ( W ( CR ) -W ( CR-3 ) ) < TOL  THEN  GOTO  10130  ELSE  GOTO  10010 

10128  ’  (Omit  for  4  rdgs) IF  ABS (W (CR) -W (CR-4) ) <TOL  THEN  GOTO  10130  ELSE  GOTO  1 
10 

10130  SD=SQR (ABS ( (SUM2/NSPL) - (AVER*AVER) ) ) 

10170  PRINT  :  RETURN 


12000  FOR  LK=49 152 !  TO  49267! 

12010  READ  VL: POKE  LK, VL: NEXT  LK 
12020  RETURN 

12030  REM  <««  MAIN  SECTION  (63  BYTES  LONG)  >»» 


12040 

DATA 

62,0 

REM 

LD  A , OOH 

12050 

DATA 

6,0 

REM 

LD  B, BCNT 

12060 

DATA 

184 

REM 

CP  A,  B 

12070 

DATA 

Am  ^  W 

REM 

JR  N2, SAMPLE 

12080 

DATA 

14,  3 

REM 

LD  C, 03H 

12090 

DATA 

205,63, 192 

REM 

CALL  SCREG 

12100 

DATA 

14, 117 

REM 

LD  C, 75H 

12110 

DATA 

205,63, 192 

REM 

CALL  SCREG 

12120 

DATA 

,0,0 

REM 

LD  HL, OOOOH 

12130 

DATA 

62 , 255 

REM 

LD  A, FFH 

12140 

DATA 

35 

REM 

INC  HL 

12150 

DATA 

188 

REM 

CP  A,  H 

12160 

DATA 

32,252 

REM 

JR  NZ ,  LOOF’l 

12170 

DATA 

14,21 

REM 

LD  C, 15H 

12180 

DATA 

205,63, 192 

REM 

CALL  SCREG 

12190 

DATA 

201 

REM 

RET 

12200 

DATA 

Ow,  0,  19^> 

REM 

LD  HL, A200H 

12210 

DATA 

14,85 

REM 

LD  C, 55H 

12220 

DATA 

205,63, 192 

REM 

CALL  SCREG 

12230 

DATA 

205,84, 192 

REM 

CALL  RSIN 

12240 

DATA 

113 

REM 

LD  (HL) , C 

12250 

DATA 

35 

REM 

INC  HL 

12260 

DATA 

205, 84, 192 

REM 

CALL  RSIN 

12270 

DATA 

113 

REM 

LD  (HL) , C 

12280 

DATA 

35 

REM 

INC  HL 

12290 

DATA 

205,84, 192 

REM 

CALL  RSIN 

12300 

DATA 

113 

REM 

LD  ( HL ) , C 

12310 

DATA 

:  REM 

INC  HL 

12320 

DATA 

16. 239 

:  REM 

DJNZ  LOOF'2 

12330 

DATA 

14,21 

:  REM 

LD  C, 15H 

12340 

DATA 

205, 63,192 

:  REM 

CALL  SCREG 

12350 

DATA 

201 

:  REM 

RET 

12360 

REM  < 

<<<<  SUBROUTINE 

SCREG  (21  BYTES  LONG 

12370 

DATA 

243 

:  REM 

DI 

i23eo 

DATA 

62,0 

:  REM 

LD  A,  00H 

12390 

DATA 

211,0 

:  REM 

OUT  00H 

12400 

DATA 

50,8,239 

:  REM 

LD  ( EFOSH ) , A 

12410 

DATA 

121 

:  REM 

LD  A,  C 

12420 

DATA 

50,0,42 

:  REM 

LD  (2A00H) , A 

12430 

DATA 

62,  1 

:  REM 

LD  A , 0 1 H 

12440 

DATA 

211,1 

:  REM 

OUT  0 1 H 

12450 

DATA 

50,8. 239 

:  REM 

LD  (EFOSH) , A 

1 2460 

DATA 

251 

:  REM 

El 

12470 

DATA 

201 

:  REM 

RET 

1 2480 

REM  < 

<<<<  SUBROUTINE 

RSIN 

(27  BYTES  LONG) 

12490 

DATA 

243 

:  REM 

DI 

12500 

DATA 

62,0 

:  REM 

LD  A , 00H 

12510 

DATm 

211,0 

:  REM 

OUT  OOH 

12520 

DATA 

50,8,239 

:  REM 

LD  (EF08H) , A 

1 2530 

DATA 

58,0, 42 

:  REM 

LD  A, ( 2A00H ) 

12540 

DATA 

31 

:  REM 

RRA 

12550 

DATA 

48,250 

:  REM 

JR  NC, TEST 

12560 

DATA 

58, 1,42 

:  REM 

LD  A, (2A01H) 

12570 

DATA 

79 

:  REM 

LD  C ,  A 

12580 

DATA 

62,  1 

:  REM 

LD  A, 01H 

12590 

DATA 

211,1 

:  REM 

OUT  01H 

12600 

DATA 

50,8,239 

:  REM 

LD  (EF08H) , A 

12610 

DATA 

251 

:  REM 

El 

12620 

DATA 

201 

:  REM 

RET 

12630 

DATA 

o 

‘M 

o 

o 

Ml 

o 

o 

•M 

o 

'M 

o 

o 

12640 

END 

APPENDIX  F 


FORMULAS  FOR  RMS-GRADIENT  METHOD 


For  brevity,  let  RMS,  our  Figure  of  Merit,  be  called  R, 
and  let  the  control -def 1 ecti on  array  be  called  X;  i,e.,  X  is 
vector  of  M  components,  X-  (i  =  1,2,...,M). 

When  we  make  M  runs,  measure  the  effects  on  RMS  of  the  M 
control  increments,  and  divide  these  effects  by  the  amounts  o 
the  control  increments,  we  have  determined  the  gradient, 
of  R,  viz.  grad  F:,  in  the  M-di  mensi  onal  X-space.  Its 

components  are  3f?/<9X  *  •  We  know  that  its  direction  in 

u 

X-space  is  the  direction  of  most  rapid  change  of  R. 

Let  S.i  denote  an  increment  of  X  in  X-space  (i.e.,  an 
array  of  control  increments).  The  general  formula  for  the 
effect  on  R  is 

S  R  =  &X  .  grad  R  (FI) 

If  £x  is  a  unit  vector  in  the  direction  of  grad  R,  we 


have 


<$X  = 


grad  R  /  Jgrad  R! 


Sr  =  grad  R  .  grad  R  /  Jgrad  R! 


=  !grad  R! 


Suppose  we  went  an  increment  of  R  equal  to  — kR,  where  k 
is  arbitrary  (a  relaxation  constant,  say),  and  R  is  the 
present  value  of  RMS.  We  want  Sr  to  be  -kR/ Igrad  R!  times  the 
vale  in  Eq. (F3) ,  so  we  must  use  a  £x  equal  not  to  the  value  in 
Eq. (F2)  but  equal  to 

-kR/!gradR!  times  grad  R  /  Igrad  R! 


which  is 


or 


-kR  grad  R  /  Igrad  R ! 

/£,  3R  W 

-kR  grad  R /2  (  g* .  ) 

.i-l  J 


CF4) 

(F5) 


This  is  a  control -i ncrement  array  whose  members  are 


APPENDIX  6 


s 

I 

& 

I 

£ 


I 

I 

K 

£ 

1 


§ 


$ 

i 

t  *, 

¥  * 

4  V 

* 


FORMULAS  FOR  "DOWELL’S  METHOD" 

Following  the  notation  used  earlier  in  this  report,  viz. 
in  the  section  entitled  "PROCEDURES",  let  the  control  matrix 
of  the  component  f  (which  denotes  Ve>  be  called  C.  Let  the 
analogous  control  matrix  of  component  g  (which  denotes  )  be 
called  D. 

Then  the  estimated  effects  of  any  array  of  control- 


setting  increments. 

say 

&X,  are 

£f  =  C  * 

(Gl) 

-  D  * 

8x 

(G2) 

and  £(Df)  =  (new  Df )  -  (old  Df > 

=  f  Cg  +  <£ g  1  -  f  -  Sf 

-  fCg]  +  f  (G3) 

In  the  same  local-linear  approximation  that  has  been  used 
earlier,  this  becomes 

&<Df )  *  f  C  5  gU  -  (ff  <B4> 

Here  we  return  to  our  more  descriptive  notation,  as  in 
Equation  (5),  and  substitute  for  $ f  and  8 9  from  Equations 
(Gl)  and  (G2)  above: 

<£(  Df )  or  £\DVe)  =  E*Binv*D*  §X-C*  &X  (G5) 
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—  ( iE  $  B  in  v  *  D  —  C)  ^  X 


Therefore,  to  achieve  a  change  of  DV  equal  to  -kDV  , 
say,  we  should  introduce  the  following  array  of  control- 
setting  increments: 

Sx  or  £<C.S.)  =  (C  -  E  *  B  inV  *  D)“1  *  kDVe  (G7 


This  is  the  counterpart  of  Equation  (6)  when  the 
procedure  we  call  “Dowell’s"  is  used.  It  is  clear  that  our 
regular  procedure  (Equation  (6)  >  results  from  neglecting  the 
second  term  inside  the  parentheses  in  Equation  (S7) .  It  is  a 
"secondary"  term,  probably  smaller  than  the  first  term  in 
low-speed  flow.  Its  neglect  might  be  expected  to  slaw  down  the 
convergence  of  the  iterative  process. 
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Matrices  Comparison 
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Figure  5.6 

Matrices  Comparison 
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Matrices  Comparison 
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Figure  5.10 

Matrices  Comparison 
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Matrix  for  Asymmetric  Run 
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Figure  5.24 

Matrix  for  Asymmetric  Runs 
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Matrix  for  Asymmetric  Runs 
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Matrix  for  Asymmetric  Runs 
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RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  4.S2  W  =  3.02 

RESULTANT  =  5.68795  at  32.0718  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  34 

RUN  No.  1  DATE  220285 


CONTROL  SETTINGS 


MATCHING  DISCREPANCY 


Lt.  Control  setting 

RC  Number 


Field  DVe 

Point 


10.6 
8.73 
7.65 
11.51 
11.39 
4.6 
6 . 36 
7.11 
6.21 
8.07 
34 
1 

220285 

4.82 

3.02 


1 

-1.0735 

2 

. 107226 

3 

.417379 

4 

. 436232 

5 

-.111687 

6 

-.919112 

7 

-.434061 

8 

— . 465648 

9 

-.387998 

10 

. 0572258 

11 

-. 133002 

12 

-.0610143 

13 

. 266227 

14 

. 157678 

15 

. 0446806 

16 


0 


16 


.0462749 


RMS 


RMS  ■  .437749 
percent  =  7.69608 


N 

* 

STREAM  VECTOR:  U  =  4.82  W 

=  3.02 

RESULTANT  = 

5.68795  at  32. 

0718  degrees 

FLAPS:  Position  4,  50  p.s 

.  i  . 

EXPERIMENT  No.  34 

* 

RUN 

No. 

2  DATE  220295 

Iteration  from  Exp 

.34  Run  1  V 

CONTROL 

SETTINGS 

MATCHING  DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

Point 

1 

8.95 

1 

-.807276 

2 

8.85 

-. 190745 

3 

8.11 

3 

. 427777 

4 

11.53 

4 

. 500862 

3 

11.52 

5 

. 126891 

6 

4.03 

6 

-.3681 

7 

6.68 

7 

-.324096 

8 

6.98 

8 

-. 441587 

9 

6.07 

9 

-.613425 

10 

8.04 

10 

-.0880244 

1 1 

34 

1 1 

-. 118548 

12 

r> 

12 

-.0139762 

13 

220285 

13 

.20721 

14 

4.82 

14 

. 172059 

15 

3.02 

15 

.0427341 

16 

0 

16 

. 0790849 

RMS  =  .358162 

RMS  percent  =  6.29685 

'%5 
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RESULTS  OF  ITERATION 
******************** 


I 


I 


STREAM  VECTOR: 


I 


FLAPS:  Position 

EXPERIMENT  No. 


RUN  No. 


DATE 


£ 


CONTROL  SETTINGS 


»  * 

Control 

setting 

s 

Number 

1 

9 

1 

<n> 

8.93 

■ 

3 

8.62 

4 

11.47 

k", 

5 

11.57 

'*• 

6 

■j .  66 

7 

7.02 

| 

8 

6.92 

* 

9 

5.96 

■ 

10 

8.02 

11 

34 

V 

12 

3 

i; 

13 

220285 

14 

4.82 

!* 

15 

3.02 

vl 

16 

0 

I 


I 


U  =  4. 82  W  «  3.02 

RESULTANT  =  5.68795  at  32.0718  degrees 
4,  50  p . s. i . 

34 


220285  Iteration  from  Exp. 34  Run  2 


MATCHING  DISCREPANCY 
DVe 


— . 865358 
-.211643 
.362402 
. 432963 
. 0485407 
-.270934 
-.3958 
-.375062 
-.486957 
.0140088 
-. 121978 
-.0465122 
. 32857 
. 179963 
.0700258 
.0586127 

RMS  =  .343563 
RMS  percent  =  6.04019 


Field-' 
Poi  nt 
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4 

5 

6 

7 

8 

9 

10 
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12 

13 

14 

15 

16 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  4.82  W  =  3.02 


RESULTANT  =  5.68795 


FLAPS:  Position  4,  50  p.s.i 


EXPERIMENT  No.  34 


at  32.0718  degrees 


RUN  No. 


DATE 


10385 


CONTROL  SETTINGS 


Control 

Number 


setting 


Iteration  from  Exp. 34  Run  4 


MATCHING  DISCREPANCY 

Field  DVe 

Point 


r?  1 


1 

9 

1 

-.919696 

Jm 

8.99 

o 

.0182368 

3 

10 

.217487 

4 

11.37 

4 

. 172068 

5 

11.63 

5 

-.0710959 

6 

2.93 

6 

.0921686 

7 

7.62 

7 

-.264858 

8 

6.8 

8 

-. 175738 

9 

5.82 

9 

-.368971 

10 

7.93 

10 

. 0776987 

11 

34 

11 

-. 143279 

12 

5 

12 

-. 121887 

13 

10385 

13 

.263451 

14 

4.82 

14 

. 13062 

15 

3.02 

15 

.0722919 

16 

0 

16 

. 078866 

RMS  *  .286481 
RMS  percent  =  5.03662 
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>•  “ 


«  .%  H "«  ,,*■  '■>  '»  '«  »'■  i,*>>  y|| 

« fjC oV*  - "  -  “  of 1  o  otvy\o  o?T 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  4.82  W  *  3.02 


RESULTANT  =  5.68795 


32.0718  degrees 


FLAPS:  Position  4,  50  p.s.i 


EXPERIMENT  No.  34 


RUN  No. 


DATE 


10385 


Iteration  from  Exp. 34  Run  5 


CONTROL  SETTINGS 


MATCHING  DISCREPANCY 


f'.  Control 
Number 


setting 


Field 

Point 


9 

9.  13 
10.35 
11.42 
11.75 

2.79 
7.65 

6.79 
5.78 
7.9 
34 

6 

10385 

4.82 

3.02 

O 


-.637391 
.19635 
—  .01  ■-•7 -j4o 
-.0974589 
-.256909 
. 0442748 
-.216988 
-. 172237 
-.260444 
_  . 26443 _ 
-. 124024 
-. 126574 
.284325 
.  0654967 
. 10503 
. 10085 


RMS  =  .233815 
RMS  percent  =  4.1107 


L‘ 


RESULTS  QF  ITERATION 
*****#**#*:*#*#*##*** 


STREAM  VECTOR:  U  =4.85  W  =  3 

RESULTANT  =  5.70285 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  34 


at  31.7415  degrees 


RUN  No. 


DATE 


10385 


Control 

Number 


CONTROL  SETTINGS 
1  setting 


2 

9.  13 

n 

4a 

£ 

sj 

4 

10.35 

11.42 

3 

4 

5 

11.75 

5 

6 

2.79 

6 

1 

7 

7.65 

7 

8 

6.79 

8 

9 

5.78 

9 

*" 

%  _ 

10 

7.9 

10 

$ 

11 

34 

1 1 

12 

6 

12 

Ml 

13 

1 0385 

13 

14 

4.85 

14 

15 

3.00 

15 

16 

0 

16 

Results  of  Search  for  best  RMS  % 

MATCHING  DISCREPANCY 

Field  DVe 

Point 


-.558286 
. 134801 
— . 0326639 
-. 107459 
-.257705 
.0474156 
-.270763 
-.245989 
-.346432 
. 17929 
-. 117779 
-. 123709 
. 29088 1 
. 0746809 
. 107499 
. 107334 


RMS  =  .229068 
RMS  percent  =  4.01673 


RESULTS  OF  ITERATION 
******* ******* ****** 


STREAM  VECTOR:  U  =  4.82  W  =  2.96 

RESULTANT  =  5.65632  at  31.5567  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  34 


RUN  No.  7 


date  10385  Results  of  Search  for  best  RMS  %  j 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

Poi  nt 

1 

9 

1 

-.589163 

2 

9.  12 

. 0606863 

T 

■-» 

1 0 . 38 

O 

-.0391496 

4 

11.5 

4 

-. 102472 

5 

11.78 

5 

-.237822 

6 

2.  72 

6 

-. 128468 

7 

7.51 

7 

-.309947 

8 

6.81 

8 

-.230426 

9 

5.74 

9 

-.367119 

10 

7.84 

10 

. 0758724 

1  1 

34 

1 1 

-. 141251 

12 

7 

12 

-. 1 18173 

13 

10385 

13 

.21104 

14 

4.82 

14 

. 130005 

15 

2.96 

15 

. 0465975 

16 

0 

16 

. 132739 

RMS  =  .229229 
RMS  percent  =  4.05261 
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RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  5.4  W  «  4.33 

RESULTANT  =  6.92163  at  38.7273  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  62 


RUN  No.  1  DATE  80785 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

Poi  nt 

1 

18 

1 

-.846241 

2 

9.99 

2 

. 299997 

3 

9.73 

3 

.246764 

4 

11.67 

4 

. 489093 

5 

12.43 

5 

-.405422 

6 

3.28 

6 

-.971464 

7 

10.29 

7 

.0179687 

8 

7.73 

8 

-1.13102 

9 

6.69 

9 

-.690904 

10 

7.71 

10 

-. 147961 

11 

62 

11 

-.270863 

12 

1 

12 

-. 133211 

13 

80785 

13 

. 0040642 

14 

5.4 

14 

-.0153173 

15 

4.33 

15 

. 0629309 

16 

0 

'  16 

-. 127118 

RMS  =  .506494 
RMS  percent  =  7.31756 


i 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  *  5.4  W  = 

RESULTANT  = 
FLAPS:  Position  4,  50  p.i 
EXPERIMENT  No.  62 


=  4.33 
=  6. 9216v 
>.  i  . 


at  38.7273  degrees 


^  RUN  No . 

2  DATE290785 

Iteration  from 

Exp. 62  Run  1 

i 

■  CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

/•  Control 
%■’  Number 

setting 

Field 

Point 

DVe 

1 

18.5 

1 

-. 194106 

■  2 

10.  17 

2 

-.252291 

"  3 

9.62 

3 

-.0489209 

4 

11.79 

4 

. 174577 

«  5 

12.46 

5 

— . 607667 

w  6 

3.08 

6 

-1.07875 

7 

10.34 

7 

-.0461287 

■  8 

7.75 

8 

-1. 18553 

B  9 

6.63 

9 

-.700273 

10 

7.88 

10 

-.236071 

11 

62 

11 

-.311262 

12 

2 

12 

-. 155531 

'v  13 

290785 

13 

. 137144 

14 

5.4 

14 

. 0296507 

V  15 

4.33 

15 

.0425143 

!v  16 

0 

16 

-.0853491 

RMS  =  .48546 
RMS  percent  =  7.0i: 


i66 


k  =  .15 


f 


b 

4 


V 


v 

V 

V' 

v' 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  5.6  W  =  4.49 

RESULTANT  =  7.17775  at  33.725  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  62 


RUN  No. 

2  DATE  290785 

Results  of 

Search  for  best 

CONTROL 

SETTINGS 

MATCHING  DISCREPANCY 

Control 

setting 

Fi  el  d 

DVe 

Number 

Point 

1 

18.5 

1 

-.699503 

2 

10.  17 

O 

.370143 

3 

9.62 

T 

. 358727 

4 

11.79 

4 

. 722863 

5 

12.46 

5 

.0242482 

6 

3.08 

6 

-.455729 

7 

10.34 

7 

. 12921 

8 

7.75 

8 

-1.06226 

9 

6.63 

9 

— . 636305 

10 

7.88 

10 

-.201586 

11 

62 

11 

-.248048 

12 

2 

12 

-. 139338 

13 

290785 

13 

. 180322 

14 

5.6 

14 

.0693018 

15 

4.49 

15 

.0618856 

16 

0 

16 

-.0415922 

RMS  =  .44723 
RMS  percent  =  6.23078 


RESULTS  OF  ITERATION 
*** ###*#*** ********* 


STREAM  VECTOR:  U  *  5.6  W  *  4.49 

RESULTANT  =  7.17775  at  33.725  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  62 


run  No.  3  date 310 785  Iteration  from  Exp. 62  Run  2 


CONTROL  SETTINGS  MATCHING  DISCREPANCY 


Control 

Number 

setting 

Field 

Point 

DVe 

1 

19.25 

1 

.21428 

2 

10.24 

2 

-.0509019 

3 

9.5 

3 

.0256528 

4 

11.82 

4 

.348692 

5 

12.38 

5 

-.385988 

6 

2.96 

6 

-.717946 

7 

10.29 

7 

.0746186 

8 

7.78 

8 

-1.03706 

9 

6.59 

9 

-.603746 

10 

8.02 

10 

-.206904 

11 

62 

1 1 

— • 25o368 

12 

3 

12 

-.0177063 

13 

310785 

13 

. 157974 

14 

5.6 

14 

-.0996144 

IS 

4.49 

15 

.0607711 

16 

0 

16 

-. 153824 

RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  5.7  W  *  4.57 

RESULTANT  =  7.30581  at  38.7239  degress 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  62 


run  No.  3  date  310785  Results  of  Search  for  best  RMS  % 


CONTROL  SETTINGS 

Control  setting 

Number 


MATCHING  DISCREPANCY 

Field  DVe 

Poi  nt 


1 

19.25 

1 

-. 0377684 

2 

10.24 

2 

. 259664 

3 

9.5 

-r 

*-» 

.229141 

4 

11.82 

4 

. 622447 

5 

12.38 

5 

-.0704317 

6 

2.96 

6 

-.406814 

7 

10.29 

7 

. 161951 

8 

7.78 

8 

-.975814 

9 

6.59 

9 

-.572163 

10 

8.02 

10 

-. 190041 

11 

62 

1 1 

-.222274 

12 

3 

12 

-9.60817E 

13 

310785 

13 

. 179563 

14 

5.7 

14 

-.0797752 

15 

4.57 

15 

.0704551 

16 

0 

16 

-. 131946 

RMS  =  .365005 
RMS  percent  =  4.9961 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  5.7  W 

RESULTANT 
FLAPS:  Position  4,  50  p 
EXPERIMENT  No.  62 


RUN  No. 


DATE  310785 


k 


Control 
k  Number 

1 
2 

3 

4 

5 

6 
7 

e 

9 

10 
n 
12 

13 

14 

15 

16 


CONTROL  SETTINGS 
setting 


20 

10.38 
9.87 
11.8 

12.38 
2.85 
10.  16 
7.8 

6.56 
8 . 05 
62 

4 

310785 

5.7 

4.57 
0 


=  4.57 

=  7.30581  at  38.7239  degrees 
.  s.  i  . 


Iteration  from  Exp. 62  Run  3  k  =  .15 


MATCHING  DISCREPANCY 


Fiel  d 

Poi  nt 

DVe 

1 

. 433659 

9 

.25014 

T 

w 

. 053035 

4 

. 476066 

5 

-. 194941 

6 

187924 

7 

. 264955 

8 

-.929725 

9 

-.69516 

10 

-.249723 

11 

-.204999 

12 

.0421891 

13 

. 159813 

14 

.0144481 

15 

.0692455 

16 

-.0785838 

RMS  =  .363501 
RMS  percent  =  4.9755 
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RESULTS  OF  ITERATION 
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STREAM  VECTOR:  U  *=  5.7  W  =  4.57 

RESULTANT  =  7.30531  at  33.7239  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  62 


RUN  No. 

5  DATE  20885 

Iteration 

from  Exp. 62  Run  4 

CONTROL 

SETTINGS 

MATCHING  DISCREPANCY 

Control 

Number 

setting 

Field 

Point 

DVe 

1 

20.9 

1 

-.306789 

2 

10.54 

2 

.  34o68v> 

3 

10.43 

O 

.096802 

4 

11.75 

4 

. 352709 

5 

12.4 

5 

-.291663 

6 

2.8 

6 

-.590209 

7 

9.91 

7 

.247234 

8 

7.81 

8 

-.404249 

9 

6.52 

9 

-. 129743 

10 

8.04 

10 

. 122896 

11 

62 

11 

-.20693 

12 

5 

12 

-.037821 

13 

20885 

13 

. 156659 

14 

5.7 

14 

.0322912 

15 

4.57 

15 

.0542122 

16 

0 

16 

— . 056 1 554 

RMS  =  .263609 
RMS  percent  =  3.6082 


RESULTS  OF  ITERATION 
******************** 


S 


& 

v* 

& 


STREAM  VECTOR:  U  =  5.7  W  =  4.57 

RESULTANT  =  7.30581  at  38.7239  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  62 


RUN  No.  6  DATE  20885  Iteration  from  Exp. 62  Run  5 


£ 


CONTROL  SETTINGS 


MATCHING  DISCREPANCY 


1 

8 

I 

jSS 

T 


Control 

Number 

setting 

Field 

Poi  nt 

DVe 

1 

20.9 

1 

. 142659 

10.4? 

2 

. 55993 

3 

10. 33 

3 

-. 149465 

4 

11.73 

4 

.0348871 

5 

12.27 

5 

-.569635 

6 

2.64 

6 

—  •  33520 1 

7 

9.86 

7 

. 366732 

8 

7.82 

8 

-.519915 

9 

6.52 

9 

-. 154074 

10 

8.09 

10 

. 131184 

11 

62 

11 

-.307006 

12 

6 

12 

. 043^078 

13 

20885 

13 

. 183093 

14 

5.7 

14 

.  1762 

15 

4.57 

15 

. 103248 

16 

0 

16 

-. 116921 

f: 

•  * 

U- 


RMS  =  .29843? 

RMS  percent  =  4.08495 


-V 


.v 

.V 


« 
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RESULTS  OF  ITERATION 


STREAM  VECTOR:  U  =  5.7  W  =  4.57 

RESULTANT  =  7.30581  at  38.7239  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  62 


RUN  No.  8  date  90885  Iteration  from  Exp. 62  Run  7 


CONTROL  SETTINGS 

Control  setting 

Number 


1 

n 

3 

4 

5 

6 

7 

8 

9 

21.3646 
10. 4154 
9.93014 
11.7372 
12. 1124 
2.55138 
9.79566 
7.80316 
6.48743 

10 

8.24247 

11 

62 

12 

8 

13 

90885 

14 

5.7 

15 

4.57 

16 

0 

MATCHING  DISCREPANCY 

Field  DVe 

Point 


1 

-.650964 

-.038782 

3 

.219081 

4 

.414383 

5 

-. 198584 

6 

.216497 

7 

-.0936011 

8 

-. 546989 

9 

-. 194427 

10 

.0518292 

11 

-.28317 

12 

6.06763E 

13 

. 208772 

14 

. 1 1 335 

15 

.0529896 

16 

-. 138106 

RMS  =  .278144 
RMS  percent  =  3.80715 


RESULTS  OF  ITERATION 
**^***************** 


STREAM  VECTOR:  U  =  5.6  W  =  4.55 

RESULTANT  =  7.21544  at  39.0967  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  62 


run  No.  s  date  90885  Results  of  Search  for  best  RMS  % 


CONTROL  SETTINGS  MATCHING  DISCREPANCY 


Control 

setting 

Field 

DVe 

Number 

Point 

1 

21.3646 

1 

-.609897 

10.4154 

n 

-. 138367 

sj 

9.93014 

. 123651 

4 

1 1 . 7372 

4 

. 266254 

5 

12. 1124 

5 

— . 383963 

6 

2  •  55 1  o3 

6 

. 0283628 

7 

9.79566 

7 

-.0718762 

8 

7.80316 

8 

-. 482612 

9 

6.48743 

9 

-.0958326 

10 

8.24247 

10 

. 157964 

11 

62 

11 

•  1 0355 

12 

8 

12 

-2.62515E 

13 

90885 

13 

. 187074 

14 

5.6 

14 

. 089 1 032 

15 

4 . 55 

15 

. 0439004 

16 

0 

16 

-. 159876 

RMS  =  .2579 

37 

RMS  percent  = 

3.57479 

57479 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =5.7  W  =  4.57 

RESULTANT  =  7.30581  at  38.7239  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  62 


run  No.  9  date  260885  Iteration  from  Exp. 62  Run  8 


CONTROL  SETTINGS 


MATCHING  DISCREPANCY 


Control  setting 

Number 


Field  DVe 

Point 


1 

20.9 

1 

-.269358 

10.  34 

'“S 

-.534495 

3 

9.44 

3 

-. 143168 

4 

11.75 

4 

-.0739921 

5 

11.98 

5 

-.679997 

6 

2.48 

6 

. 358658 

7 

10.03 

7 

-.323569 

8 

7.84 

8 

-.400527 

9 

6.48 

9 

7.98999E 

10 

8 . 36 

10 

.811223 

1 1 

62 

11 

-.345491 

12 

9 

12 

- . 0806786 

13 

260885 

1  *-> 

. 224805 

14 

5.7 

14 

-.265732 

15 

4.57 

15 

-. 147428 

16 

0 

16 

— .  1 63 1 65 

RMS  =  .370295 
RMS  percent  =  5.06849 


RESULTS  OF  ITERATION 
ttt ******  *********** 


STREAM  VECTOR:  U  =  5.35  W  =  4.6 

RESULTANT  =  7.44194  at  38.1317  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  62 


V. 

V. 


RUN  No. 


DATE 


26038: 


Results  of  Search  for  best  RMS  % 


CONTROL  SETTINGS 


Control 

Number 


§  3 

4 


6 

7 

8 
o 

10 
1 1 
12 

13 

14 

15 

16 


setting 


20.9 

10.34 

9.44 

11.75 

11.98 

2.  43 

10.03 

7.84 
6.48 
8.36 
62 

9 

260885 

5. 85 
4.6 
O 


MATCHING  DISCREPANCY 


Field 
Poi  nt 


4 


u 

6 

7 

8 

9 

10 
1 1 
12 

13 

14 

15 

16 


DVe 


-.330958 
-.385118 
-1 .52362E-03 
. 143201 
- . 40 1 928 
. 640859 
-.356157 
-.497094 
-. 139903 
. 65202 
-.304713 
-.0676394 
. 257351 
-.229361 
-.  133794 
-. 130511 


RMS  =  .346582 
RMS  percent  =  4.65715 


a 

n 


Vi 


a 


m 


,v 


.V 

jj 


m 


Figure  6.3 


RESULTS  OF  ITERATION 
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STREAM  VECTOR:  u  =  6  W  =  4 

RESULTANT  =7.2111  at  33.6926  degrees 

FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  76 


RUN  No.  1  DATE  60186 


CONTROL  SETTINGS 


MATCHING  DISCREPANCY 


Control  setting 

Number 


Field  DVe 

Point 


1 


4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

1 6 

17 

18 

19 

20 
21 


24 


26 

27 

28 

29 

30 

31 


20.6 

10.24 

11.6 

11.74 


1  1.93196 

2  -2.81352 

^  .  IcOvi 

4  .41151 


12.85 

3 

-.926918 

2.06 

6 

-.520968 

7.08 

7 

-.711694 

6.65 

8 

-.670019 

5.29 

9 

-.96393 

8.27 

10 

-.912625 

9.01 

11 

-.772983 

12.4 

12 

-1.03816 

12.  19 

13 

1. 10604 

2.49 

14 

1.54951 

5 . 93 

15 

1.46971 

6.55 

16 

1 . 53936 

6.  17 

17 

. 256473 

7.95 

18 

.410514 

76 

19 

.341187 

1 

20 

. 0094422 

60186 

21 

— . 273545 

6 

r>r> 

-.0404702 

4 

A>W 

1. 10328 

0 

24 

-.0112939 

0 

25 

. 0433099 

0 

26 

-.0353547 

0 

27 

- . 1 30003 

0 

28 

-.0182608 

0 

29 

. 0578872 

0 

30 

.104713 

0  ? 

31 

. 123459 

0 

32 

-.0573147 

RMS  =  1.07893 


RMS  percent  =  14.9623 


RESULTS 

OF 

ITERATION 

******:|^************ 

STREAM  VECTOR: 

U  = 

6  W  =  4 

RESULTANT  =  7. 

2111  at  33. 

6926  degrees 

FLAPS:  Position 

4, 

50  p . s . i . 

EXPERIMENT  No. 

76 

RUN 

No. 

2  DATE 

60186  Iteration  from 

Exp. 76  Run  1 

CONTROL 

SETTINGS 

MATCHING  DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

Point 

1 

22.6 

1 

1 . 87304 

10.29 

-2. 19969 

3 

12.22 

3 

-2.49464 

4 

11.6 

4 

-.505959 

S 

13.29 

5 

-.790615 

6 

1.64 

6 

-.349789 

7 

6.98 

7 

-.554302 

8 

6.62 

8 

-.304679 

9 

5.  13 

9 

-.783578 

10 

8. 33 

10 

-. 60555 

11 

8.  17 

11 

— . 255966 

12 

12.47 

12 

-.642662 

13 

11.97 

13 

.901547 

14 

2.99 

14 

1.32611 

15 

5.43 

15 

1. 15437 

16 

6.47 

16 

.7446 

17 

6.48 

17 

. 0755932 

18 

7.91 

18 

-. 101481 

19 

76 

19 

-.302102 

20 

r> 

20 

-.528102 

21 

60186 

21 

-.327858 

6 

r>r> 

-. 175263 

23 

4 

23 

.715771 

24 

0 

24 

-.0452944 

25 

0 

25 

.0918714 

26 

0 

26 

-. 139337 

27 

0 

27 

- . 1 63005 

28 

0 

28 

. 136013 

29 

0 

29 

.324128 

30 

0 

30 

. 069305 1 

31 

0 

31 

. 157549 

■j  2 

0 

32 

-. 101026 

RMS  =  .850511 

RMS 

percent  =  11 

.7945 

113 

1.--  L>  .v  „-V  i.%  A 

JV.% 
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.*>•  .*•  »>  /•  ; 

RESULTS  OF  ITERATION 
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STREAM  VECTOR:  U  *  6  W  =  4 

RESULTANT  =7.2111  at  33.6926  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  76 


run  No.  3  date  70186  Iteration  from  Exp. 76  Run  2 


CONTROL  SETTINGS 


MATCHING  DISCREPANCY 


Control 

Number 

setting 

Field 

Point 

DVe 

1 

24.5 

1 

.914148 

jL 

10.  26 

n 

4_ 

-1. 10769 

3 

12.59 

3 

-1.67753 

4 

11.49 

4 

-.702318 

5 

13.55 

5 

-.573965 

6 

1.31 

6 

. 0668495 

7 

6.97 

7 

-.0184752 

8 

6.55 

8 

. 100844 

9 

5.  14 

9 

-.395746 

10 

8.39 

10 

-.0466141 

11 

7.37 

11 

. 262995 

12 

12.48 

12 

-.217792 

13 

11.84 

13 

. 500967 

14 

3.46 

14 

. 655602 

15 

5.28 

15 

.0357181 

16 

6.49 

16 

-.265804 

17 

6.51 

17 

. 154772 

18 

7.93 

IS 

— . 0639054 

19 

76 

19 

-.350879 

20 

3 

20 

-.63788 

21 

80186 

21 

-.262463 

22 

6 

r>2. 

-.0394212 

23 

4 

23 

. 355884 

24 

0 

24 

-.0901826 

25 

0 

25 

. 0740584 

26 

0 

26 

-.223045 

27 

0 

27 

-. 196313 

28 

o 

28 

. 174619 

29 

0 

29 

.  156418 

30 

0 

30 

.0491077 

31 

0 

31 

.213751 

32 

0 

32 

-.074993 

RMS  =  .490798 
RMS  percent  =  6.80614 


RESULTS  OF  ITERATION 
:M****************** 


STREAM  VECTOR:  U  =  6  W  =  4 

RESULTANT  *  7.2111  at  33.6926  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  76 


run  No.  4  DATE  801S6  Iteration  from  Exp. 76  Run  3 


CONTROL  SETTINGS  MATCHING  DISCREPANCY 

Control  setting  Field  DVe 

Number  Point 


1 

25.6 

1 

1 16143 

o 

4- 

10.25 

2 

-1.53448 

3 

12.79 

3 

-1.70159 

4 

11.56 

4 

-1.28354 

5 

1  .  5c* 

5 

-.0632382 

6 

1.23 

6 

.693812 

7 

6.84 

7 

.596102 

8 

6.48 

8 

.748137 

9 

5.3 

9 

-.490838 

10 

8.34 

10 

-.208208 

11 

7.8 

11 

. 0649803 

12 

12.47 

12 

164201 

13 

11.86 

13 

. 734603 

14 

3.62 

14 

. 903849 

15 

5.  12 

15 

.257761 

16 

6.5 

16 

— .  43oc*88 

17 

6.45 

17 

-.266788 

18 

7.99 

18 

-.622334 

19 

76 

19 

-.756774 

20 

4 

20 

-1. 19644 

21 

80186 

21 

-. 167287 

22 

6 

22 

-. 154001 

23 

4 

23 

.226355 

24 

0 

24 

-.0389358 

25 

0 

25 

.0254461 

26 

0 

26 

-. 125643 

27 

0 

27 

-. 127765 

28 

0 

28 

.27181 

29 

0 

29 

. 362588 

30 

0 

30 

.0913814 

31 

0 

31 

.216355 

32 

0 

32 

-. 100583 

RMS  =  .640025 
RMS  percent  =  S. 87554 


k 


.15 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  6.12  W  =  4.03 

RESULTANT  =  7.32771  at  33.3672  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  76 


RUN  No. 

4  DATE 

80186 

Results  of 

Search  for  best 

CONTROL 

SETTINGS 

MATCHING  DISCREPANCY 

Control 

setting 

Fi  el  d 

DVe 

Number 

Point 

1 

25.6 

1 

.699647 

2 

10.25 

2 

-.504814 

3 

12.79 

3 

-.885841 

4 

11.56 

4 

-.253876 

5 

13.53 

5 

-.582857 

6 

1.23 

6 

. 128292 

7 

6.84 

7 

.0759266 

8 

6.48 

8 

.680393 

9 

j  ■  -j 

9 

. 1 4006 

10 

8.34 

10 

. 488357 

11 

7.8 

11 

.719259 

12 

12.47 

12 

. Oo 

13 

11.86 

13 

.214993 

14 

3.  62 

14 

. 338338 

15 

5.  12 

15 

-.262403 

16 

6.5 

16 

-.501118 

17 

6.45 

17 

.364122 

IS 

7.99 

18 

.0742512 

19 

76 

19 

-. 102472 

20 

4 

20 

-.99656 

21 

80186 

21 

-.342978 

22 

6.  12 

22 

-.161 933 

23 

4.03 

23 

.117049 

24 

0 

24 

-.219238 

oe 

0 

25 

.0181102 

26 

0 

26 

-.234883 

27 

0 

27 

-.303452 

28 

o 

28 

. 263881 

29 

0 

29 

.253281 

30 

0 

30 

-.0889087 

31 

0 

31 

.209017 

32 

0 

32 

-.209822 

RMS  =  .410796 
RMS  percent  =  5. 


60607 


RESULTS  OF  ITERATION 
************** ****** 


STREAM  VECTOR:  U  =  6  W  =  4 

RESULTANT  =7.2111  at  33.6926  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  76 


run  No.  5  DATE  100186  Iteration  from  Exp. 76  Run  4 


CONTROL  SETTINGS  MATCHING  DISCREPANCY 


Control 

Number 

setting 

Field 

Point 

DVe 

1 

26.4 

1 

1 . 24042 

2 

10.  15 

2 

-1.78095 

3 

12.  14 

3 

-2.58105 

4 

11.77 

4 

-1.19607 

5 

13.29 

5 

-.526091 

6 

1.42 

6 

.0679979 

7 

6.95 

7 

-.0840388 

8 

6.61 

8 

-.0303798 

9 

5.44 

9 

-.609038 

10 

8.27 

10 

-. 403283 

11 

7.32 

11 

-.222102 

12 

12.31 

12 

-.337194 

13 

11.88 

13 

.968516 

14 

3.79 

14 

1.23392 

15 

5.34 

15 

.595147 

16 

6.51 

16 

-.0963131 

17 

6.22 

17 

-.0674332 

18 

8.  13 

18 

-. 230224 

19 

76 

19 

-.362096 

20 

5 

20 

-.719651 

21 

100186 

21 

-.214692 

22 

6 

Tin 

-. 141538 

23 

4 

23 

.351733 

24 

0 

24 

-.0405744 

25 

0 

25 

. 0275096 

26 

0 

26 

-. 154647 

27 

0 

27 

- . 07536 1 8 

28 

0 

28 

. 233966 

29 

0  ? 

29 

.529115 

30 

0 

30 

. 13101 

31 

0 

31 

. 143975 

32 

0 

32 

RMS  =  .74834 

-.0808527 

RMS  percent  =  10.3776 


RESULTS  OF  ITERATION 
tt*ttt**tt*tt*t***tt 


STREAM  VECTOR:  U  =  6.12  W  =  4.03 

RESULTANT  =  7.32771  at  33.3672  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  76 


RUN  No.  6  DATE  100186  Iteration  from  Exp. 76  Run  4 


CONTROL  SETTINGS  MATCHING  DISCREPANCY 


Control 

setting 

Field 

DVe 

Number 

Point 

1 

26.5 

1 

-.370616 

o 

10.37 

*■> 

-.57602 

13.48 

T 

>-* 

-.769646 

4 

11.66 

4 

-.70944 

5 

13.53 

5 

-.228761 

6 

1.  16 

6 

.366531 

7 

6.62 

7 

. 196559 

8 

6.39 

8 

. 575802 

9 

5.58 

9 

— . 0608375 

10 

8.27 

10 

. 102511 

11 

7.3 

11 

. 169874 

12 

12.54 

12 

-.212846 

13 

11.9 

13 

.0750564 

14 

3.68 

14 

. 13083 

15 

4.64 

15 

-.529618 

16 

6.46 

16 

-.787823 

17 

6.52 

.  17 

. 3364 1 3 

18 

8.  1 

18 

-.0313194 

19 

76 

19 

-.234776 

20 

6 

20 

-.707018 

21 

100186 

21 

-.266289 

22 

6.  12 

22 

-.0162528 

23 

4.03 

23 

. 174903 

24 

0 

24 

-.0823276 

25 

0 

25 

.0402406 

26 

0 

26 

-.228064 

27 

0 

27 

-.285685 

28 

0 

28 

.217149 

29 

0 

29 

. 234485 

30 

0 

30 

-.0400586 

31 

0 

‘ 

31 

.222953 

32 

0 

32 

-. 1 19611 

RMS  =  .362916 
RMS  percent  =  4.95265 


RESULTS  OF  ITERATION 
*^****************** 


STREAM  VECTOR:  U  =  6.15  W  =  4.05 


RESULTANT  =  7.36376 


FLAPS:  Position  4,  50  p.s.i. 


EXPERIMENT  No.  76 


at  33.3688  degrees 


RUN  No. 


date  100186  Results  of  Search  for  best  RMS  % 


CONTROL  SETTINGS 


MATCHING  DISCREPANCY 


Control 

Number 


setting 


Field 
Poi  nt 


1 

26.5 

10.37 

-T 

W 

13.48 

4 

11.66 

5 

13.53 

6 

1.16 

7 

6.62 

8 

6.39 

9 

5.58 

10 

8.27 

1 1 

7.3 

12 

12.54 

13 

11.9 

14 

3.68 

15 

4.64 

16 

6.  46 

17 

6.52 

18 

8.  1 

19 

76 

20 

6 

21 

100186 

6.  15 

4.05 

24 

0 

25 

0 

26 

0 

27 

0 

28 

0 

29 

0 

30 

0 

31 

0 

32 

0 

1 

-.211227 

•n, 

-.274041 

3 

-.610268 

4 

-.407462 

5 

-.335481 

6 

. 252456 

7 

. 0944578 

8 

. 586395 

9 

. 120074 

10 

. 303959 

11 

.361389 

12 

-. 135353 

13 

— • 0316613 

14 

.0167572 

15 

-.631716 

16 

-.777227 

17 

.517328 

18 

. 170135 

19 

-.0432557 

20 

-.629518 

21 

-.309253 

22 

-.01836 

23 

. 147562 

24 

-. 128363 

25 

. 0385306 

26 

-.25536 

27 

-.328648 

28 

. 2 1 5043 

29 

.207144 

30 

- . 0860908 

31 

.221242 

RMS  =  . 336255 
RMS  percent  =  4.56634 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  6.  12  W  =  4.03 

RESULTANT  =  7.32771  at  33.3672  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  76 


RUN  No. 


date  130186  Iteration  from  Exp. 76  Run  6 


Control 

Number 


CONTROL  SETTINGS 
setting 


MATCHING  DISCREPANCY 


% 


Field 

Point 


DVe 


1 

26.8 

1 

. 348866 

n 

10.38 

2 

-.914142 

3 

13.32 

3 

-1.29477 

4 

11.86 

4 

-.753692 

5 

13.4 

5 

-.444226 

6 

1.27 

6 

.0901634 

7 

6.51 

7 

1. 10274E 

8 

6.43 

8 

. 10621 

9 

5.7 

9 

-. 191545 

10 

8.23 

10 

-.0494778 

11 

6.98 

11 

.0128038 

12 

12.54 

12 

-.0927146 

13 

11.97 

1 3 

. 326044 

14 

3.62 

14 

.359368 

15 

4.41 

15 

-.217405 

16 

6.49 

16 

-.827755 

17 

6.  46 

17 

•  2531 18 

18 

8.  18 

18 

-. 178262 

1? 

76 

19 

-.403909 

20 

7 

20 

— . 701 338 

21 

130186 

21 

-.245163 

22 

6.  12 

22 

-.0720376 

23 

4.03 

23 

.302214 

24 

0 

24 

-.0153799 

25 

0 

25 

. 0309835 

26 

0 

26 

-. 153599 

27 

0 

27 

-. 158814 

28 

o 

28 

. 383758 

29 

0 

29 

. 1 82627 

30 

0 

30 

-.011 058 1 

31 

0 

31 

. 143041 

32 

0 

32 

4. 03623E- 

m 


RMS  =  .418072 
RMS  percent  =  5.70536 


RESULTS  OF  ITERATION 

♦  lit*##:*******#**##### 


STREAM  VECTOR:  U  =  6. 16  W  =  4.09 

RESULTANT  =  7.39417  at  33.5S51  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  76 


RUN  No.  7  DATE  130186  Results  of  Search  for  best  RMS  % 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

Point 

1 

26.8 

1 

. 44256 

2 

10.  38 

o 

-.392686 

sJ 

13.  32 

3 

-1.20109 

4 

11.86 

4 

-.232235 

5 

13.4 

5 

-.524692 

6 

1.27 

6 

.0108778 

7 

6.51 

7 

-.06151 

8 

6 .  4o 

8 

. 193742 

9 

5.7 

9 

.111495 

10 

8.23 

10 

.291931 

11 

6. 98 

11 

. 342668 

12 

12.54 

12 

.0840168 

13 

11.97 

13 

.245582 

14 

3.62 

14 

. 280087 

15 

4.41 

15 

-.27902 

16 

6.49 

16 

-.740216 

17 

6.46 

17 

.556162 

18 

8.  18 

IS 

. 163155 

19 

76 

19 

-.0740362 

20 

7 

20 

-.524597 

21 

130186 

21 

-.299889 

22 

6.  16 

22 

-.075178 

23 

4.09 

.265722 

24 

0 

24 

-.0793182 

25 

0 

25 

. 0290346 

26 

0 

26 

-. 189955 

27 

0 

27 

-.213539 

28 

o 

28 

. 3806 1 9 

29 

0 

29 

. 146134 

30 

0 

30 

-.0749917 

31 

0 

31 

.141091 

0 

32 

-.0323191 

RMS  =  .362907 
RMS  percent  =  4.90801 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  6.12  W  =  4.03 


RESULTANT  =  7.32771 


at  33.3672  degrees 


FLAPS:  Position  4,  50  p.s.i. 


EXPERIMENT  No.  76 


RUN  No. 


DATE 


130186 


CONTROL  SETTINGS 


ontrol 

umber 


setting 


26.8 
10.38 
13.  32 
11.86 
13.4 

I. 27 
6.51 
6.  43 
5.7 
8.23 
6.98 
12.54 

II. 97 
3.62 
4.41 
6.  49 
6.46 
8.  18 
76 

8 

1 30 1 86 
6.  12 
4.03 
0 
0 
0 
0 
O 
0 
0 
o 
0 


Repeat  of  Exp.  76  Run  7 


MATCHING  DISCREPANCY 


Field 

Point 


.343417 
-.720713 
-1.15019 
-.693664 
-.465668 
. 062929 
-.0373104 
.296885 
-. 134171 
.024671 
.0885149 
-.0567751 
. 236645 
.239661 
— . 36527 6 
-.727437 
. 305904 
-. 1 14867 
-.339279 
-.673779 
-.263767 
-.0747336 
.288156 
— . 0366095 
.0288844 
-. 166727 
-. 183299 
.3777 
. 169858 
-.0326707 
. 140021 
-9. 6306 IE— O' 


RMS  =  .381093 
RMS  percent  =  5.20071 
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RESULTS  OF  ITERATION 
*ttt*t**tt ********** 

STREAM  VECTOR:  U  =  6.15  W  =  4.08 

RESULTANT  =  7.38031  at  33.5634  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  76 


RUN  No. 


8 


DATE  130186 


{ 

CONTROL  SETTINGS 

'J 

^  Control  setting 

Number 


i 

M 

k 

* 

I 


V 

V 


4 

5 

6 

7 

8 

9 

10 
1 1 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


24 

25 

26 

27 

28 

29 

30 

31 


26.8 
10.38 
1  3  •  --<2 
11.86 
13.4 

I. 27 
6.51 
6.43 
5.7 
8.23 
6.98 
12.54 

II. 97 
3.62 
4.41 
6.49 
6.46 
8.  18 
76 

8 

130186 
6.  15 
4.08 

0 

0 

0 

0 

0 

0 

0 

0 

0 


Results  of  Search  for  best  RMS  % 


MATCHING  DISCREPANCY 


Field  DVe 

Poi  nt 


1 

. 395866 

r> 

-.311801 

-1.09775 

4 

-.284746 

5 

-.516745 

6 

.0143855 

7 

-.07235 

8 

. 373544 

9 

. 102384 

10 

.291651 

1 1 

.347091 

12 

. 0867847 

13 

.  185571 

14 

.  191121 

15 

- . 4003 1 2 

16 

-.650772 

17 

.542462 

18 

. 152119 

19 

-.0806962 

20 

-.53021 1 

21 

-.304429 

-.0771385 

23 

.260781 

24 

-.0849473 

.0274723 

26 

-. 193988 

27 

— . 223959 

28 

.375297 

29 

. 142483 

30 

- . 08 1 005 

31 

. 138608 

-.0368916 

RMS  =  .348481 
RMS  percent  =  4.72177 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  6.12  W  =  4.03 

RESULTANT  =  7.32771 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  76 


33.3672  degrees 


RUN  No. 


date  150186  Iteration  from  Exp. 76  Run  8 


CONTROL  SETTINGS 


MATCHING  DISCREPANCY 


ontrol 

umber 

setting 

Field 

Poi  nt 

DVe 

1 

27.3 

1 

.265616 

n 

10.39 

•-y 

-.271352 

w 

13.42 

3 

-.83825 

4 

12.01 

4 

-.877068 

5 

13.36 

5 

-.378671 

6 

1.21 

6 

. 097093 

7 

6.4 

7 

-.0972522 

8 

6.42 

8 

. 438776 

9 

5.8 

9 

. 144892 

10 

8.  15 

10 

.317891 

11 

6.58 

11 

.414971 

12 

12.54 

12 

- 1 . 29689E-03 

13 

12.05 

13 

. 0419285 

14 

3.67 

14 

3. 1 6059E-03 

15 

4.  13 

15 

/  '-n  c* -r  T 

16 

6.5 

16 

-.74894 

17 

6.4 

17 

. 495 1 1 4 

18 

8.25 

18 

. 0256223 

1? 

76 

19 

-. 199397 

20 

9 

20 

-.640831 

21 

150186 

21 

-. 290732 

6.  12 

22 

-.229881 

n? 

4.03 

OT 

.15158 

24 

0 

24 

-.11 7303 

-L.O 

0 

C 

.0425142 

26 

0 

26 

-.31 1502 

27 

0 

27 

-.351 34 

28 

0 

28 

. 397899 

29 

0 

29 

.  17709? 

30 

0 

30 

-.0887928 

31 

0 

31 

.217527 

32 

0 

-. 171004 

RMS  =  .37855? 

RMS  percent  =  5.16613 
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RESULTS  OF  ITERATION 
******  ************** 


STREAM  VECTOR:  U  =  6. 12  W  =  4.06 


RESULTANT  =  7.34425 


FLAPS:  Position  4,  50  p.s.i 


EXPERIMENT  No.  76 


33.5627  dearees 


RUN  No. 


DATE 


150186 


CONTROL  SETTINGS 


Control 

Number 


setting 


27.3 

10.3? 

13.42 
12.01 
13.36 
1.21 

6.4 

6.42 
5.8 
8.  15 
6 . 58 
12.54 
12.05 
3.67 
4.  18 

6.5 
6.4 
8.25 
76 

9 

150186 
6.  12 
4.06 
0 
0 
0 
0 
O 
0 
0 
0 
0 


Results  of  Search  for  best  RMS  % 


MATCHING  DISCREPANCY 

Field  DVe 

Point 

1  . 158675 

2  -.164911 

3  -.945191 

4  -.770127 

5  - . 323027 

6  . 16263 

7  -.0301897 

8  . 504843 

9  . 200537 

10  .383423 

1 1  . 482033 

12  .0647706 

13  .0975729 

14  .0736939 

15  -.55816? 

16  -.68287 

17  .550758 

18  .0911554 

19  -.132334 

20  -.574763 

21  -.28843 

22  -.23017? 

23  .151545 

24  -.119606 

25  .042812 

26  -.311468 

27  -.349038 

28  . 397602 

29  .177064 

30  - . 09 1 095 

31  .217825 

32  -.170969 


RMS  =  .373751 
RMS  percent  =  5.08903 
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RESULTS  OF  ITERATION 
*t**t*t***t*t**tx*t* 

STREAM  VECTOR:  U  =  6  W  =  4 

RESULTANT  =  7.2111  at  33.6926  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  79 


RUN  No. 


DATE 


502S6 


Control 

Number 


CONTROL  SETTINGS 
setting 


MATCHING  DISCREPANCY 


5 

$ 


10 

11 

12 

13 

14 

15 

16 


16.3 
10.  IS 
10.  1 
12.  15 
12.52 
2 . 25 
7.66 
6. SI 
5.94 
8.  13 
79 
1 

50286 

6 

4 

0 


Field 

Point 


DVe 


1 

-1.27901 

o 

*■» 

.440601 

3 

. 558702 

4 

. 723035 

5 

. 245504 

6 

-.0582475 

7 

-.213245 

S 

— .  1 3538-. 

9 

—  .  1  4uO%j'U 

10 

-. 103974 

11 

-.249632 

12 

-.0795374 

13 

. 162493 

14 

. 168597 

15 

. 110476 

16 

7. 2855  IE 

RMS  =  .431241 
RMS  percent  =  5.9S024 


RESULTS  OF  ITERATION 
*  *  *  .****#*; *********** 


I 

£3 


STREAM  VECTOR:  U  =  6  W  =  4 

RESULTANT  =  7.2111  at  33.6926  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  79 


TF, 

iV 

RUN  No. 

2  DATE 

70286 

Iteration  from 

Exp. 79  Run  1 

!N- 

Is 

CONTROL 

SETTINGS 

MATCHING  DISCREPANCY 

1 

Control 

setting 

Field 

DVe 

L*\ 

Number 

Point 

' 

1 

14.74 

1 

-.800182 

L% 

2 

10.09 

-.544487 

| 

w 

9.56 

-T 

.366994 

■ 

4 

12.  1 

4 

. 620643 

5 

12.54 

5 

. 0430775 

£ 

6 

2.41 

6 

-.0582329 

■  , 

7 

8 .  05 

7 

™ .  o  1  o  1  o  9 

8 

6.91 

8 

— . 6156 

■ 

9 

5.86 

9 

-.517861 

I 

10 

8.  17 

10 

-.521491 

Or 

11 

79 

11 

— .  loo-j-42 

12 

2 

12 

-.0703102 

fc? 

13 

70286 

13 

. 102645 

» 

14 

6 

14 

. 160989 

15 

4 

15 

. 115465 

<5 

16 

0 

16 

-.0140474 

RMS  =  .412083 
RMS  percent  =  5.71457 
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RESULTS  OF  ITERATION 


STREAM  VECTOR:  u  =  6  W  =  4 

RESULTANT  =  7.2111 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  79 


at 


6926  degrees 


RUN  No. 


DATE 


702S6 


Iteration  from  Exp. 79  Run  2 


t' 


Control 

Number 


CONTROL  SETTINGS 
setting 


MATCHING  DISCREPANCY 


U 

6 

7 

8 

9 

10 
1 1 
12 

13 

14 

15 

16 


Fi  el  d 
Point 


DVe 


15.5 

10.51 

10.07 

12.25 

12.84 

2.29 

8.37 

7.01 

5.81 

8.28 

79 


1 

r> 


70286 

6 

4 

0 


6 

7 

8 

9 

10 
1 1 
12 

13 

14 

15 

16 


195676 
- . 530984 
•  1  o6^>o4 
.469373 
-.202079 
-1.  18665E— O- 
-.463691 
— . 659859 
-.529316 
-.516517 
-. 151159 
-.0470836 
. 1 46009 
. 121091 
. 135255 
-.0397562 


RMS  =  .342626 
RMS  percent  =  4.75137 
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RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  5.S  W  =  4. 085 


RESULTANT  =  7.09417  at  35.16  degrees 


FLAPS:  Position  4,  50  p.s.i 

EXPERIMENT  No.  79 


RUN  No. 


DATE 


70286 


CONTROL  SETTINGS 


Control 

Number 


setting 


10.51 
10.07 
12.25 
12.84 
2.29 
8.37 
7.01 
5.81 
8.  28 
79 
O' 

70286 

5.8 

4.085 

0 


Results  of  Search  for  best  RMS  % 


MATCHING  DISCREPANCY 

Field  DVe 

Point 

1  -.554154 

2  -.29179 

.  175305 

4  . 435654 

5  -.292571 

6  -.15499 

7  -.193712 

8  -.269948 

9  -.0802561 

10  -.0121288 

11  -.196139 

12  -.0660992 

13  .1 03865 

14  .0636177 

15  .117956 

16  -.081616 

RMS  =  .23948 
RMS  percent  =  3.37573 
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RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  6  W  =  4 

RESULTANT  =7.2111  at  33.6926  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  79 


RUN  No.  4  DATE  70286  Iteration  from  Exp. 79  Run  3 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

Number 

setting 

Field 

Point 

DVe 

1 

16.5 

1 

-.978442 

1 1 . 03 

-.445638 

O 

10.  93 

. 365099 

4 

12.  41 

4 

. 66099 

5 

13.  16 

5 

.084311' 

6 

2.05 

6 

.21119 

7 

S.  38 

7 

— . 406503 

a 

7.03 

8 

-.595732 

9 

5.  79 

9 

-.508798 

10 

8.39 

10 

-.444517 

1 1 

79 

11 

-. 143645 

12 

4 

12 

-.084854: 

13 

70286 

13 

. 155099 

14 

6 

14 

•  1 625c*o 

15 

4 

15 

.  087083' 

16 

0 

16 

-.098895' 

RMS  =  .422059 
RMS  percent  =  5. 3529 
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RESULTS  OF  ITERATION 
*********  *********** 


STREAM  VECTOR:  U  =  5.8  W  =  4.085 

RESULTANT  =  7.09417  at  35.16  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  79 


run  No.  4  date  70286  Results  of  Search  for  best  RMS  % 


CONTROL 

SETTINGS 

Control 

Number 

setting 

1 

16. 5 

2 

11.03 

10.93 

4 

12.41 

5 

13.  16 

6 

2.05 

7 

8.38 

8 

7.03 

9 

5.79 

10 

8.39 

11 

79 

12 

4 

13 

70286 

14 

5.8 

15 

4.085 

16 

0 

MATCHING  DISCREPANCY 


Field 

Point 

DVe 

1 

-1.33692 

n 

-.206444 

3 

. 40407 

4 

.62727 

5 

-6. 17526E 

6 

. 0573889 

7 

-. 136523 

8 

-.20582 

9 

-.0597374 

10 

.0598714 

11 

-. 188621 

12 

-. 10387 

1 3 

.112956 

14 

. 105059 

15 

.0697844 

16 

-. 140757 

RMS  =  .399456 
RMS  percent  =  5.63077 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  5.8  W  =  4.0e5 

RESULTANT  =  7.09417 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  79 


at  35. 16  degrees 


RUN  No. 


DATE 


70286 


CONTROL  SETTINGS 


Control 

K  Number 

k 

setting 

1 

15.75 

10.58 

1  3 

10.  16 

"  4 

12.26 

5 

12.91 

£  6 

7 

8.39 

8 

7.04 

5.79 

|  10 

8.28 

79 

».  12 

5 

ft  13 

100286 

a  14 

5.8 

15 

4.085 

W  16 

0 

Iteration  from  Exp. 79  Run  3 

MATCHING  DISCREPANCY 

Field  DVe 

Point 

1  -.947109 

2  -.273098 

3  .375127 

4  . 695357 

5  .0856213 

6  .149337 

7  -.0927312 

8  -.0443029 

9  .0674641 

10  . 197754 

11  -.180847 

12  -.15422 

13  .119786 

14  .0929612 

15  .0416573 

16  -.135764 

RMS  =  . 33348 
RMS  percent  =  4.70076 
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RESULTS  OF  ITERATION 
********** ********** 


STREAM  VECTOR:  U  =  5.7  W  =  3.?? 

fS  RESULTANT  —  6.95774  at  34.9946  degrees 

FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  79 


RUN  No. 

5  DATE 

100286 

Results  of 

Search  for  best 

CONTROL 

SETTINGS 

MATCHING  DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

F'oi  nt 

1 

15.75 

1 

-.645738 

4m 

10.58 

-.634112 

3 

10.  16 

3 

. 146173 

4 

12.26 

4 

. 39207 

5 

12.91 

5 

-.247964 

6 

2 .  *?2 

6 

-.209483 

7 

8.39 

7 

-.206181 

8 

7.04 

8 

-. 135774 

9 

5.79 

9 

3. 65514E- 

10 

8.28 

10 

. 168032 

11 

79 

11 

—.21 3385 

12 

5 

12 

-. 162347 

13 

100286 

13 

.0989419 

14 

5.7 

14 

.0742737 

15 

3.99 

15 

.0316274 

16 

0 

16 

-. 156922 

RMS  =  .285213 
RMS  percent  =  4.09922 


RESULTS  OF  ITERATION 
**************** 


STREAM  VECTOR:  U  =  5.8  W  =  4.085 

RESULTANT  =  7.09417  at  35.16  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  79 


run  No.  6  DATE  100286  Iteration  from  Exp. 79  Run  5 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

:ontrol 

lumber 

setting 

Field 

Point 

DVe 

1 

15.6 

1 

-.345092 

10.  48 

rj 

-.297454 

3 

9.86 

sJ 

. 105489 

4 

12.21 

4 

. 468233 

5 

12.88 

5 

-.0824493 

6 

2.33 

6 

. 0860336 

7 

8.61 

7 

B.74807E 

8 

7.  19 

8 

-.0444072 

9 

5 . 75 

9 

—3. 22543E 

10 

8.24 

10 

.0112419 

1 1 

79 

1 1 

-.234525 

12 

6 

12 

-. 102781 

13 

100286 

13 

. 183556 

14 

5.8 

14 

.0959734 

15 

4.085 

15 

.0403286 

16 

0 

16 

-. 106463 

RMS  =  . 189696 
RMS  percent  =  2.67397 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  5.9  W  =  4.085 

RESULTANT  =  7.09417  at  35.16  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  79 


run  No.  7  DATE  170286  Iteration  from  Exp. 79  Run  6 


CONTROL  SETTINGS  MATCHING  DISCREPANCY 


Control 

setting 

Fi  eld 

DVe 

Number 

Point 

1 

15.75 

1 

-.589862 

10.46 

-.372613 

3 

9.79 

3 

.252179 

4 

12.2 

4 

.617525 

5 

12.85 

5 

6.57074E 

6 

2.32 

6 

. 150168 

7 

8.57 

7 

-. 119603 

8 

7.27 

8 

-.258021 

9 

5.74 

9 

-.0435012 

10 

8.24 

10 

.0231954 

1 1 

79 

11 

-.236806 

12 

7 

12 

-. 16452 

13 

170286 

13 

. 152135 

14 

5.8 

14 

■  0Sc.852o. 

15 

4.085 

15 

.0162397 

16 

0 

16 

-. 129659 

RMS  =  .270186 

RMS  percent  =  3. 

80857 

RESULTS  OF  ITERATION 
************. *  ******* 


STREAM  VECTOR:  U  =  5.8  W  =  4.085 

RESULTANT  —  7.09417  at  35. IS  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  79 


RUN  No.  8  DATE  21028S  Repeat  of  Exp.79  Run  6 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

Point 

1 

15.6 

1 

-.391283 

*-> 

10.  48 

2 

-.580943 

-T 

W 

9.86 

T 

w 

.  1  3  7  /  55 

4 

12.21 

4 

. 483906 

5 

12.88 

5 

-. 126402 

S 

.  jo 

6 

-.0553834 

7 

8.61 

7 

-.0575841 

8 

7.  19 

8 

-. 196711 

9 

5.75 

9 

-.0397111 

10 

8.24 

10 

. 0566793 

11 

79 

11 

-.20895 

12 

8 

12 

-. 128327 

13 

210286 

13 

. 151052 

14 

5.8 

14 

. 0890202 

15 

4.085 

15 

.014661 

16 

0 

16 

-.08536 

RMS  =  .238217 
RMS  percent  =  3.35792 
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Figure  6,5 


RESULTS^ OF^ ITERATION 


STREAM  VECTOR:  U  =  1.8  W  =  1.563 

RESULTANT  =  2.38718  at  41.0625  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  85 


RUN  No.  1  DATE  if 0486 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

Number 

setting 

Field 

Point 

DVe 

1 

15.5 

1 

-.116501 

'I 

8. 17183 

o 

.544412 

T 

6. 00858 

O 

-. 12777 

4 

10.5035 

4 

-.231903 

5 

11.028 

5 

-.474919 

6 

3.86468 

6 

-.241853 

7 

8.79956 

7 

.251761 

8 

7.9893 

8 

. 497649 

9 

6. 19908 

9 

.502178 

10 

7.31773 

10 

.  ■-< 2^547 

1 1 

85 

11 

-.207143 

12 

1 

12 

-.0957222 

13 

280486 

13 

.0886521 

14 

1.8 

14 

.0107302 

15 

1.568 

15 

. 0304997 

16 

0 

16 

. 34282 

RMS  =  .307257 
RMS  percent  =  12.8711 


RESULTS  OF  ITERATION 
*4^***************** 


STREAM  VECTOR:  U  =  1.82  W  =  1.585 

RESULTANT  =  2.41343  at  41.055  degress 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  85 


run  No.  l  date  280486  Results  of  Search  for  best  RMS  % 


CONTROL  SETTINGS  MATCHING  DISCREPANCY 


Control  setting  Field  DVe 

Number  Point 


1 

15.5 

1 

-. 170165 

2 

8. 17188 

2 

.610006 

3 

6.00858 

3 

-.0853957 

4 

10.5035 

4 

-. 175139 

5 

1 1 . 028 

5 

-.412167 

6 

3.86468 

6 

-. 173966 

7 

8.79956 

7 

.271034 

8 

7.9893 

8 

.512004 

9 

6. 19908 

9 

.510974 

10 

7.31773 

10 

.325614 

11 

85 

1 1 

-.200773 

12 

1 

12 

-.0940777 

13 

280486 

13 

. 092824 

14 

1.82 

14 

.0146059 

15 

1.S85 

15 

.0324867 

16 

0 

16 

.347049 

RMS  =  .308363 
RMS  percent  =  12.777 


RESULTS  QF  ITERATION 
****#######*##**#*#* 


STREAM  VECTOR:  U  =  1 . 82  W  =  1.585 

RESULTANT  =  2.41343  at  41.055  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  85 


run  No.  2  DATE 300486  Iteration  from  Exp. 85  Run  1 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

Number 

setting 

Field 

Point 

DVe 

1 

11.4 

1 

“ . 435637 

8.27 

2 

.510571 

3 

*j .  iik j 

O 

.0235796 

4 

10.61 

4 

. 0968493 

5 

11.68 

5 

-.0259511 

6 

3.24 

6 

.0981959 

7 

8.83 

7 

. 250926 

8 

7.96 

8 

.394471 

9 

6 .  54 

9 

.384797 

10 

7.45 

10 

.280164 

1  1 

85 

1 1 

-. 15901 

12 

2 

12 

-. 131575 

13 

300486 

13 

.0776174 

14 

1.82 

14 

5. 92517E- 

15 

1 . 585 

15 

-.0125584 

16 

0 

16 

.290259 

RMS  =  .256041 
RMS  percent  =  10.609 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  1.75  W  =  1.524 

RESULTANT  =  2.32058  at  41.0542  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  85 


RUN  NO.  2  DATE  300486  Results  of  Search  for  best  RMS  % 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

Point 

1 

11.4 

1 

-.242045 

2 

8.27 

2 

. 275229 

T 

w 

5.25 

3 

-. 12771 

4 

10.61 

•  4 

-. 105101 

5 

11.68 

5 

-.249041 

6 

3.24 

6 

-. 14279 

7 

8  •  8  ^ 

7 

. 180489 

8 

7.96 

8 

. 340792 

9 

6.54 

9 

. 35055 

10 

7.45 

10 

. 269546 

11 

85 

11 

-. 181424 

12 

2 

12 

-. 137314 

13 

300486 

13 

.0630184 

14 

1.75 

14 

-.0075192 

15 

1.524 

15 

-.0195295 

16 

0 

16 

.275457 

RMS  =  .212193 
RMS  percent  =  9. 14398 


RESULTS  OF  ITERATION 


STREAM  VECTOR:  U  =  1.75  W  =  1.524 

RESULTANT  =  2.32058  at  41.0542  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  85 


run  No.  3  date  300486  Iteration  from  Exp. 85  Run  2  k  =  .20 


CONTROL  SETTINGS  MATCHING  DISCREPANCY 


Control 

Number 

setting 

Field 

Point 

DVe 

1 

10 

1 

nnnT 

2 

8.26944 

2 

.592106 

3 

6.07842 

3 

.0647479 

4 

10.8797 

4 

.22101 

3 

11.9547 

5 

. 137909 

6 

2.83911 

6 

. 206503 

7 

8.85002 

7 

. 255825 

8 

7.95207 

8 

.328259 

9 

6.92137 

9 

.319278 

10 

7.57455 

10 

. 322066 

1 1 

85 

1 1 

-. 181326 

12 

3 

12 

-.0917266 

13 

300486 

13 

.0611748 

14 

1.75 

14 

. 0183705 

15 

1.524 

15 

. 0602798 

16 

0 

16 

. 180844 

RMS  =  .250694 
RMS  percent  =  10.8031 
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RESULTS  OF  ITERATION 
****** ************** 


STREAM  VECTOR:  U  =  1.65  W  =  1.44 


RESULTANT  =  2.19  at  41.1151  degrees 


FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  85 


RUN  No. 


DATE  20586 


Iteration  from  Exp. 85  Run  3 


CONTROL  SETTINGS 


Control 

Number 


setting 


10 

0 . 

4.54 

11.06 

12.26 

2.29 

8.82 

7.98 

6.91 

7.44 
85 

4 

20586 

1.65 

1.44 
0 


MATCHING  DISCREPANCY 

Field  DVe 

Point 


-.0497232 
6 . 30045E— 03 

_  riTr'^'T/i 

-. 128973 
-. 171625 
-.244519 
.0631197 
. 19496 
. 182288 
. 209005 
— . 176553 
-. 102482 
. 0265805 
-.0293676 
6. 90059E— 03 
.26868 


RMS  =  . 157592 
RMS  percent  =  7. 19599 


RESULTS  OF  ITERATION 
##*###*#******#**#*# 


8 


3 


STREAM  VECTOR:  U  *  1.69  U  *  1.475 

RESULTANT  =  2.24315  at  41.1169  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  S5 


run  No.  4  date  205S6  Results  of  Search  for  best  RMS  % 


CONTROL  SETTINGS  MATCHING  DISCREPANCY 


Control 

Number 

V 

setting 

Field 

Point 

DVe 

1 

10 

1 

-.  16085 

>,  2 

8.33 

.  141284 

3 

4.54 

O 

-.  145523 

4 

11.06 

4 

-.013274 

_  5 

12.26 

5 

.  04  j84-->5 

||  6 

2.2? 

6 

-. 10651? 

"  7 

8.82 

7 

.  108629 

8 

7.98 

8 

. 225933 

o  9 

6.91 

9 

.202159 

v-  10 

7.44 

10 

.215368 

11 

85 

1 1 

-. 163735 

m  12 

4 

12 

-.0992041 

1 3 

20586 

13 

. 0349225 

1 4 

1 . 69 

14 

-.0216956 

15 

1.475 

15 

.0108855 

\;  1 6 

is 

0 

16 

RMS  »  .147016 

.27713? 

t 

RMS  percent  =  6.55 

399 

£ 


151 

N 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  1.69  M  =  1.475 

RESULTANT  =  2.24315  at  41.1169  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  85 


run  No.  5  DATE  20586  Iteration  from  Exp. 85  Run  4 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

Number 

setting 

Field 

Point 

DVe 

1 

10 

1 

-. 1240G6 

2 

8 . 33 1 86 

2 

. 0882299 

5.57215 

— p  ' 

•J 

-.0755432 

4 

11.4224 

4 

. 0444258 

5 

12.2534 

5 

.0403712 

6 

1 . 96293 

6 

4.43399E 

7 

8.88627 

7 

.0872111 

8 

7.96552 

8 

. 128843 

9 

7.22249 

9 

. 189288 

10 

7.58023 

10 

. 131057 

1 1 

85 

11 

-. 186712 

12 

5 

12 

-.0791567 

13 

20586 

13 

.0450741 

14 

1 . 69 

14 

-1.6449SE 

15 

1.475 

15 

. 0237484 

16 

0 

16 

.252978 

RMS  =  . 116577 
RMS  percent  =  5. 


19704 


RESULTS  OF  ITERATION 
t***tt******tt**t*t* 


STREAM  VECTOR:  U  =  1 . 67  W  =  1.4575 

RESULTANT  =  2.21658  at  41.116  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  85 


run  No.  5  date  20586  Results  of  Search  for  best  RMS  % 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

Number 

setting 

Field 

Point 

DVe 

1 

10 

1 

-.0688394 

2 

8.38186 

r> 

a. 

.0210545 

-> 

w 

5.57215 

-. 118735 

4 

1 1 . 4224 

4 

-.0132353 

5 

12.2534 

5 

-.0233293 

6 

1 . 96293 

6 

-.0643806 

7 

8.88627 

7 

.0671198 

8 

7.96552 

8 

.113544 

9 

7.22249 

9 

. 179543 

10 

7.58023 

10 

. 128061 

11 

85 

11 

-. 193114 

12 

5 

12 

-.0807967 

13 

20586 

13 

. 0409029 

14 

1.67 

14 

—4. 00707E 

15 

1 . 4575 

15 

.0217569 

16 

0 

16 

.248748 

RMS  =  .111253 
RMS  percent  =  5.01915 


RESULTS  OF  ITERATION 
It******************* 


STREAM  VECTOR:  U  =  1.74  W  =  1.44 


RESULTANT  =  2.25858  at  39.6136  degrees 


FLAPS:  Position  4,  50  p.s.i 

EXPERIMENT  No.  85 


RUN  No. 


DATE 


20586  , 


CONTROL  SETTINGS 


Control 

Number 


setting 


10 

8.38186 

5.57215 

1 1 . 4224 

12.2534 

1 . 96293 

8.88627 

7.96552 

7.22249 

7.58023 

85 

5 

20586 

1.74 

1.44 

0 


Results  of  Search  for  best  RMS  % 
(continued)  Note  Stream  Angle. 

MATCHING  DISCREPANCY 

Field  DVe 

F'oi  nt 

1  .0135629 

2  -.019598 

3  -.110114 

4  .0242319 

5  .0341787 

6  .0147121 

7  -5.11 166E-03 

8  2. 74072E-03 

9  . 0482082 

10  -.0232134 

11  -.176472 

12  -.0742648 

13  .0556326 

14  .0152082 

15  .0279352 

16  .26342 

RMS  =  .0893433 
RMS  percent  =  3.95572 


mmw 
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RESULTS  OF  ITERATION 
*tt*t*t*t*******t**t 


STREAM  VECTOR:  U  =  1.67  W  =  1.4 s 

RESULT  ANT  =  2. CISCO  <s  *  J;  .  lc-i'  deo^ees 

FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  S5 


RUN  No.  6  DATE  50586  Iteration  from  Exp. 85  Run  5 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

Point 

1 

10 

1 

. 0793736 

8.53 

-. 0457516 

■»' 

5.57 

-. 101274 

4 

11.74 

4 

. 0422947 

5 

12.  35 

5 

. 11598 

6 

1 . 56 

6 

.0732114 

7 

8.85 

7 

-3.32529E 

8 

7.93 

8 

. 0867657 

9 

7.  46 

9 

.0792918 

10 

7.69 

10 

. 0430581 

11 

85 

1 1 

-.21767 

12 

6 

12 

—  •  07d>/ 648 

13 

50586 

13 

.0361148 

14 

1.67 

14 

.018466 

15 

1.46 

15 

.032166 

16 

0 

16 

. 387o9 / 

RMS  =  . 1276 
RMS  percent  =  5.75236 


RESULTS  OF  ITERATION 
******************** 


STREAM  VECTOR:  U  =  1.67  W  =  1.46 


RESULTANT  =  2.21S22  at  41.1647  degrees 


FLAPS:  Position  4,  50  p.s.i 

EXPERIMENT  No.  85 


RUN  No. 


DATE  70586 


CONTROL  SETTINGS 


Control 

Number 


setting 


10 

8.45 

«j  .  * — 1 6 

1 1 . 56 
12.3 
1.76 
8.S7 
7.95 
7.34 
7.64 
85 
8 

705S6 

1.67 

1.46 

0 


Iteration  from  Exp. 85  Run  5 


MATCHING  DISCREPANCY 

Field  DVe 

Point 

I  -.0506968 

2-  -.0760766 

3  -.0411339 

4  . 130697 

5  .168979 

6  -7.73967E— O: 

7  .13605 

8  . 0609092 

9  .0716079 

10  .046942 

II  -.172297 

12  -.14297 

13  7.81066E— OZ 

14  .0385913 

1 5  . 0680 1 ? 1 

16  .509341 

RMS  =  . 158317 
RMS  percent  =  7.13712 


RESULTS  OF  ITERATION 
*X*****t*tt*%+*t**** 


STREAM  VECTOR:  U  =  1.67  W  =  1.46 

RESULTANT  —  2.21822  at  41.1647  degrees 
FLAPS:  Position  4,  50  p.s.i. 

EXPERIMENT  No.  85 


run  No.  7  date  5<35e6  Iteration  from  Exp. 85  Run  6  k 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

Point 

1 

10 

1 

-.010317 

8.75 

-.0897504 

4.86 

•J* 

-.0716949 

4 

11.94 

4 

•  u6o  1 45-o 

5 

12.  46 

5 

. 133764 

6 

1 . 0*0 

6 

-. 0186475 

7 

8.84 

7 

. 163413 

8 

7.91 

S 

-.0296856 

9 

7.  55 

9 

. 0277688 

10 

7.79 

10 

-. 0609595 

1 1 

85 

11 

-. 140828 

12 

7 

12 

-. 14654 

13 

50586 

13 

.0101 033 

14 

1.67 

14 

. 0499643 

15 

1.46 

15 

.0731423 

16 

0 

16 

.541955 

RMS  =  .16019 
RMS  percent  =  7.22155 
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Figure  7  continued 


REPEATABILITY  OF  DATA  : 

RESULTS  OF  REPEATED 

RUNS 

************************************************* 

EXPERIMENT  No.  68 

MODEL  ANGLE 

-1 

3  deg. 

FLAPS:  Position  4,  50  p.s.i. 

CONTROL 

SETTINGS 

Control 

Number 

Setting 

1 

16.9 

10.78 

O 

11.05 

4 

12.24 

5 

12.87 

6 

2.4 

7 

7.81 

8 

7.  18 

9 

6 

10 

8.  12 

1 1 

68 

12 

5 

13 

141085 

14 

6 

15 

4 

COMPARISON 

OF 

MEASURED 

VELOCITY  COMPONENTS 

First  Run  No. 

5  ,  Date 

141085 

Second  Run  No.  8  , 

Date  161085 

Total 

Vt 

I 

Total 

Ve 

1st  Run 

2nd  Run 

1st  Run 

2nd  Run 

4.23919 

4.21201 

1 

6.44238 

6.42154 

4.01946 

3.94985 

2 

5. 80013 

5.60667 

6.45808 

6.48029 

3 

4. 12322 

4.05132 

7. 09294 

7.08124 

4 

3.26854 

3. 15614 

7.04003 

•  7. 12825 

5 

2. 87443 

2.82268 

6.22331 

6. 27334 

6 

3. 33719 

3 . 33 1 3 

5 .  4— *  4 

5.45318 

7 

4.00201 

4.00532 

5. 96974 

5.97476 

8 

3. 6578 

3. 66596 

5.84734 

5.84811 

9 

3. 39676 

3.29016 

6.41814 

6.5486 

10 

3. 67782 

3 . 55034 

7.00129 

6.99625 

11 

6.85492 

6.86298 

6.20825 

6.38871 

12 

6. 20107 

6.2938 

5.40469 

5. 18764 

13 

5.0409 

5 . 07738 

5.85429 

5.80224 

14 

5.72838 

5.69822 

5 . 70306 

5.67647 

15 

5.73208 

5. 67806 

6.25031 

6.25839 

16 

6.25176 

6.20425 

Average  Absolute  Vt  Difference  =  .05 

71029 

Average  Absolute  Ve  Difference  =  .0606714 
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REPEATABILITY  OF  DATA  :  RESULTS  OF  REPEATED  RUNS 
******* ******************************************* 


EXPERIMENT  No.  71  &  79 

MODEL  ANGLE  -13  deg.  FLAPS:  Position  4,  50  p.s.i. 

CONTROL  SETTINGS 

Control 

Number  Setting 


1 

16.3 

2 

10.  18 

3 

10. 1 

4 

12.  15 

5 

12.52 

6 

2.25 

7 

7.66 

8 

6.81 

9 

5.94 

10 

8.  13 

11 

71 

12 

7 

13 

41 185 

14 

6 

15 

4 

COMPARISON  OF 

MEASURED 

VELOCITY  COMPONENTS: 

First  Run  No. 

7  ,  Date 

41185 

Second  Run  No. 

1  ,  Date 

(Experiment  71 ) 

(Experiment  79) 

Total 

Vt 

I 

Total 

Ve 

1st  Run 

2nU  Run 

1st  Run 

2nd  Run 

4.28228 

4.09406 

1 

6.4629 

6.42253 

3.99854 

3.85721 

2 

5.97534 

5.70471 

6.53492 

6.5314 

3 

3.84941 

3.79966 

7. 17134 

7. 12269 

4 

3.09096 

3. 16774 

6.93882 

7.09076 

5 

2.71163 

2.80407 

5.98967 

6.09619 

6 

2.98414 

2.7574 

5.72478 

5.42862 

7 

3.99896 

3.99021 

5.84648 

5.89751 

8 

3.49858 

3.69706 

5.9892 

5.72944 

9 

3. 45947 

3.45324 

6.57948 

6.72709 

10 

3.68209 

3.8341  1 

6.93922 

6.94898 

1 1 

6.86336 

6.95904 

6.4404 

6.57778 

12 

6.41552 

6.49035 

4.90309 

5.20067 

13 

5.00919 

5.25744 

5.88147 

5.69875 

14 

5.81445 

5.62685 

5.64628 

5.66514 

15 

5.71598 

5. 6703 

6. 18861 

6. 18011 

’’  16 

6. 11359 

6.08999 

Average  absolute  Vt  difference  =  .  12S097 
Average  absolute  Ve  difference  *  .112365 
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REPEATABILITY  OF  DATA  :  RESULTS  OF  REPEATED  RUNS 
************************************************** 


EXPERIMENT  No. 
MODEL  ANGLE  -13  deg. 

CONTROL  SETTINGS 


79 


FLAPS:  Position  4,  50  p.s.i 


Control 

Number 


Setti no 


4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 


15 . 6 
10.48 
9.86 
12.21 
12.88 
2.33 
8.61 
7.  19 
5.75 
8.24 
79 
6 

100286 

5.8 

4.085 


COMPARISON  OF  MEASURED  VELOCITY  COMPONENTS: 


First  Run  No.  6  ,  Date  100286 


Second  Run  No.  8  ,  Date  210286 


Total  Vt 


Total  Ve 


1st  Run 


2nd  Run 


1st  Run 


2nd  Run 


4.31279 

4.36856 

1 

6.28777 

6.32178 

4. 18014 

4. 18853 

2 

5.81149 

5.8611 

6.2954 

6.29914 

3 

4.01818 

4.02078 

6.94463 

6.94917 

4 

3. 1583 

3. 1465 

6.94712 

6.95626 

5 

2.80563 

2.82171 

5.77942 

5.80797 

6 

2.37799 

2.41674 

5.28446 

5.24913 

7 

3.828 

3. 78655 

5.78452 

5.81904 

8 

3.55441 

3.61096 

5.64362 

5.65182 

9 

3.33741 

3.32761 

6.6468 

6.59196 

10 

3 . 80663 

3.71753 

6.88781 

6.78102 

11 

6.81279 

6.70137 

6.23427 

6.07107 

12 

6.22176 

6.11424 

4.92517 

4. 79086 

13 

4.97767 

4.90358 

5.83399 

5.72269 

14 

5.74627 

5.67487 

5.57275 

5.54071 

15 

5.63728 

5.6194 

6.0517 

5.95527 

16 

6 . 05365 

5 . 95435 

Average  absolute  Vt  difference  —  .0554438 
Average  absolute  Ve  difference  =  .0519622 
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RESULTS  OF  RUNS  AT  CONSTANT  CONTROL  SETTINGS 
*******  *  *  ****  **  ********  *******  ******** 


&  VARYING  ANGLE  OF  ATTACK 
*  *******  *****  **  ****** 


STREAM  VECTOR:  U  =  4.7  W  =  2.94 


RESULTANT  =  5.54379  at  32  degrees 


FLAPS:  Position  4,  50  p.s.i. 


ANGLE  OF  ATTACK  =  6  deg. 


EXPERIMENT  No. 


RUN  No.  1 


DATE 


30685 


CONTROL  SETTINGS 


MATCHING  DISCREPANCY 


Control 

Number 


setting 


Field 

Point 


9.  1 

9.25 

10.41 

11.5 

11.85 

2.79 

7.81 

6.78 

5.76 

7.94 
53 

1 

30685 
4.  7 

2.94 
0 


-.301896 
.0107004 
-.111419 
-. 189839 
-.356 
. 284407 
-. 154613 
-.0463801 
-.0963289 
.0346619 
-. 176361 
-. 103039 
.0946593 
. 104255 
.0826186 
. 0586202 


RMS  =  . 16876 
RMS  percent  ■  3.04413 


6ft 


l 


ft 


RESULTS  OF  RUNS 

AT  CONSTANT 

CONTROL  SETTINGS 

*******  **  **** 

**  ******** 

*******  ******** 

&  VARYING  ANGLE  OF 

ATTACK 

*  *******  *****  ** 

****** 

STREAM  VECTOR: 

U  =  4.65  W 

s  2.9 

RESULTANT  = 

=  5.48019 

at  32  degrees 

FLAPS: 

Position  4,  50 

p. s. i . 

ANGLE  OF 

ATTACK  =  11  deg. 

EXPERIMENT  No. 

52 

RUN 

No.  1  DATE  30685 

CONTROL  SETTINGS 

MATCHING  DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

Point 

1 

9.  1 

1 

-.568961 

2 

9.25 

.0377719 

3 

10.41 

3 

6. 76S54E-03 

4 

11.5 

4 

-.0755723 

5 

11.85 

5 

-.259827 

6 

2.79 

6 

. 295727 

7 

7.81 

7 

-.226295 

8 

6.78 

8 

-. 133467 

9 

5.76 

9 

-. 174384 

10 

7.94 

10 

— 1 . 03 1 66E— 03 

1 1 

52 

1 1 

-. 153294 

12 

1 

12 

-.0765105 

13 

30685 

13 

. 180485 

14 

4.65 

14 

. 0904809 

IS 

2.9 

15 

. 0800023 

16 

0 

16 

. 0237688 

RMS  =  .2 

03467 

RMS  percent 

=  3.71276 

m 
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RESULTS  OF  RUNS  AT  CONSTANT 
*******  **  ****  **  ******** 

?<  VARYING  ANGLE  OF 
*  *******  *****  ** 

STREAM  VECTOR:  U  »  4.7  W 

CONTROL  SETTINGS 
*******  ******** 

ATTACK 

****** 

=  2.94 

RESULTANT  > 

=  5.54379  at  32 

degrees 

FLAPS: 

Position  4,  50  p.s.i. 

ANGLE  OF  ATTACK 

-  16  deg. 

EXPERIMENT  No.  51 

RUN  No.  1  DATE  310585 

CONTROL  SETTINGS 

MATCHING  DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

1 

9.  1 

Poi  nt 

1 

— . 353755 

afc. 

9.25 

2 

. 12608 

T 

10.41 

3 

-.079341 

4 

11.5 

4 

-. 141969 

5 

11.85 

5 

-.291962 

6 

2.79 

6 

.314784 

7 

7.81 

7 

-. 185051 

8 

6.78 

8 

— . 135731 

9 

5.76 

9 

-.246977 

10 

7.94 

10 

. 0675498 

11 

51 

1 1 

-. 157955 

12 

1 

12 

-.0761504 

13 

3 1 0585 

13 

. 246903 

14 

5 

14 

. 101027 

15 

2.98 

15 

.0743992 

16 

0 

16 

. 0858678 

RMS  =  . 191192 

RMS  percent  *  3.44877 

RESULTS  OF  RUNS  AT  CONSTANT  CONTROL  SETTINGS 
*******  **  ****  **  ********  *******  ******** 


&  VARYING  ANGLE  OF  ATTACK 
*  *******  *****  **  ****** 


FLAPS: 


STREAM  VECTOR: 

Position  4,  50 

EXPERIMENT  No. 


U  *  4.67  W  =  2.918 

RESULTANT  =  5.50669  at  32  degrees 
p.s.i.  ANGLE  OF  ATTACK  =21  deg. 

50 


RUN  No.  2  DATE  310585 


CONTROL  SETTINGS  MATCHING  DISCREPANCY 


Control 

Number 

setting 

Field 

Point 

DVe 

1 

9 

1 

-.443813 

2 

9.25 

2 

-.0728233 

3 

10.41 

3 

-.044828 

4 

11.5 

4 

-. 127905 

5 

11.85 

5 

-.265883 

6 

2.79 

6 

. 189212 

7 

7.81 

7 

-.237221 

8 

6.78 

8 

-.21007 

9 

5.76 

9 

-.357874 

10 

7.94 

10 

-.0560596 

11 

50 

11 

-. 15325 

12 

2 

12 

-. 142676 

13 

3 1 0585 

13 

.245655 

14 

4.67 

14 

. 121542 

15 

2.918 

15 

. 126694 

16 

0 

16 

. 159922 

RMS  =  .212035 
RMS  percent  =  3.85051 


RESULTS  OF  RUNS  AT  CONSTANT  CONTROL  SETTINGS 
*******  **  ****  **  ********  *******  ******** 


S<  VARYING  ANGLE  OF  ATTACK 
*  *******  *****  **  ****** 

STREAM  VECTOR;  U  =  6.28  w  =  4.42 

RESULTANT  =  7.67951  at  35.14  degrees 
FLAPS:  Position  4,  50  p.s.i.  ANGLE  OF  ATTACK  =  4.14  deg. 

EXPERIMENT  No.  91 


RUN  No.  4  DATE  140786 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

Point 

1 

1  kJ  .  6 

1 

-.301606 

10.48 

—9. 27832E 

O 

9.86 

3 

.316667 

4 

12.21 

4 

. 74759 

5 

12.88 

5 

.326153 

6 

2.33 

6 

-1. 11889 

7 

8.61 

7 

. 132387 

8 

7.  19 

8 

-.415298 

9 

5.  75 

9 

-.018726 

10 

8.24 

10 

-.205758 

1 1 

91 

11 

-.0876777 

12 

4 

12 

-.0438776 

13 

1 40786 

13 

. 102822 

14 

6 

14 

. 14044 

15 

4.42 

15 

. 0386385 

16 

0 

16 

-. 110281 

RMS  =  .386875 
RMS  percent  =  5.03776 


RESULTS  OF  RUNS  AT  CONSTANT 
*******  **  ****  *x  ******** 

&  VARYING  ANGLE  OF 
*  *******  *****  ** 

STREAM  VECTOR:  U  =  6.28  W 

RESULTANT  = 

CONTROL  SETTINGS 
*******  ******** 

ATTACK 

****** 

=  4.47 

=  7.70839  at  35. 

45  degrees 

FLAPS: 

Position  4,  50  p.s.i. 

ANGLE  OF  ATTACK 

=  10.45  deg 

EXPERIMENT  No.  91 

RUN  No.  3  DATE  140786 

CONTROL  SETTINGS 

MATCHING  DISCREPANCY 

Control 

setting 

Field 

DVe 

Number 

1 

15.6 

Point 

1 

-.362643 

2 

10.48 

2 

-.0339967 

— * 

o 

9.86 

3 

.  328435 

4 

12.21 

4 

. 766749 

5 

12.88 

5 

. 154448 

6 

2  .  o  o 

6 

-1.09042 

7 

8.61 

7 

. 182402 

a 

7.  19 

8 

— . 315458 

9 

5.75 

9 

-. 193219 

10 

8.24 

10 

-.288987 

i  l 

91 

1 1 

-.0935399 

12 

3 

12 

.0170164 

13 

140786 

13 

.  182854 

14 

6.2.8 

14 

. 1 37 1 33 

15 

4.4 

15 

.0621131 

16 

0 

16 

-.116777 

RMS  =  .385171 

RMS  percent  =  4.99677 
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RESULTS  OF  RUNS  AT  CONSTANT  CONTROL  SETTINGS 
*******  **  ****  **  ********  *******  ******** 


?<  VARYING  ANGLE  OF  ATTACK 
*  *******  *****  *.*  ****** 


STREAM  VECTOR:  U  =  6.16  W  =  4.39 

RESULTANT  =  7.56424  at  35. 
FLAPS:  Position  4,  50  p.s.i.  ANGLE  OF  ATTACK 

EXPERIMENT  No.  91 


5  degrees 
=  16.5  deg 


RUN  No. 


DATE 


1 107S6 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

Number 

setting 

Field 

F'oi  nt 

DVe 

1 

15.6 

1 

-. 564145 

r* 

10.  43 

-. 171625 

9.86 

. 171623 

4 

12.21 

4 

. 532479 

5 

12.88 

5 

0304471 

6 

n  rr-r 
■  -J-J 

6 

-.678878 

7 

8.61 

7 

. 0884154 

8 

7.  19 

8 

-.315222 

9 

5 . 75 

9 

142459 

10 

8.24 

10 

-.262817 

1 1 

91 

1 1 

-. 165648 

12 

r> 

12 

-. 047596 

13 

1 10736 

13 

. 325178 

14 

6.  16 

14 

. 113085 

15 

4.39 

15 

. 114835 

16 

0 

16 

-. 172657 

RMS  =  .307234 
RMS  percent  =  4.06167 
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RESULTS  OF  RUNS  AT  CONSTANT  CONTROL  SETTINGS 
*******  **  ****  **  ********  *******  ******** 


S<  VARYING  ANGLE  OF  ATTACK 
*  *******  *****  **  ****** 


FLAPS: 


STREAM  VECTOR: 

Position  4,  50 

EXPERIMENT  No. 


U  =  6  W  =  4.42 

RESULTANT  =  7.45227  at  36.33  degrees 
p.s.i.  ANGLE  OF  ATTACK  =  23.38  deg. 

91 


RUN  No.  1  DATE  110786 


CONTROL 

SETTINGS 

MATCHING 

DISCREPANCY 

Control 

Number 

setting 

Field 

Point 

DVe 

1 

15.6 

1 

-.250883 

2 

10.48 

2 

-.336079 

3 

9.86 

■j 

. 1 6634 1 

4 

12.21 

4 

.5817 

5 

12.88 

5 

-5. 9849  IE 

6 

2.33 

6 

-. 475758 

7 

8.61 

7 

.256896 

8 

7.  19 

8 

-.20276 

9 

5.75 

9 

-.0634727 

10 

8.24 

10 

-.278764 

11 

91 

1 1 

-.20282 

12 

1 

12 

.0542531 

13 

110786 

13 

. 294798 

14 

6 

14 

■  .113819 

15 

4.42 

15 

. 0827493 

16 

0 

16 

-. 111046 

RMS  =  .264608 
RMS  percent  =  3.5507 


EFFECTS  OF  MATCHING  ERROR 
*******  «*  ********  ***** 


Experiment  Run  No.  6 
No.  34 

RMS  Error 


DVe’s  (input) 


Speed  5.7 
4.027. 


Model  Angle 
-11 


-.558 

-.258 

-.346 

.291 

.  135 

.047 

.  179 

.075 

-.033 

-.271 

-.  118 

.  107 

-.  107 

-.246 

-.  124 

.  107 

P 

Vx 

Vz 

P 

1 

-.0557379 

-.020418 

1 

2 

-.0537508 

-.020287 

2 

~T 

N-* 

-.0471657 

-.0212467 

3 

4 

-.0819244 

. 155086 

4 

5 

-.0405932 

-.0644855 

5 

6 

-.0520729 

-.0814909 

6 

7 

-.0335878 

-.0971938 

7 

I 


I 
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EFFECTS  OF  MATCHING  ERROR 
*******  **  ********  ***** 


Experi ment 
No.  62 


Run  No.  8  Speed  7.22 

RMS  Error  =  3.57% 


Model  Angle 

-.18 


DVe's  (input) 


384 

-.096 

.  187 

028 

.  158 

.089 

072 

-.31 

.044 

483 

— •  00c. 

-.  16 

Vs 

P 

0946182 

3.66S17E-03 

1 

0916128 

2.04809E-03 

2 

0820689 

— 5. 68907E-03 

T 

w 

140832 

. 126375 

4 

0875125 

-.0404226 

5 

138242 

-.0399161 

6 

105478 

— .  07S0c.c. 

7 

w 

m 

[v  1 

I*.** 

V,\ 


EFFECTS  OF  MATCH I NB  ERROR 
*******  **  ********  ***** 


175 


E:< peri  men t 

Run  No.  6 

Speed  7.36 

Model 

No.  76 

RMS  Error  = 

4.57% 

— 

DVe's  (input) 

-.211 

.  12 

.517 

.039 

-.274 

.304 

.  17 

-.255 

-.61 

.361 

-.043 

-.329 

-.407 

-.  135 

— .  63 

.215 

-.  335 

-.032 

-.309 

.207 

.017 

-.018 

-.086 

.094 

_  /  -rn 

.  148 

.221 

.586 

-.777 

-.  128 

-.  147 

P 

V>; 

Vs 

P 

1 

-.0595572 

.0186814 

1 

O 

-.0727993 

. 0258054 

O 

si> 

-.0892194 

.0290155 

O 

4 

-. 100622 

.0273918 

4 

5 

-.0991806 

.0186742 

5 

6 

-. 179233 

. 0392906 

6 

7 

-.0183678 

8. 43228E— 04 

7 

8 

-.0659455 

-.0287526 

8 

9 

.0254231 

. 0544307 

9 

EFFECTS  of  matching  error 
*******  **  ********  ***** 


periment 

Run  No .  6 

Speed  7.21 

Model 

No.  79 

RMS  Error  = 

2.67% 

DVe's  (input) 

345 

-.082 

-.003 

.  184 

297 

.086 

.011 

.096 

105 

.009 

— .  235 

.04 

468 

-.044 

-.  103 

-.  106 

Angle 


P 


V:< 


Vs 


P 


1 


4 

5 

6 
7 


-.0695064 

-.067178 

-.0594161 

-.117406 

-.0539953 

-.0993245 

-.0430545 


.0395217  1 
.0366136  2 
. 0255524  3 
. 0385572  4 
.0282769  5 
.0562356  6 
. 0203902  7 


